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XXII. BAND SECHSTES (SCHLUSS-)HEFT 1954 


Strémung durch Schlitz- und Lochblenden bei kleinen Reynolds-Zahlen*. 
Von W. Wuest. 


1, Einleitung. Kleine Reynoldssche Zahlen treten auBer bei der Strémung in sehr zihen 
Medien auch bei sehr geringen Gasdichten auf. Bei sehr weitgehender Verdiinnung des Gases 
ist schlieBlich die freie Weglinge groB gegen den mittleren Wandabstand, d.h. die Strémung 
hat molekularen Charakter (molekulare Gasstrémung). Im Grenzfall unendlich grofer freier 
Weglange ergeben sich besonders einfache Verhaltnisse. Fiir diesen Grenzfall sind daher ver- 
schiedene Strémungsprobleme, darunter auch der Durchflu8 durch Lochblenden, schon sehr 
frihzeitig gelést worden (Knudsen 1909). Bei endlicher freier Weglange treten Abweichungen 
infolge von Zweierstéfen zwischen den Molekiilen auf. Dieses Problem ist neuerdings von 
Hiby und Pahl * untersucht worden, auch fiir den Fall der Lochblende. Bei weiterer Steigerung 
der Gasdichte nimmt schlieBlich die Strémung molaren Charakter an, d.h. die Gesetze der 
Kontinuumsphysik sind anwendbar. Die Reynoldssche Zahl kann dabei aber immer noch sehr 
klein sein. Die Strémung durch lange Rohre und Spaltkandle geht allmahlich in die Poiseuille- 
strémung iiber, das entsprechende Gegenstiick der Strémung durch Schlitz- und Lochblenden 
ist dagegen bisher nur wenig untersucht worden. Zu erwahnen ist hier lediglich eine experi- 
mentelle Arbeit von F. C. Johannsen *, der die laminare Strémung durch Lochblenden bei sehr 
kleinen Reynoldszahlen untersucht hat. Im folgenden wird der Versuch gemacht, diese Zu- 
sammenhange auch theoretisch klarzustellen. Da das ebene Problem der Schlitzblende leichter 
beherrscht werden kann als das rotationssymmetrische Problem der Lochblende, beginnen wir 
die Untersuchung mit dem erstgenannten. Um die Probleme einer Berechnung zuganglich zu 
machen, werden fernerhin folgende Beschrankungen hinsichtlich der geometrischen Ab- 
messungen und der Strémungsform vorausgesetzt: 

1. Spaltweite bzw. Durchmesser der Offnung klein gegen die Abmessungen des Rohres, in 
dem sich die Blende befindet. 

2. Wandstarke der Blende klein gegen die Offnung. 

3. Strémung inkompressibel. Soweit es sich also um eine Gasstrémung handelt, soll die 
Druckdifferenz vor und hinter der Blende klein gegen den absoluten Druck sein. 


2. Die Schlitzblende; Grenzfall sehr kleiner Reynoldszahl. Der Grenzfall sehr kleiner 
Reynoldszahl ist dadurch gekennzeichnet, da8 nur Druck- und Reibungskrafte beriicksichtigt 
werden, wahrend Tragheitskrafte vernachlassigt werden. Die Navier-Stokesschen Gleichungen 
nehmen dann folgende Form an: 


grad p = rotrotb , diy b= 0". (1.1) 


Fir das folgende brauchen wir diese Differentialgleichungen in krummlinigen isometrischen 
Koordinaten: Wenn wir mit H(é,7) den ortlich veranderlichen Mafistabsfaktor des ebenen, 
orthogonalen und isometrischen Koordinatensystems bezeichnen, so gehen mit 


OE hee ea a ; 
== 0 ta | Bn = 3E ag ) 
die Gleichungen (1.1) ttber in 
One dw op =e do OH ws) O(H un) es () les 
3 am’ an BE 0& on i (ce) 
woraus man ableiten kanu 
Ap ==); Ay =O. (1.4) 


* Vorgetragen auf der Gamm-Tagung Aachen 1953. 
1 J. W. Hiby u. M. Pohl, Z. Naturforschung 7a, (1952) S. 533, 542. 
2 F.G. Johannsen, Proc. Roy. Soc. London Ser. A 126, (1930) S. 231. 
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Fiir das vorliegende Problem empfiehlt es sich, zu Kegelschnittkoordinaten iiberzugehen, die 
aus konfokalen Ellipsen und Hyperbeln bestehen (Abb. 1). Der Zusammenhang zwischen der 
(x, y)-Strémungsebene und den neuen Koordinaten £, 7 ist dabei gegeben durch 


x = co) Ecosy, y=cOmésinn, 
und der Ma®stabsfaktor H ist 


Hee “+ (Gof 2 & “ee 2 ny : (1.5) 


Da die Strémungsebene fir beliebige Werte €, 7 un- 
endlichfach iiberdeckt wird, miissen wir noch eine Ver- 
einbarung iiber eine eindeutige Zuordnung machen. 
Zunichst soll 7 nur Werte zwischen 0 und z durch- 
laufen. Ferner wollen wir annehmen, da wir beim 
Durchgang durch den Schlitz in ein anderes Riemann- 
sches Blatt der (x, y)-Ebene gelangen und zwar so, 
daB negativen Werten von y negative € entsprechen 
(vgl. Abb.1). Der Potentialgleichung (4) geniigen be- 
liebige Funktionen der komplexen Variablen € = & 
+in. Es zeigt sich, daB auch die Randbedingungen 
u, =u, = 0 fir 7 = 0 und 7 erfulle werden, wenn 
man folgenden Lésungsansatz macht: 


(ts. 2 ce 
p=4C,R (890) =" gsr —eaay 18) 


Es folgt dann namlich 


Abb. 1. Kegelschnitt-Koordinaten fiir die Berechnung sin 2 1] 
der Strémung durch eine Schlitzblende. W C, Ee ? (1.7) 
Coj 2 — cos 2H 


so da man zusammen mit der Kontinuitatsgleichung folgendes Gleichungssystem fiir die Ge- 
schwindigkeiten erhalt: 


(Hus) OAH uy) _ c? O(Hus) , O(H un) 
a aia C, 5 Sin 2h aE a on Diese (yin (1.8) 
Die Lésung lautet offenbar unter Beachtung der Randbedingungen 
er FE: me Cie 
Up On us = gees a (cos 27 — 1). (1.9) 
Die Konstante C, kann man noch auf die durchflieBende Flissigkeitsmenge zuriickfihren: 
0= [ud = Cre Geog 2 —1)d mse OLR 
eer 4 fi 1h Uae oa (1.10) 
6 
Die Geschwindigkeitsverteilung ist dann 
Be 
Me dx Va 2E—cos 27 (I cos 2 n) (1.11) 
mit der Spaltweite d = 2c. 
und der Druckabfall Ap ist gegeben durch 
Asis sce 32 
AD Pee Pie ee AC a (1.12) 


Fir grofen Abstand r vom Spalt kann man naherungsweise (pj 2 & = 2 r/c? setzen und 
im Vergleich damit cos 2 1 vernachlassigt, so da8 (1.11) iibergeht in 


2Q 
Us =—~ (1 — cos 27m). (1.13) 
Im Spalt selbst ist dagegen 


sin 7] » (1.14) 
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d.h. also elliptisch verteilt, da 


sinn = yi i : 

c 
3. Schlitzblende bei maBigen Reynoldszahlen. Wenn man sich auf reine Radialstrémungen 
beschrankt, d. h. also die unmittelbare Umgebung des Spaltes selbst ausschlieBt, ist nach dem 
Vorgang von G. B, Jeffery? und G. Hamel? eine exakte Lésung der Navier-Stokesschen Differen- 
tialgleichungen auch bei Beriicksichtigung der Tragheitsglieder miglich. Bezeichnet man mit 


wu die Radialgeschwindigkeit, so nehmen die Differentialgleichungen in Polarkoordinaten r, # 
folgende Form an: : | 


du ep. (ou 1 O?u 1 du u 
Paar (eee Lege (a 5 r2 Op TanOn r : ey) 
1 op- 2 Ou ‘ 
ord rap’ (2.2) 
O(u r) P 
a =O (2.3) 
Aus der letzten Gleichung folgt dann 
ur = F(d) (2.4) 
und damit aus (2.2) 
2 Bi 
p = 5 F(d) + h(r), (2.5) 


wobei h zunachst noch eine beliebige Funktion von r ist. Setzt man diese Ergebnisse in (2.1) 
ein: 


Fv LAF — a pe : (2.6) 
so ist offenbar zu schlieBen, daB 
u C 
eae f e 
hr) = oe + konst. 


Die Randhbedingungen lauten 
=) fir 0 = s. 
n/2 
Ferner gilt fiir die durchstrémende Menge Q = [{ F dé. Da es sich um eine flachenhafte Menge 


handelt, kann man |Q|/y = R, also gleich einer fairies Zahl setzen. Mit f, = F/Q geht 
(2.6) tuber in 
fit4fit RA =G. (2.7) 
Im Grenziibergang R— 0 erhalt man offenbar bei Beachtung der Randbedingungen 
u =" (1 + cos 29). (2.8) 


Dies ist aber, abgesehen von einer anderen Wahl des Nullpunktes von @ gegeniiber 7 identisch 
mit (1.13). 
Fir den Druckverlauf gilt 


P= ine cos 2 # +- const. (2.9) 
Dabei ist also Q auf der Einstrémseite negativ, auf der Ausstrémseite positiv; die Konstante 
ist auf beiden Seiten verschieden, jedoch unbestimmt. Vergleichen wir den Druckverlauf nach 
Formel (2.9) mit demjenigen nach der strengen Formel (1.6), und zwar auf der Mittellinie, d. h. 
fir & =0 bzw. 7 =2/2. Unter Beriicksichtigung von y = c Gin é auf der Mittellinie geht 


(1.6) itiber in 


eee ee (2.10) 
me Vit aly 


1G, B. , Phil. Mag. Ser. 6, 29 (1915) S. 455. 
a C UnChIe EER aati Vereinigung 25 (1916) S. 34, engl. Ubersetzung NACA TM 1342 (1953). 
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Durch asymptotische Entwicklung fiir groBe y ergibt sich 


Paes ce Ee (2.11) 


Dh Mee 
Man erkennt daraus, da8 zumindest auf der Mittellinie (2.9) den Anfang einer asymptotischen 
Entwicklung der strengen Lésung darstellt. Durch Vergleich zwischen (2.9) und (2.11) kann 
auch die Konstante bestimmt werden. 

In Abb. 2 ist der Druckverlauf in der Mittelachse nach der strengen Formel und nach der 
asymptotischen Lésung aufgetragen. Man entnimmt aus der Auftragung, daB die asympto- 
tische Lésung etwa bis zu einer Entfernung vom Schlitz gleich der Schlitzweite giiltig ist. 


Ale 


Abb, 2. Druckverlauf in der Mittelebene (7 = 2/2) bei der Strémung durch eine Schlitzblende fir R = 0. 
Die ausgezogene Kurve entspricht der exakten Lésung (2.10), die gestrichelte Kurve der asymptotischen Naherungsformel (2.11). 


Um eine Lésung fiir allgemeine R zu erhalten, wollen wir zu f = Rf, ibergehen, womit 
(2.7) folgende Form annimmt: 


fl +4ftf—GQ=0. (2:0 


In der Mitte (d. h. fir % = 0) ist aus Symmetriegriinden f’ = 0. und f=f,,. Durch weitere 


Integration von (2.12) erhalt man 
if 


o= df (2.13) 
/D +2 C,f — 2/3 f? —4f? 


7/2 i 
mit der Bedingung f fdd = +R. 
—n/2 
Wenn fy =f, = fe die Wurzeln des Nenners sind, kénnen wir (2,13) in der Form schreiben 


= if 
2 df 
i 2.14 


Durch Vergleich mit (2.13) gilt ferner 


fo +h +fe=— 6. (2.15) 


Als Losung kommt nur der Wertbereich f, >f =f, in Betracht, wobei f, > 0 sein muB, um die 
Bedingungen an der Wand erfiillen zu kénnen. Eine der Wurzeln muf dabei gleich dem Wert fi, 
eS der Mitte sein und zwar gilt im Fall des Einstrémens f,, = f,, fiir den Fall des Ausstrémens 
m = fo. 
Wie schon Hamel gezeigt hat, ist ein reines Ausstrémen in eine Halbebene (Winkel 180°) 
nur fiir den Grenzfall R = 0 méglich, fiir den bereits eine exakte Lisung angegeben wurde. 
Fiir endliche Reynoldszahl sind also auf der Ausstrémseite nur radiale Strémungen mit teil- 
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weisem Kinstrémen an der Wand méglich. In Abb.3 ist dieser Lésungstyp mit denim folgenden 
gebrauchten Bezeichnungen schematisch dargestellt. 


Durch die Substitution /f, — f = t erhalt man das unvollstandige elliptische Integral erster 
Gattung 


Hierbei ist 


Ui oon (ere ie 

C=fo —fp=2fo tft. 
e wird also hier voriibergehend in 
einer anderen Bedeutung gebraucht. 


Dabei ist die Bedingung (vgl. Abb. 3) 
v, a D, =F 


! 
- % Q iE 30° 60° 50° 
Abb. 3. Geschwindigkeitsverteilung bei einer radialen Abb. 4. Geschwindigkeitsverteilung auf der Einstrémseite (gestrichelt) 
Strémung. und Ausstrémseite (ausgezogen) einer Schlitzblende bei verschiedenen 
Reynoldszahlen. 


zu erfiillen oder nach Einfiihrung des unvollstandigen elliptischen Integrals F(a, py) in der bei 
Jahnke-Emde, Funktionstafeln, tabellierten Form 


F(a, $1) + Fag) =< =, (2.17) 


wobei 


; a : c if 
we the >, =n/2 ae eG) a 
i= are sin(*), Dy == 7/2. Po are sin(£ S 


Die Gleichung (2.16) muB durch schrittweise Naherung aufgelést werden. Zu einem ge- 
gebenen f, wird so das zugehdérige f, bestimmt. Nach dessen Ermittlung kénnen leicht Zwi- 
schenwerte berechnet werden, wobei zu einem gegebenen f der Winkelabstand ? (Abb. 3) 
folgendermaBen bestimmt wird: 


6 
mit 
As Oe 
gp = arc sin eo ; 
Auf diese Weise sind die ausgezogenen Kurven in Abb. 4 ermittelt worden. Die Kurve f, = 0 


2/3 
entspricht der Verteilung cos* 9. Durch Integration { fd erhalt man die Reynoldszahl R. 
—n/2 
Wenn man R iber f, auftragt (Abb. 5), erkennt man, dafs nur Reynoldszahlen bis R = 2,85 mit 
radialen Stromlinien in geniigendem Abstand vom Schlitz méglich sind. 
Auf der Einstrémseite geht die Berechnung ganz entsprechend vor sich. Nur beschrankt 
sich hier der Wertevorrat von f auf den mit #, bezeichneten Bereich. Entsprechend gilt hier die 
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v7 Bedingung 0 =2/2. In Abb. 4 sind 
gestrichelt einige Geschwindigkeits- 
verteilungen auf der Einstrémseite 
aufgetragen. 

Man erkennt dabei, da® die Ge- 
schwindigkeitsverteilung auf der Ein- 
strémseite sich mit zunehmender Rey- 
noldsscher Zahl einer ,,Grenzschicht- 
lésung** nahert, bei der der Geschwin- 
digkeitsabfall sich auf den wandnahen 
Teil beschrankt. Dies wird noch deut- 
licher, wenn man die Druckverteilung 
auftragt, die nach (2.5) gegeben ist 
durch: 


‘ (e= zits + (F = ze konst. 


% 


wobei 


i ne. oan 
? a Y ioe C= Sgt + Side - fo fo) - (2.20) 


Abb. 5. Abhangigkeit der Reynoldsschen Zahl von der gréBten Ausstrém- 
geschwindigkei it bei Annabme von radialen Stromlinien. 


Die auf diese Weise fiir die Ein- 
strémseite berechneten Druckvertei- 
lungen (Abb. 6) lassen erkennen, daf fiir R = 0 der Druck entsprechend der Funktion cos 2 0 
verteilt ist und sich mit wachsender Reynoldsscher Zahl immer mehr einer gleichformigen 
Verteilung nahert, wie es von der Grenzschichttheorie vorausgesetzt wird. 


fees Die Grenzschichtlésung kann 
10 erie: Se - tbrigens in einfacher Weise durch 
| ‘s Spezialisierung der Lésung (2.14) 
Bi, gewonnen werden. Wenn sich der 
Geschwindigkeitsabfall nur auf die 
Wandschicht beschrankt, mu im 
.< mittleren Teil f’ =f’? =0 sein. Dies 
rs fihrt aber auf die Bedingung f, =f, 
Ee | ; : und somit wegen (2.15) auf 
a | jo = inet: 
Die Lésung (2.14) nimmt damit fol- 
Jee Se an: 


| 6/2 = =|; 2 

, la zalie\ sates : (ff /—CP+ of) 

(2.21) 
Dieses Integral kann aber durch die 
Substitution y = j= (27 6. 
leicht berechnet werden. Wenn 
man die Integrationskonstante so 
-95t # ieee 2 ; ‘ bestimmt, daB am Rande, d. h. fiir 
%, =2/2f =0 wird, erhalt man 
die Lésung 


ie (9 — 8) 


( 49 (thct 8 
vi ae tg PALES 
405 ae ae oo ine ECT aes ei) 2. fica 8 


Abb. 6. Druckverlauf auf der Einstrémseite einer Schlitzblende in Abhingigkeit —_— Gy tq : 2f, + 6 es i (2 22) 
von der Re-Zahl 3f, +6 heh , 
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oder nach f aufgelést 


f= (3 f, + 6) Tg? IV/- a +42) (J — B) + Ar tg ee a 4 (Qipao) a .(2.23) 


3f; ae 
Mit wachsender Reynoldsscher Zahl nimmt f, gro®e negative Werte an, so daB die Zahl 2 als 
Summand daneben vernachlassigt werden kann. Setzt man noch f, = — A/y, so daB auBerhalb 
der Grenzschicht die Geschwindigkeit durch u = — A/r gegeben ist, so erhalt man die erstmals 
_von K. Polhausen? abgeleitete Grenzschichtlisung : 
: Ae a A | y 
f=- 3 — Bg? IVs (% — 0d) 4 1,146] +2 aS (2.24) 


Man kann (2.19) benutzen, um den Druckverlust auf der Einstrémseite angendhert zu be- 
rechnen. Fiir R—>» 0 ist der Druckverlust auf der Einstrémseite gleich dem halben Gesamt- 
druckverlust, also nach 
(1.12) gegeben durch 

16 Q ¥ i 7 


APey == alt « | we 


Vou 


Entsprechend kann man YA 
fir R + 0 formal setzen : 
4, 160 [2 Ge 
pPe= 7 al 2R u u 10} f = Ww a 
= Ap.y G(R). (2.25) VA Ke 
Die Funktion 5 4 
G(R) = Ap./APeo 
ist in Abb. 7 aufgetragen. 
Sie nahert sich fiir groBe R_ 
asymptotisch der Geraden fit 
cf z Q 
G,(R) ~~ 39 Rw~ 39 a > 7 a dee gil Ge + 
woraus sich umgekehrt (Jarg) r 
ergibt a 2 2 a0. 60 100 200 500 


2 
Aperso0o™~ > = 0. (2.26) Abb. 7. Drucksprung bei der Schlitzblende und Lochblende [Formeln (2.25); (2.29); (3.18)] 


N 


Dies entspricht dem sich aus der Bernoullischen Gleichung ohne Beriicksichtigung der Strahl- 
kontraktion ergebenden Wert. 

Wesentlich schwierigere Verhaltnisse liegen auf der Ausstrémseite vor. Die vorhergehende 
Untersuchung ergab lediglich unter welchen Bedingungen als asymptotische Lésungen Radial- 
strémungen moéglich sind, ohne daB aber gesichert ist, da} solche sich iberhaupt einstellen. 

Die Jaminare Diffusor-Strémung ist wohl erstmalig von H. Blasius? untersucht worden und 
zwar fir einen exponentiell erweiterten Diffusor. Die an einem gleichartigen Diffusor durch- 
gefiihrten Messungen von G. N. Patterson® ergaben allerdings erheblich niedrigere Werte fiir 
die Durchflu8menge, bei der die Strémung abzureien beginnt. Neuerdings ist das gleiche 
Problem fir den keilférmigen Diffusor von M. Abramowitsch* und S. M. Targ® mit einem 
héheren Grad der Naherung untersucht worden, Wenn auch die Rechnung sich in diesen 
beiden Fallen auf verhaltnismaBig kleine Diffusorwinkel bezieht, so kénnen doch die Ergeb- 
nisse wenigstens als grobe Naherung auch fiir den Fall der Strémung hinter einer Schlitzblende 
angesehen werden. Fiir den hier besonders interessierenden Druckabfall erhalt Targ die Be- 
ziehung 

1 1 04 


= ee pe 0 
o U2 2 1—TqVa R/Vo R R 


« steht hier als Naherung fiir sin « und ist demnach im vorliegenden Fall durch sina = 1 zu 


Piso Pom (2.27) 


1 Polhausen, Z. Angew. Math. Mech. 1 (1921) S, 252. 
2 Blasius, Z. Math. Mech. Phys 58 (1910) S. 225. 
3 Patterson, Canad, J. Research 12 (1935) S. 676. 
4 Abramowitsch, J. Math. Phys. 28. (1949). 5.1. 


5 Targ, Osnovnye zadaé teorii laminarnyi teéenii (Hauptprobleme der Theorie d. laminaren Stré- 


mungen.) Moskau 1951. 
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ersetzen. Die Reynoldssche Zahl R bezieht sich auf die mittlere Geschwindigkeit im Schlitz 
und die Schlitzbreite und ist demnach R = Q/y. Um den Anschlu8 an den Grenzfall R = 0 zu 


finden, fiir den eine strenge Lésung vorliegt, setzen wir 


16 
Po — Po = et G(R), (2.28) 


wobei 
I 1 1 2, 


=— —|] - ——— 
He 2 1—3g/R/VR RB 


Diese Funktion ist in Abb. 8 aufgetragen. Fir kleine Werte von R liefert die Formel jeden- 
falls zu kleine Werte, weil G,(0) = 1 sein mu8. Fiir R = 8,3 wird G, = 0, d.h. auf der Aus- 
strémseite fallt der Druck nicht mehr ab. Die Zahigkeitsspannungen iberwiegen nicht mehr 

gegeniiber den Tragheitskraften. Bereits bei 

G(R) geringer Uberschreitung dieses Wertes, namlich 
08\— fir R = 9,87, lost sich nach Targ die Strémung 
in groBer Entfernung vom Schlitz ab und bei 
weiter wachsender Reynoldsscher Zah] wandert 
die Ablésestelle immer naiher an den Schlitz. 

Fir den Gesamtdruckabfall beim Durch- 


G2(R) 


06 


o4 strémen der Schlitzblende kann man folgen- 
den Ansatz machen: 
Be 

a2 Ap == Su G(R), (2.29) 


R wobeiG, = 4(6,+6,) ist. Nach der strengen 
2  Lésung ist dabei G,(0) =1 wahrend fiir groBe 
R entsprechend der Kirchhoffschen Theorie 


-02 der AusfluBzahl' folgende asymptotische Be- 
Abb.8. Drucksprung auf der Ausstrémseite einer Schlitzblende ziehung gilt : 
nach Targ. Der Kreis kennzeichnet beginnende Strémungs- Tat D2 
aaa : G,( RB) ~ 22 Rese Re ee 
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Fir mafBige Werte von R sind die im vorhergehenden nach der Formel von Targ errechneten 
Werte in Abb. 7 eingetragen. Vergleichbare Messungen liegen nicht vor. Doch lat ein Ver- 
gleich mit den experimentellen Ergebnissen bei der Lochblende vermuten, da G, héher liegen 
muB und etwa den strichpunktierten Verlauf hat. 

Fir den ausgebildeten ebenen laminaren Strahl ist die Geschwindigkeitsverteilung von 
W.G. Bickley? berechnet worden. Die Stromlinien sind dabei gekriimmt, Druck und Impuls 
bleiben langs des Strahles konstant. Auf den Druckabfall beim Durchstrémen des Schlitzes 
kénnen hieraus also keine Folgerungen gezogen werden. Auch ist der Umschlag des laminaren 


Strahles zur Turbulenz, der von C. C. Lin? und Mitarbeitern untersucht wurde, offenbar ohne 
Einflu8 auf den Druckabfall. 


4, Die Lochblende; Grenzfall sehr kleiner Reynoldszahl. Zur Beschreibung der rotations- 
symmetrischen Strémung durch eine Lochblende gehen wir von dem rechtwinkligen (x, y, z)- 
Koordinatensystem zu einem krummlinigen, orthogonalen (£, 7, ¢)-Koordinatensystem iiber, 
wohbei zwischen beiden folgenden Relationen bestehen mégen: 


w=f(En)eosS, y=s(é,n), «= f(E,7) sing. (3.1) 


Es ist also angenommen, da die y-Achse die Symmetrieachse bildet. Die értlich verander- 
lichen Mafstabsfaktoren der (€, 7, €)-Koordinaten sind dann 


i Vfi + 8} 7 Hi, = Vf3 + 85 Hy=f, (3.2) 
Im Grenzfall R->0, d.h. bei Vernachlassigung der Tragheitsglieder, nehmen die Navier-Stokes- 


' Vgl. W. Kaufmann, Angewandte Hydromechanik I. Bd. Berlin 1934. 
> Bickley, Phil. Mag. Ser. 7, 23 (1937) S. 727. 
weeded. Foote u. C, C. Lin, Quart. appl. Math. 8, (1950) S. 265. P. Chiarulli, Stability of two-dimen- 
sional velocity distributions of the half-jet type, Thesis Brown University 1949. 
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schen Differentialgleichungen bei Rotationssymmetrie folgende Form an: 


i op = Uu iL ‘ 0 ree 0( H, Un) ; o( Hy Ws) A 
H, 0& H, H; 0n|H, H, 0g an > 

1 ap _,,_1_ 2 | Hy (A Hyus)  (Hyus)\] Aue H, Hy) | (ue H,H,) (3.3) 
A, on H, H, 0g Jak fale on on ae + — an —(0., 


Wir untersuchen zunichst, ob auBerhalb der nichsten Umgebung des Loches radiale Stré- 
mungen miglich sind und nehmen daher zunichst ein riumliches Polarkoordinatensystem an. 
Dabei ist in (3.1) f(§, ) = sin und g(£, 4) = € cos 7 zu setzen. Man findet damit 


1: ok ee Hee. H, =€&siny . (3.4) 
Setzt man noch wu, = 0 voraus, so gehen die beiden ersten Gleichungen (3.3) tiber in 
Op 1 ON/ae Que 
bE SS ® sin n on (sin 1) fi) ° (3.5a) 
Op us 
an = OE an . (3.5b) 


Aus der Kontinuitatsgleichung findet man ferner 
ue €? sin ny = F,(y) sin n , (3.6) 


wobei also (7) eine noch zu bestimmende Funktion von 7 ist. Durch Einsetzen von (3.6) in 
(3.5b) findet man 
Op 2u 


an i 3 F(7) 


und durch Integration 
9) 
2 bb 


x Py(n) + h(E) ; (3.7) 


ieee 
h(€) bestimmt man zu h(é) ao +C,. 


Einsetzen in (3.5a) liefert schlieBlich folgende Differentialgleichung fiir F'\(7): 


6F,+6C+ Fi + cotgy Fi =0. (3.8) 
Bei Beriicksichtigung der Randbedingungen F’, = 0 fiir 7 = + 2/2 lautet die Lésung 
r= — 3 Ccos? 7 . (3.9) 
Die Konstante C bestimmt man aus der DurchfluBmenge 
a[2 
Q = f 2m & sin ug dy = Dep Os 
0 
Damit wird also die Geschwindigkeitsverteilung 
b= ae . cos? (3.10) 
und der Druck 
p= is (1 — 3 cos? 7) + konst. (3a) 


Die vorstehende Untersuchung hat gezeigt, daB die Stromlinien auch im Fall der Lochblende 
bei sehr kleinen Reynoldsschen Zahlen radial verlaufen. Die Konstante in Formel (3.11) ist 
auf beiden Seiten der Lochblende verschieden, doch liefert diese nur in einer gewissen Entfer- 
nung vom Loch giltige Theorie keine Angaben iiber den Drucksprung. Es ist naheliegend, 
analog zum frither behandelten ebenen Fall die Stromlinienflachen als Rotations- Hyperboloide 
anzunehmen. Wir gehen also zu einem Koordinatensystem tiber, das durch Rotation der Kegel- 
schnittkoordinaten in Abb. 1 um die y-Achse entsteht. Es ergeben sich somit folgende Zu- 
ordnungen: 


x =cojEsinycosl, y =c Ging cosy, 2 = cp) ésinn sing , (3:12) 
rie. Gy hs ey" gf ey H, =cGo| é sinn. (3.13) 


(Abweichend von Abb. | entspricht 7 = 0 der Mittelachse.) 
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Setzt man u, = 0 voraus, so mub us H, A; = F,(y) sein, wobei wir F (7) daraus be- 
stimmen, da} shee groBe € wieder die alte Lésung (3.10) erhalten miissen. Die Gleichungen 


gehen tiber in 


if / 7 
Lee ay la — cot Fa)). | 


ag Hoey (3.14) 
op J 0 1 , ~y 
af = are fe 7 (Hi (F; — cotg 7 ry).| 
Nach Einsetzen von Fy — cotgy F, = — 3 Q/m cos 7 sin? 7 und Integration erhdlt man aus 
der zweiten Gleichung (3.14) fiir den Druck p 
; 3 Q 12 ; Gin ae es By i [is 
p= & [ware testy tHE): ee 


Die Funktion h(€) miBte man nun so bestimmen, daB die erste der Gleichung (3.14) erfullt 
wird. Es zeigt sich, daB dies aber nicht streng mioglich ist. Unsere Annahme, daf die Strom- 
linienflachen Rotationshyperboloide sind, ist 
also nur angenahert richtig. Wir kénnen aber 
eine angenaherte Lésung erhalten, wenn wir 
die Funktion h(&) so wahlen, daB der Fehler 
moglichst klein ist. Wenn das Fehlerquadrat 
der Abweichung zwischen der rechten Seite 


| 


40 
Y 
a 
15 
G5 
0 7 @ J = == -0,50 
Abb. 9. Hilfsfunktion h(é) zur Berechnung des Druckverlaufes Abb. 10. Druckverlauf in der Mittelachse einer Lochblende bei 
bei der Lochblende vgl. Formel (3.15)]. R= 0 (gestricheltasymptotische Lésung mit radialen Stromlinien). 


der ersten Gleichung (3.14) und der nach & differentiierten Gleichung (3.15), itber 7 integriert, 

ein Minimum werden soll, erhalt man 

7 Gof? § Coj92 E+ 7 (3 Goj2E + Sin? 4) Gop2E + 
HONS Coj F (Coj* DE + Gof? DE + 3/16) (ag 

Die Funktionen h’(é) und h(é) sind in Abb. 9 dargestellt. Da h(oo) — h(— oo) = 2,10 dees 

erhalt man fiir den Drucksprung bei der Lochblende 


— 24 Qu ; tg Q Lu 
(AW Oe met 2,10 = 50,4 re (3.17) 
Experimentell ist von Johannsen‘ fiir den Grenzfall R = 0 gefunden worden 
eet Qu 
Ap = 66,8 = j (3.17a) 


Dieser Unterschied ist vermutlich nicht nur auf den Naherungscharakter der Rechnung zu- 
ruckzufihren, sondern wird auch damit zusammenhingen, das bei den Versuchen eine Blende 
mit endlicher Dicke und konisch zugescharfter Offnung verwendet wurde. Hierdurch ergibt 
sich offenbar eine Erhéhung des DurchfluBwiderstandes gegeniiber dem der Rechnung zu- 
grundegelegten Idealfall. 

Wie die Auftragung in Abb. 10 zeigt, vollzieht sich der Druckabfall bei der Lochblende auf 
einer kirzeren Strecke als bei der Schlitzblende. Fiir den Fall endlicher kleiner Reynoldsscher 


Ny 
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Zahlen kénnen bei der Lochblende noch weniger Aussagen als im ebenen Fall gemacht werden. 
Radiale Strémungen in Analogie zu den Jefferey-Hamelschen Strémungen sind im rotations- 
symmetrischen Fall nicht méglich. Eine Lisung durch Reihenentwicklungen ist fiir das Dif- 
fusorproblem von NV. A. Sljozkin! angegeben worden. Diese Reihen konvergieren jedoch nur 
fiir kleine Diffusorwinkel geniigend schnell und kénnen also auf das vorliegende Problem der 
Lochblende nicht angewandt werden. Andererseits liegen aber fiir die Lochblende die experi- 
mentellen Ergebnisse von Johannsen vor, nach denen man fiir den Druckabfall setzen kann 


Ap = 06,8 24 GARye (3.18) 


wobei G, in Abb. 7 emgetragen ist. Diese Werte sind an einer in ein Rohr eingebauten Blende 
erhalten worden und gelten fiir den Fall, da® die Blendenéffnung klein gegen den Lochdurch- 
messer ist. 

Nach den Beobachtungen von Johannsen ist der Stromlinienverlauf bis zu Reynoldsschen 
Zahlen R = uy d/y = 10 zu beiden Seiten der Blende symmetrisch, bei gréBerer Reynoldszahl 
beginnt die Strémung in Ubereinstimmung mit den theoretischen Betrachtungen im ebenen 
Fall sich von der Wand abzulésen. Bei R = 30 zeigt der Strahl schon deutlich abgeliste Struk- 
tur, ist aber weiterhin laminar. Der laminare Strahl beginnt etwa bei R = 150 instabil zu 
werden. Bei Uberschreitung dieses Wertes bilden sich am Strahlrand Wirbelrippeln, die bei 
weiterer VergréBerung der Reynoldszahl in ausgepragte Wirbelringe itibergehen in Analogie zum 
ebenen Fall der Kdrmanschen WirbelstraBbe. Die Wirbelringe rufen stromabwarts Turbulenz 
hervor. Uber den ZahigkeitseinfluB bei gréBeren Reynoldsschen Zahlen liegen Untersuchungen 
von M. Hansen? vor. 


5. Zusammenfassung. Bei der Strémung durch Schlitz- und Lochblenden ist fiir sehr 
kleine Reynoldsche Zahlen (R < 2, wobei R auf die Schlitzbreite bzw. den Lochdurchmesser 
bhezogen wird) der Druckfall proportional der durchflieBenden Menge. Erst bei gréBeren Rey- 
noldsschen Zahlen bildet sich eine quadratische Abhangigkeit heraus. Fiir das ebenen Problem 
und den Grenzfall R = 0 wird eine exakte Liésung der Navier-Stokesschen Differentialglei- 
chungen mit auf der Ein- und Ausstrémseite symmetrischen Stromlinien hergeleitet und eine 
Formel fiir den Druckabfall aufgestellt. Fir endliche Reynoldssche Zahlen existieren nur Nahe- 
rungslésungen, die lediglich eine Abschaétzung des Druckverlaufes in Abhangigkeit von der 
Reynoldsschen Zahl liefern. Auf der Ausstrémseite lést sich die Strémung etwa bei R > x? ab 
und bildet einen laminaren Strahl. Ahnliche Verhaltnisse gelten auch fiir die Lochblende. Auch 
hier wird fiir den Fall R = 0 eine Liésung abgeleitet, die aber nur naherungsweise giiltig ist. 
Der rechnerisch ermittelte Wert fiir den Druckabfall bei sehr kleiner Reynoldsscher Zahl liegt 
etwa 25% unter dem von Johannsen (1930) gemessenen Wert. Der Unterschied ist vermutlich 
auf die endliche Dicke und konische Offnung der Lochblende bei den genannten Versuchen zu- 
rickzufihren. 


Aus dem Max-Planck-Institut fiir Str6mungsforschung, Gottingen. 
(Eingegangen am 19, November 1953). 
Anschrift des Verfassers: Dr. W. Wuest, Gottingen, Bunsenstrafe 10. 


1 Sljozkin, Math. Sbornik 42 (1935) 5.43, 
2 Hansen. Uber das AusfluBproblem, VDI-Forschungshett 438 (1949/50). 
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Zur Berechnung ebener und drehsymmetrischer Grenzschichten 
mit kontinuierlicher Absaugung. 


Von K. Wieghardt. 


1. Einleitung. Der Strémungswiderstand eines Kérpers kann bekanntlich durch Absaugen 
der Grenzschicht oft wesentlich verringert werden. Einerseits kann dadurch etwa bei einem 
Kérper mit stumpfer Hinterkante die Ablésung der Grenzschicht verzégert und damit der 
Druckwiderstand verkleinert werden, wie in den bekannten Strémungsaufnahmen von L. 
Prandtl aus dem Jahr 19041. Andererseits kann aber auch der Reibungswiderstand strom- 
linienformiger Kérper mit vernachlassigbarem Druckwiderstand noch dadurch gesenkt werden, 
daB die laminare Grenzschicht durch Absaugen stabilisiert und somit nicht turbulent wird; 
denn besonders bei hohen Reynolds-Zahlen ist die laminare Reibung viel kleiner als die der 
turbulenten Reibungsschicht bei unbeeinfluBter Strémung. 

Bereits vor etwa zehn Jahren war aus theoretischen Arbeiten von H. Schlichting?® und 
J. Pretsch* zu ersehen, daB auf diese Weise — auch bei Beriicksichtigung der Absaugeleistung — 
der Widerstand um die Halfte und mehr verringert werden kénnte. Trotzdem sind nur wenige 
Versuche bekannt geworden, diese groBen Méglichkeiten technisch zu verwirklichen. Nach den 
ersten Versuchen schien es so, als ob sich bei groBen Reynolds-Zahlen (iiber 6- 10°) diese 
Methode der Widerstandssenkung wegen der aufzubringenden Absaugeleistung nicht lohnte>. 
Spatere Windkanalversuche ©? an einem Tragfliigel aus gesinterter Bronze bestatigten aber 
die optimistische, theoretische Voraussage doch zum Teil. Die Ergebnisse fiir diesen Fligel 
(bei Re = 6- 10% und 20- 10°) sind in den Abb. 3 und 4 miteingetragen. Danach ist eine 
Widerstandssenkung von 50% und mehr auch bei Einbeziehung der Absaugleistung durchaus 
zu verwirklichen; allerdings la8t der unregelmaBige Gang der Widerstandsbeiwerte mit der 
Reynolds-Zahl noch auf Ungenauigkeiten der komplizierten Versuche schlieBen. 

Dies alles bezieht sich auf kontinuierliche Absaugung durch eine porése Oberflache. 
Experimentell eingehender untersucht worden ist die Absaugung aus mehreren Einzelschlitzen 
in mehr oder weniger weiten Abstanden voneinander. Dieser theoretisch wesentlich kompli- 
ziertere Fall der diskontinuierlichen Absaugegrenzschicht ist erst vor kurzem von W. Wuest® 
der Rechnung zuganglich gemacht worden. Es scheint, daB bei dieser technisch leichter zu 
verwirklichenden diskontinuierlichen Absaugung nicht so hohe Widerstandeinsparungen er- 
wartet werden kénnen, wie bei der theoretisch einfacheren aber praktisch schwierigeren kon- 
tinuierlichen Absaugung. 

Im folgenden wird ein Naherungsverfahren zur Berechnung der stationaren, inkompres- 
siblen laminaren Grenzschicht bei kontinuierlicher Absaugung dargestellt, das sowohl fiir ebene 
wie fir drehsymmetrische Kérper gilt. Ein Verfahren fir ebene Kérper hat H. Schlichting® 
bereits 1945 entwickelt. Die von ihm speziell zur Darstellung der kontinuierlichen Absaugung 
aufgestellte Klasse von Grenzschichtprofilen wird fiir das vorliegende Verfahren itbernommen. 
Methodisch ist dieses eine Erweiterung eines Verfahrens fiir ebene Grenzschichten ohne Ab- 
saugung von A. Walz®, insofern als auch hier Impuls- und Energiesatz von einer einpara- 
metrigen Profilklasse erfillt werden. Es ist vor allem zur direkten Berechnung der wirtschaft- 
lich giinstigsten Absaugung geeignet. 


1 L. Prandtl, Uber Fliissigkeitsbewegung bei sehr kleiner Reibung. Verh. d. II. Intern. Math. Kongr, 
Heidelberg 1904. 

° HA. Schlichting, Luftfahrtforschung 19 (1942), S. 293. 

3 H. Schlichting, Ing.-Arch. 16 (1948), S. 201. 

* J. Pretsch, Die Leistungsersparnis durch Grenzschichtbeeinflussung beim Schleppen einer ebenen 
Platte. UM 3048 (1943). 

: A. E. von Doenhoff und L. K. Loftin, Journ. aero Sci. 16 (1949), S. 729, 
_ § D.L. Burrows, A. L. Braslow and N. Tetervin, Experimental and theoretical studies of area suction 
for the control of the laminar boundary layer on a porous bronze NACA 64 A 010 airfoil, NACA Techn, 
Note 1905 (1949). 

’ A. L. Braslow and F. Visconti, Further experimental studies of area suction for the control of the 


laminar boundary layer on a porous bronze NACA 64 A 010 airfoil, NACA Techn. Note 2112 (1950). 
8 W. Wuest, Ing.-Arch. 21 (1953), S. 90. 


9 A. Walz, Ing.-Arch. 16 (1948). S, 243. 
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2. Niherungsverfahren. Es wird wieder die mathematisch einschneidende Annahme ge- 
macht, dafi§ auch bei kontinuierlicher Absaugung die Geschwindigkeitsprofile durch eine ein- 
parametrige Kurvenschar angenahert werden kénnen. Dann ist die Grenzschicht an jeder 
Stelle durch einen Formparameter und ein MaB fiir die Grenzschichtdicke bereits vollstandig 
heschrieben. Als Parameter werde benutzt 


mit 
0, = Verdrangungsdicke Il (2 -- Z| dy (2) 
0 
und 
0. = Impulsverlustdicke = fo (1 = a) dy , (3) 


wo y den Wandabstand, u die Geschwindigkeit innerhalb und U diejenige auSerhalb der 
Grenzschicht bezeichnen. Ein dimensionsloses MaB fiir die Dicke der Schicht ist die drtliche 
Reynolds-Zahl der Grenzschicht 


ctor (a 


mit »y = kinematische Zahigkeit. 

Zur Berechnung dieser zwei GréBen, H und re, abhangig von der Bogenlange s vom Stau- 
punkt, geniigen auch hier Impuls- und Energiegleichung, wenn auBer der Geschwindigkeits- 
verteilung aufSerhalb der Schicht U(s) auch die értliche Kérperradius r,(s) bei Drehkérpern 
und die Absaugeschwindigkeit v,(s) gegeben sind. 

Der Impulssatz lautet bekanntlich 


: i @ 0, dU v 
EO) seek (To de) + (2 A) ds aii (5) 


9 a ae a 
0 U? re Tapas 


und analog dazu der Energiesatz 
BaD ee MIE d 
Can Sree s 


(U* 1,63) — = : (6) 
Dabei bedeuten 1, die Wandschubspannung, 0 die Dichte, ¢ die dimensionslose Wandtangente 


des Profils 
A(u/U) 


E—E= a(y/0s) ° (7) 
und D eine dimensionslose Dissipationsfunktion. 
ioe ip r (du/U\? 7 y 8 
sa | ae as @) 
0 


mit r= Abstand von der Achse des Drehkérpers, und 6, = Energieverlustdicke: 


o-feh-ele : 


Bei einer einparametrigen Profilklasse bestehen nun eindeutige Beziehungen zwischen den 


verschiedenen, dimensionslosen Parametern wie H, e, D oder 


0 
Hy = 5 - (10) 


Man kann sie daher alle als eindeutige Funktionen eines Parameters, z. B. von H, auffassen: 


¢(H), D(H) und Hy, (H). 


Dann lassen sich aus den zunachst noch allgemein giiltigen Gleichungen (5) und (6) sofort 


-gwei Differentialgleichungen fiir die Unbekannten H und re ableiten, die in dimensionslosen 
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GréBen geschrieben lauten 


dH 1 d(U/U,) U Re eA CrT Re Uo ll) 
mince ailhin cE SOE ise: tne ( 

dre H+1aqU/O) , 1 dro/R) , 1, U Res py % 12 
do a feces enue | GES) Drs ina rel gee U (12) 


Darin bedeuten R eine charakteristische Linge, z. B. der gréBte Radius des Drehkérpers, 
cg =s/R die dimensionslose Bogenlange vom Staupunkt aus gemessen, U, die Anstrém- 


geschwindigkeit und 
Re — UR . (13) 


f, g und h sind Funktionen des Profilparameters H, die fiir die zugrunde gelegte Profilklasse 
28 


einmal berechnet werden miissen: 


(2) aw 
gH) = — 22 58M | | (14) 
Gz = == oe (15) 
und 
ois 


N= (16) 


dH * 

Zur schrittweisen Integration des Systems (11) und (12) denken wir uns, ahnlich wie in! 
alle auftretenden Funktionen der Bogenlange iiber das jeweilige Intervall linearisiert. Die 
Schrittlange As von der Bogenlange s, bis s,; muB so klein gewahlt werden, da®B alle Funk- 
tionen in diesem Intervall geniigend genau durch Geradenstiicke angenahert werden kénnen. 
Dann lassen sich Differentialquotienten durch Differenzenquotienten ersetzen, wie z.B. in 


ya ty me gl a8 a 
: =a) 
und 
ana Didtn 45 oh sale (18) 
nS — [(To)s = (Toe | 


Z 


Fiir die gesuchten Zuwichse von Are und AH ergibt sich dann 


—{(H +1)y + @m}re, + « wi Re Ao + “° Re Ao 
‘ U, rey U, 
Are= ee ———- =F — (19) 
Tan Say epee eae 
4 2 Ug re, 
und 
/ R 
Sra on ea Frey eine Rs 
: re? 
ae I lst) ee ce 
e e ° Cab 
» g— — Ao —— h— — A 
: zZ fr me inenes 3 ai 1h $ 


Der Querstrich ttber U, v, und re bedeutet Mittelbildung im Intervall also z.B. U = fed (U,+U,). 


Der Punkt bezeichnet Differentiation nach H, zB. f = df(H)/dH. Die Funktionen f, g und h 
sowie ihre Ableitungen sind in Tabelle 1 angegeben. Dabei wurden die Absaugeprofile von 
H. Schlichting® zugrunde gelegt, die zwischen dem Blasius-Profil (v, =0) und dem asympto- 
tischen Absaugeprofil interpolieren. 


1 Siehe FuBnote 9 S. 368. 
2 Siehe FuBnote 3 S. 368. 
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Tabelle 1. Hilfsfunktionen berechnet aus den Absaugeprofilen von H. Schlichting. 


ee ee re NF 
# ; ; a ee ete cs 
aaa | | eee CAP Bae é 

2 60°) 12-227, 60 «ht 3.19 Reet Lets 0,7 0.500 
2.05 a7 —1,48 ove Se) = 0.9 0,467 
21. OS oe 8589 Solel lua t05 L.1 0,437 
2.15 10.3 127% 3.49 8.9 1,05 | 1,3 0,409 
22 Lis —1,17 3,63 LO;0e i eG ba wr 1S 0,383 
2,25 12,4 -1,06 | 3,79 Us a feat eh 1,65 0,359 
2s 13,6 —0,95 | 3,96 1250 | abel 1,75 0,337 
2,35 14,9 —0,84 | 4,14 13,4. Lal 1,85 | 0,316 
2,4 16.3 —0,73 | 4,33 14,7 1,15 2,0 0,297 
2,45 17,8 —0,61 4,54. 16,2 125 2,25 0,279 
R5i 19.5 —0,48 1,78 183%, 105 2,6 0,262 
2,55 2s —0,34 | 5,06 20,6 | 1,4 3,0 0,246 
2,6 23.7 —0,20 | 5,38 23,0 | 1,45 3,00) eoU,ba) 
2,65 26,1 —0.06 73 25,5. | 1,5 3165 O27 
2.7 28.7 -0.10 6.11 7A aaa een 3,95 | 0,204 
25 oy 0.27 6,52 S100 2 is 4,25 0,192 
2,8 34,9 0.45 | 6,96 | 34,6 | 2,0 4,6 0,181 
2.85 38.7 0,67 | 7,44 40.0 | 2,4 , 4.9 0.170 


Die Formeln (19) und (20) gelten natiirlich auch fiir 1, —0, also fiir unbeeinfluBte Grenz- 
schichten. Jedoch ist dann ein etwas abgedndertes Gleichungssystem! geeigneter, das in 
Nahe des Staupunkts besser konvergiert und fiir das die Profilfunktionen f und g aus den 
Hartree-Profilen® berechnet worden sind. 


3. Stabilitatskriterium. Gewéhnlich interessiert nun weniger die Weiterentwicklung einer 
Grenzschicht bei vorgegebener Absaugung an sich, sondern die praktische Frage, wie weit 
der Widerstand eines Kérpers durch Laminarhaltung der Grenzschicht durch Absaugen 

_uberhaupt verringert werden kann. Nach J. Pretsch® wird der Reibungswiderstand, auch bei 
Einbeziehung der Absaugeleistung dann ein Minimum, wenn die oértliche Absaugegeschwindig- 
keit tberall so bemessen wird, da®B das laminare Geschwindigkeitsprofil der Gefahr des Um- 
schlags in die turbulente Form noch gerade widersteht. Die geringste dazu erforderliche 
Absaugung ist dann von Fall zu Fall erst zu berechnen. Zu den bisherigen Unbekannten re 
und H tritt also noch eine dritte — namlich v,(s) — hinzu, zu deren Berechnung neben Im- 
puls- und Energiegleichung als dritte Gleichung eine Stabilitatsbedingung benétigt wird. 
Ein Kriterium zur Bestimmung des wirklichen laminar-turbulenten Umschlagpunktes ist 
zwar noch nicht bekannt. Statt dessen rechnet man gewohnlich mit dem berechenbaren 
Indifferenzpunkt, d.h. mit dem Punkt, bis zu dem die laminare Grenzschicht gegen jede 
kleine Stérung stabil ist, und der noch stromaufwarts vor dem eigentlichen Umschlagpunkt 
liegt. Das Stabilitatskriterium hierfiir la8t sich fiir einparametrige Geschwindigkeitsprofile 
in der Form (U6,/7)krit. = Ortliche, kritische Reynolds-Zahl der Verdrangungsdicke abhangig 
von irgendeinem Profilparameter angeben. Fiir das vorliegende Verfahren, bei dem die 
Grenzschicht durch re und H charakterisiert wird, ist es praktischer 


peered: ee (21) 


als Funktion des Parameters H aufzufassen. 

In dieser Art sind die veréffentlichten Ergebnisse fir verschiedene Profilklassen in Abb. 1 
aufgetragen. Die vollen Kreise gelten nach J. Pretsch*®*® far Hartree-Profile (z.B. ebenes 
Staupunktprofil H =2,22; rotationssymmetrisches Staupunktprofil H 2,33; Platten- 
strémung H =2,59). Die Punkte fiir diese Grenzschichten an undurchlassigen Wanden liegen 
offenbar mit guter Annadherung auf einer Geraden, der eine Gleichung 


Carers (22) 


1 K, Wieghardt, Aero. Quart. V (1954), S. 25. 

2 D. R. Hartree, Proc. Cambr. Phil. Soc. 33. Part II (1937), S. 223. 

3 Siehe FuBnote 4 von S. 368. 

4 J. Pretsch, Jahrb. d, deutschen Luftfahrtforschung I (1941), 5. 58. 

5 J. Pretsch, Jahrb. d. deutschen Luftfahrtforschung I (1942). 5.54. 
6 J. Pretsch, Luftfahrtforschung 18 (1941), S. 341. 
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entspricht, mit 
d= 20,5 IO" beeen (23) 


Auch die kritische Re-Zahl fiir das asymptotische Absaugeprofil* 


fe tpl Ea (24) 
liegt noch in Nahe dieses Gesetzes. 
Etwas hihere re-Werte fiir gleiche H-Werte markieren die leeren Kreise, die Ergebnisse 
von H. Schlichting und A. Ulrich? auch ohne Absaugung wiedergeben; sie gelten fiir 
P,-Profile, d.h. fiir Geschwindigkeitspro- 
ele file, die aus gewissen Polynomen sechsten 
7 i: Grades bestehen. Etwa fir 2,3 << H <2,7 
ergibt sich in Abb.1 wieder eine Gerade 
mit den Konstanten 


Gia 29,1, | nd Se (25) 


2 Fur die hier benutzten Absaugeprofile 
» Martree-Profle[tetd] nach H. Schlichting schlieBlich haben D.L. 
EN Burrows, A.L.Braslow und N. Tetervin® 
© roratonssymmeprisches ? Fag s 

Staupunkiproti [14] die Stabilitatsrechnung durchgefihrt. Das 
Bs CITE Ergebnis ist in Abb.1 durch die ein- 
gezeichnete Kurve dargestellt. Auch diese 
Kurve wird im praktisch interessierenden 
Bereich von 2,2 <H <2,6 durch die Ge- 
rade (22) mit den Konstanten (25) genau 
genug angenahert. Diese Konstanten wer- 
den daher der weiteren Rechnung zugrunde 
gelegt. 


45 


40 


35 =i 


30 


Absaugeprotile [6] 


Die Ergebnisse von Stabilitatsrechnun- 
gen an ebenen Grenzschichten lassen sich 
nach J. Pretsch* ohne weiteres auf den 
rotationssymmetrischen Fall whbertragen, 
solange die Grenzschichtdicke klein gegen- 
uber dem Kérperradius bleibt. Daher kann 
das obige Kriterium auch auf Drehkérper 
angewandt werden. 


6 Zs 


4, Optimales Absaugen. Die Grenz- 
schichtrechnung fiir den Optimalfall der 
Absaugung kann nun folgendermafen 
durchgefiihrt werden. Zunachst wird vom 
‘-d/6, Staupunkt an, wo die Grenzschicht auch 

ohne Absaugung noch laminar ist, nach 

Abb. 1. Stabilitaétskriterium fiir verschiedene Profilklassen. einem der iblichen Verfahren gerechnet. 

Vom Indifferenzpunkt an, wo bei undurch- 
lassiger Wand die Gefahr des Turbulentwerdens beginnt, wird Absaugung angenommen, und 
zwar gerade so stark, da iiberall stromabwarts das Kriterium (22) erfillt ist. 


Fihrt man dazu diese Gleichung (22) als starre Bedingung in die Impuls- und Energie- 


gleichung (5) und (6) ein, so kann man jetzt re und v,(s) eliminieren, um eine Gleichung allein 
fir H zu erhalten: 


60 


60 Za a4 60 68 


dH __—sfh(H) dU, f,(H) dr, URED ee 
a ie Ue Tk de) a (re =e 


ee (26) 


re 


1 J. Pretsch, Jahrb. d. deutschen Luftfahrtforschung I (1942) S. 7. 

2 H. Schlichting und A. Ulrich, Jahrb. d. deutschen Luftforschung I (1942). S. 8. 
3 Siehe FuBnote 6 von S. 368. 

4 Siehe FuBbnote 6 von S. 371. 
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Daraus folgt wieder durch Linearisieren fiir kleine Schrittlangen die Rechenformel 


Re U 
—hy+ho—h a 40 
AH = peso: ee et (27) 
Ly Se pe cae Re U 
aie ne ae ahe Or Ey bh) Be Us, lo 


Nachdem so H(c) und damit auch re(o) berechnet sind, folgt die drtliche Absaugegeschwindig- 
keit aus dem Impulssatz 


“0 wt 1 dre fi MI dU Ida re (GP 9s 
U, Re do | ( U da ' TF | Reo tee: (28) 
Die Hilfsfunktionen f. fo und f, hangen dabei von H wie folet ab: 
jSpLEG ION 2 
AG 9 oe (29) 
1— H, 
NE eee ad (30) 
c¥_ 2D 
f(H) = 35 (31) 
mit 
dH. 


Fir die benutzten Absaugeprofile und die Konstante b =9 aus (25) erhalt man Funktionen, 
die sich fiir 2<H <2,8 genau genug durch lineare Ausdriicke annahern lassen, namlich durch 


= = 050121, 47 0.2745 (Hi 2,2), (33) 
fo ©  0,1210 + 0,0033 (H —2,2) , (34) 
fs ~  0,0385 (H —2). (35) 


Damit sind alle Formeln zur schrittweisen Berechnung der Absauge-Grenzschicht mit Hilfe 
des Rechenschiebers zusammengestellt. 


5. Anwendungsbeispiele. a) Plattenstrémung. Die Grenzschicht an einer langs an- 
gestrémten Platte (also bei konstanter AuBengeschwindigkeit U) mit optimaler Absaugung 
ist 1943 von J. Pretscht berechnet worden. Seine Lésung kann jetzt zur Kontrolle benutzt 
werden. 

Bei dieser einfachsten ebenen Strémung kann man die Gieichung (26) direkt integrieren; 
denn sie lautet hier 


dH U, 
re" ds fs - (36) 
Mit 
fs =¢(H—2) (35) 
und 
eee ae ee (22) 
wird daher 
H 
Ee pr Us 
Re(s) — Re(s,) = —ce ee rs tec (Re (se "i IE (37) 
H, 


Der Index 0 bezieht sich auf den Indifferenzpunkt, wo die Absaugung beginnt. Nachdem 
so ermittelt worden ist, wie der Parameter H mit wachsender Lauflange Re(s) abnimmt, 
kann das zugehdrige v,/U, aus Gl. (28) berechnet werden, die sich hier vereinfacht zu 

Ce (38) 


tai nae re 


Das Ergebnis ist in derselben Weise wie bei J. Pretsch, namlich —,/U, itiber V Re(s), 
in Abb.2 aufgetragen. Rechnet man mit den Konstanten a= 26,3 und 6 =8 fiir das Stabili- 


1 Siehe FuBnote 4 von S 368. 
26 
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tatskriterium, so wird —v,/U, tiberall zu gro. Das ist aber nach Abb.1 auch verstandlich: 
die Gerade mit diesen a- und b-Werten liegt deutlich tiefer als die Indifferenzkurve fiir die 


Abbsaugeprofile. Um diese unndtig scharfe Stabilitatshedingung zu erfiillen, mu eben starker 
abgesaugt werden, als eigentlich 

bul" erforderlich. 

Auch das andere Konstanten- 


4. | paar (a=29,1 und b=9) gibt 
| 


— 
if 


\ — JPretsthl¥] aus dem gleichen Grund noch 
se zu groBe Absaugegeschwindig- 
keiten, doch sind die Unter- 
ieoens } schiede im gezeichneten Bereich 
Paya Se — jedenfalls fiir Re(s) >107 — 
| es bedeutungslos. Andererseits 
re wird hier die asymptotische Ab- 
week saugegeschwindigkeit fiir Re(s) 
—oo wesentlich zu klein, im 
aa Ayptoter Gegensatz zur Rechnung mit 
05-10" 70° 15°10* Esjy VO” a= 26,3 und b=8. Auch das 
: ist nach Abb.1 erklarlich, da 
die obere Naherungsgerade fir 
H-> 2 wesentlich iiber den Wert 
(re)isit, fiir das asymptotische 
wt Absaugeprofil —hinausschieBt, 
wahrend die andere Gerade ganz 
gut dorthin trifft. 
Der Beiwert fiir den Rei- 
bungswiderstand der ganzen 
Platte ergibt sich aus 


2-10" 


\ °  berechner Milt refi =e s 
: 291-91 
(f \\ \ x »” ” TC Ly if =@ 


> 


{ — 
70° 70” 5-10" a 310° 


Abb. 2. Optimales Absaugen an der Platte. 


0° 
Re 
1 Tv) 
—————— -dRe 
f~ Re | @ ye 
Deas 
: Ugb/y Be 
0 YI elle 2 
70° 70° 70” 70° 70° ae ekH) dhe. (39) 
Ellipsoid 1:8 Ovold 1:8 NACA-Profiimessungen [67] Re re 
& Cy Messung [17] = Cy 0 
& op bel optimaler Absaugung +c be/ optim Abs. © Cy be/ optim. Abs. einschhebiich . : 
sunlelingmomie, “ Veaigatay Wenn sich auch in den Absauge- 
& Cy be/ opt Abs. ohne Abs-Leitung geschwindigkeiten gewisse Unter- 
Abb. 3. Vergleich mit Plattenwiderstand. schiede gegenitber der Lésung 
von J. Pretsch ergeben haben, 
; so stimmen diese cr-Werte und damit die 
jhe aw rechnerischen Ersparnisse im Widerstand 


3 


durch Laminarhalten durch Absaugen gut 
iiberein (s. Abb.3 und 4). 

b) Rotationsellipsoid 1:8. Nachdem 
sich das Verfahren im vorigen Beispiel als 
brauchbar erwiesen hatte, wurde es auf ein 
Rotationsellipsoid angewandt. Das Langen- 


Ersparnis iim Kelbungswiderstand 
S 
= 
| 


© Platte, repr =e 4 2 ee 

| : a e, rit yah verhaltnis gréBter Durchmesser/ Gesamtlange 
71 & 
WACA-Profil-Messungen [6,7] ist 2 R/L =1:8. 

© einsihleblich Absaugeleistung Als Re-Zahl wurde gewahlt: 

|a ome » ‘ 5 
Ull/p 3 j 
oe oe sa? es bezogen auf den gréBten Radius 


a= c 5 
Abb. 4. Ersparnis im Reibungswiderstand der Platte bei : Uy Riv 2 iy ‘ 
optimaler Absaugung. bezogen auf die Gesamtlange L 


P U,Liy = 8+ 10% 

ie giinstigsten Absaugegeschwindigkeiten sind in Abb.5 aufeetragen. V Indiff 

punkt der unbeeinflu8ten Grenzschicht am Kopf an bis etwa 34, is Linge jean th raaen 
( 0 
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Bits) 


konstant, und zwar 
Buel . 10-4 
ase Da kOe (40) 
Fir das hinterste Ende jedoch schreibt die Rechnung értlich sehr starkes Absaugen vor, 
wenn man auch dort nicht nur Ablésung sondern auch turbulentes Umschlagen der Reibungs- 
schicht vermeiden wollte. 

In die Berechnung des Absaugebeiwertes cg geht bei Drehkérpern auch die Kérperform 
ein. Wird bis zu einer Lauflange « lings der Achse abgesaugt, so ist der Mengenbeiwert wie 
folgt zu definieren: 

x/R 


Q ZT Gee ene S 
(=a = —— Oise Al 
c= FU, i Be ok ae 


Abb. 5. Optimale Absaugung ain Ellipsoid 1:8 bei Re = U, R/y = 5-108; U,L/vy = 8-108. 


mit Q=in der Zeiteinheit abgesaugtes Flissigkeitsvolumen, und / —gesamte Oberflache 
des Ellipsoids. Fiir das Seitenverhaltnis L/2 R=8 wird 


= 2 70 R aay e* + arcsin ele) a= 0 Syl 1K? (42) 


mit 
e = Exzentrizitat = J1 — (2 R/L). (43) 


Der in Abb.5 gezeichnete Verlauf von cg langs der Achse zeigt, daB die Absaugemenge erst 
spater als die Absaugegeschwindigkeit stark anwachst. 

Praktisch ist ein derart starkes Absaugen am hintersten Kérperende natiirlich undurch- 
fiihrbar; bei Drehkérpern wird man ja dort die bis dahin abgesaugte Fliissigkeit wieder nach 
hinten ausblasen miissen. Man mu also sowieso das Turbulentwerden und Ablisen der 
Grenzschicht dort zulassen. Zur Berechnung des Reibungswiderstandes des ganzen Kérpers 


iB, L/R 
i ( A aay ONC ees 
= ; =| 2 Eo toi X= Fm | = R(T, d= (44) 
5 USF 5 


ist aber der Anstieg der Wandschubspannung auf diesem letzten kurzen Teilstiick bedeutungs- 

los. Wesentlich kénnte nur der Formwiderstand sein, d.h. der Druckwiderstand, der sich 

daraus ergibt, daB am Heck der potentialtheoretische Druckriickgewinn (mit einem Stau- 
Jom 
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punkt auf der Achse) sich nicht einstellt, sondern statt dessen ein Nachlauf. Bei langen, schlan- 
ken Drehkérpern ohne Grenzschichtbeeinflussung ist dieser Formwiderstand jedoch Klein 
gegeniiber dem Reibungswiderstand. Der gemessene Gesamtwiderstand solcher Kérper 
hezogen auf die Kérperoberflache erweist sich namlich als nur wenig gréBer als der bezogene 
turbulente Widerstand einer gleich langen ebenen Platte. Man wird daher nur einen kleinen 
Fehler machen, wenn man auch bei laminar gehaltener Grenzschicht annimmt, daB der Ge- 
samtwiderstand fast dem berechenbaren Reibungswiderstand gleichzusetzen ist. Allerdings 
wirde hier schon ein kleiner Formwiderstand mehr ins Gewicht fallen als bei natirlicher, 
turbulenter Strémung. 

In Abb.3 wird nun der auf die Oberflache bezogene Beiwert des Reibungswiderstands des 
Ellipsoids verglichen mit dem der ebenen Platte. Der von H. Amtsberg! gemessene Beiwert 
des Gesamtwiderstands ohne Absaugen liegt nur wenig iiber dem der turbulenten Platten- 


Abb. 6. Optimale Absaugung am Ovoid 1:8 bei Re = U, R/y = 5-105; U, L/y = 8- 10°. 


strémung, und der bei optimaler Absaugung ist ebenfalls nur wenig gréBer als fir die ent- 
sprechende Plattenstrémung. Ebenso verhalten sich die rechnerischen Beiwerte fiir unbeein- 
fluBte, aber trotzdem laminar gebliebene Grenzschichtstrémungen an Ellipsoid und Platte. 


Diese Ergebnisse fiir das Ellipsoid 1:8 berechtigen wohl zu der sowieso plausiblen Verall- 
gemeinerung, das bei schlanken Umdrehungskérpern, deren Widerstand fast ausschlieBlich 
aus Reibungswiderstand besteht, durch Laminarhalten der Grenzschicht durch Absaugung 
Widerstandsgewinne von derselben GréBenordnung wie fiir die ebene Platte (s. Abb.3 und 4) 
méglich sein miSten. 


c) Ovoid (Quell-Senken-Kérper) 1:8. Als letztes Beispiel wurde noch das Ovoid 1:8 
fir die gleiche Re-Zahl U, R/y =5-10°, bzw. U, L/v =8-10° berechnet. Die Potentialstrémung 
wurde durch Annahme einer punktférmigen Quelle und gleichstarker Senke bestimmt. 

Die Aufengeschwindigkeit U steigt hier nicht wie beim Ellipsoid monoton bis zur Kérper- 
mitte an, sondern hat ein Maximum in Nahe der Kérperenden. Demzufolge ergibt sich fiir 
die optimale Absaugegeschwindigkeit kein so ausgeglichener Verlauf wie beim Ellipsoid, und 
im letzten Viertel des Kérpers schreibt die Rechnung teilweise sogar Ausblasen statt Ab- 
saugen vor (Abb.6). Das riihrt davon her, daf das Stabilitatskriterium (22) als starre Bedin- 
gung eingefiihrt worden ist und nun automatisch iberall erfillt wird. Dort, wo (potential- 
theoretisch) ein stabilisierender Druckabfall — also Ansteigen von U — vorhanden ist, ware 
die laminare Grenzschicht ohne Ausblasen stabiler als erforderlich, um das turbulente Um- 
schlagen zu vermeiden; die Grenzschicht ware dann unnétig diinn und die Wandschubspan- 
nung somit tiberfliissig groB. Sowie die Geschwindigkeit U jedoch wieder abnimmt, wird das 
starke Absaugen in Hecknihe wie beim Ellipsoid gefordert. 


' H. Amtsberg, Jahrb. d. Schiffbautechn. Gesellschaft 38 (1937), S. 177. 
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Trotz des unterschiedlichen Verlaufs der Absauge- bzw. Ausblasgeschwindigkeit lings des 
Kérpers ist der Beiwert des gesamten Reibungswiderstands wieder fast derselbe wie beim 
Ellipsoid, wenn man diesen auf die hier gréBere Oberflache bezieht (vgl. Abb.3 und Tabelle 2) 


Tabelle 2, Rechenergebnisse fiir zwei Drehkérper bei U, L/v == 8 - 10°. 


eS 


Rankinesches 
Rotationsellipsoid Ovoid (Quell- 
Senken-Kérper) 


Lange/Durchmesser L/2 R 


8 
Oberflache F 79,5 Re 06,5 R 
Volumen V Soh) ate AD Oe Ee 
Meridianflache M (groBte Wasserlinienfliche) Donan OH TRE 
Querschnittsflache (%) (Hauptspant) 3,14 R? 3,14R? 
a= M/2RL 0,785 | 0,942 
B = (2 RP | 0,785 | 0,785 
6=Vi(2 RPL | 0,524 0,711 
¢ = VIB L = 6/6 0.667 0,906 


Turbulenter Gesamtwiderstand bezogen auf die Oberfliche 
(Messung Amtsberg) 


ee : 
wus E Bye o IOS == 
Reibungswiderstand bei optimaler Absaugung bezogen auf die 
Oberflache 
Wr abe U8 F 0,67 1078 0,68, «10-3 
Laminarer Reibungswiderstand bezogen auf die Oberflache 
EEO as fio 
tom U3 I 0,513 - 10 


Turbulenter Gesamtwiderstand bezogen auf den Querschnitt 


>1Q 
Wis USS 84,0 -10- = 
Reibungswiderstand bei optimaler Absaugung bezogen auf den 
Querschnitt 
W, avs [SUR &) 17d. 10-8 Ziel One 


6. Zusammenfassung. Es wird ein verhaltnismafig einfaches Verfahren zur Berechnung 
von ebenen und drehsymmetrischen laminaren Grenzschichten mit Absaugung oder Aus- 
blasung beschrieben. Besonders geeignet ist es zur direkten Bestimmung des optimalen Falles, 
wenn iiberall am Kérper so abgesaugt werden soll, daB die laminare Grenzschicht gerade noch 
stabil ist. 


(Eingegangen am 19. November 1953.) 
Anschrift des Verfassers: Dr. K. Wieghardt, Hamburg 1, Berliner Tor 21. 
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Zur Thermodynamik des Uberschall-Geradstofdiffusors. 
Von E, Sanger. 


1. Problemstellung. Der Diffusor tibernimmt beim Uberschallflug die Rolle der Luft- 
schraube im Unterschallflug als eigentlich treibendes Organ. 

Die freifahrende Luftschraube verlor ihre Anwendungsméglichkeit praktisch bei Flug- 
Machzahlen iiber etwa 0,7. Die ummantelte Luftschraube wird in ihren mannigfaltigen Formen, 
z.B. als Kompressor der Turbinentriebwerke und Turbinenstaustrahltriebwerke noch bis zu 
Flug-Machzahlen von etwa 2 anwendbar, wobei jedoch ein wesentlicher Teil des Schubes 
schon nicht mehr an der Schraube, sondern an der als Diffusor wirkenden Ummantelung 
ibertragen wird. 

Bei weiterer Steigerung der Flug-Machzahl verschwindet die Schraube véllig aus den Luft- 
strahltriebwerken und diese degenerieren zum Staustrahltriebwerk, wobei der Diffusor allein 
den Triebwerksschub erzeugt. ‘ 

Unter den verschiedenen Formen der Uberschall-Diffusoren als Lufteinlaufe fiir Uber- 
schalltriebwerke scheint bisher hauptsichlich der Uberschall-GeradstoBdiffusor praktische 
Bedeutung zu erlangen, mit dem in 
den Jahren 1941—44 die ersten 
Unterschall-Flugschleppversuche 
mit Uberschall-Staustrahltrieb- 
werken!? in Hinblick auf die spa- 
tere Entwicklung von Uberschall- 
Staustrahl-Verkehrsflugzeugen * 
| ausgefiihrt wurden. 

A EN Querschnitts Machen, me ip Die grundsatzlichen Uberlegun- 


ie Tig fay Dritcke, kg fre, ; fe gen, die zu dieser damals noch sehr 
Qr Ox O,' Dichter kg seKi/m Ge ungewohnlich erscheinenden Diffu- 
i le Te Jemperaturen, °K h f, fahrt b 5 : 
cP Up 70,! Geschwindigke/ten, m/sek 7) SOFLOrM bes FLOR, Desc enen ae 
thy Ag, Oy! Schallgeschm. m/sek ae schen steigendem Interesse, auch 

10, /tty>1 Mp lg, Ty! Stromungsmachzahlen — 70,/z fir andere Strahltriebwerksarten, 
Hy Hy Hy! Adiabatenexponenten, — 22 


und sind deshalb in der vorliegenden 
Abb. 1, Uberschall-GeradstoBdiffusor. Arbeit zusammengestellt worden, 
zu den Grundlagen einer quanti- 

tativen, thermodynamischen Theorie des Uberschall-GeradstoSdiffusors. 

Abb.1 zeigt eine Ausfiihrungsform desselben mit einigen Kennzeichnungen der Innen- 
strémung und der verwendeten Formelzeichen. 

Die Uberschallstrémung wird in einem vor dem Diffusormund F\, liegenden senkrechten 
Verdichtungssto8 in F in Unterschallstrémung verwandelt und der Diffusor selbst daher 
mit Unterschallgeschwindigkeit angestrémt und durchstrémt. Hinter dem Strémungsquer- 
schnitt F', folgen meist Kompressoren, Aufheizung der Luft usw., die uns hier nicht inter- 
essieren. 

Ebenso soll die AuBenumstrémung des Diffusors hier unerértert bleiben, d.h. es wird nur 
die Wirkung der innerhalb des Einstrémzylinders F, flieBenden Luft betrachtet. Nur soweit 
die Schuberzeugung durch diese Innendurchstrémung in einigen Fallen auch die Aufen- 
umstrémung berithrt, wird diese mit in den Kreis der Betrachtungen gézogen. 

Die Besonderheit des Uberschall-GeradstoBdiffusors, da® er nach dem vorne liegenden 
Verdichtungssto8 schlieBlich mit Unterschallgeschwindigkeit wie die bekannten Rudenschen 
Fangdiffusoren angestrémt wird, hat zur Folge, daB sich, ahnlich wie bei diesen, der Luft- 
durchsatz bis zu gewissem Grad, trotz fester Bauquerschnitte, an die Betriebsbedingungen 
anpassen kann. Die Grenzen dieser automatischen Regelfahigkeit werden untersucht. 


1 EK. Sanger u. I. Bredt. Deutsche Luftfahrtforschung UM 3509, Ainrine 1943. — NACA, T. M. 1106, 
Washington 1947. a 

* I. Sdnger- Bredt, VDI-Forschungsheft 437, Diisseldorf 1953, 
° FE. Sanger, Z. Weltluftfahrt, 1 (1949), S, 167, 


| 
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Wahrend die Luftschraubentheorie im wesentlichen mit den Gesetzen inkompressibler 
Strémung, also der Aerodynamik, auskommt und nur in beschrinktem Umfang zu den Ge- 
setzen der kompressiblen Strémungen, also der Gasdynamik, greift, wird die Theorie der 
Uberschalldiffusoren auch die Veranderlichkeit der spezifischen Warmen und unter Um- 
standen sogar der chemischen Zusammensetzung der Luft beachten miissen, also die Gesetze 
der Thermodynamik bzw. der Thermochemie verwenden. 


Die Strémungsuntersuchung erstreckt sich auf drei deutlich getrennte, einzelne Strémungs- 
phasen: 


1. den geraden Verdichtungssto8 im Querschnitt F’, mit dem Ubergang der mit dem In- 
dex (1) bezeichneten ZustandsgréBen in jene mit dem Index (x), 

2. die isentropische Unterschallanstrémung zwischen den Querschnitten F, und Fy, wobei 
die Zustandsgréfen von (x*) nach (1’) iibergehen, und schlieSlich 

3. die adiabatische Einstrémung zwischen den Querschnitten Fy und F,, mit dem Uber- 
gang von den ZustandsgréBen des Strémungsmediums von (1’) nach (2). 

Die Gleichungen wurden derart bezeichnet, daB jene in thermochemischer Schreibweise 
nur Nummern fiihren, jene in thermodynamischer Schreibweise zur Nummer den Zusatz- 


buchstaben (a), jene in gasdynamischer Schreibweise zur Nummer den Zusatzbuchstaben (b) 
fiihren. 


2. Luft als StrOmungsmedium. Alle Strahltriebwerke arbeiten mit gasférmigen Arbeits- 
medien, die aus veranderlichen Gemischen von Photonen, Elektronen, anderen Elementar- 
teilchen, lonen, Atomen, Radikalen, Molekiilen, sogar zusitzlich Staubteilchen oder Fliissig- 
keitstrépfchen bestehen kénnen, und die nicht nur kompressibel sind, sondern wahrend der 
Strémungsvorgange chemische und physikalische Veranderungen und Verzégerungen in 
der Gleichgewichtseinstellung ihrer zahlreichen energetischen Freiheitsgrade aufweisen, end- 
liche freie Weglangen zwischen den Zusammenstifen ihrer Partikel und mit der Licht- 
geschwindigkeit vergleichbare Bewegungsgeschwindigkeiten derselben haben kénnen. 

Wieweit solche Gase im thermodynamischen Sinn als ideal zu betrachten sind, also der 
klassischen Zustandsgleichung pV=RT gehorchen, wird daher von Fall zu Fall sehr zu 
prifen sein. 

Wesentliche Voraussetzungen idealer Gase sind bekanntlich: 

1. Die Abmessungen der Gaspartikel sind unendlich klein gegeniiber ihren gegenseitigen 
Abstanden. Diese Voraussetzung ist um so besser erfillt, je geringer die Gasdichte ist. 

2. Die Gaspartikel titben aufeinander auf erhalb ihrer Zusammensto®e keine Krafte aus. 
Diese Voraussetzung ist ebenfalls um so besser erfillt, je geringer die Gasdichte ist, aber auch, 
je hoher die Gastemperatur ist, da mit wachsender Geschwindigkeit der Partikel die ablen- 
kende Wirkung der Wechselwirkungskrafte immer geringer wird. 

Wenn ideale Gase entarten, d.h. statt der Boltzmann-Statistik einer andern Statistik 
gehorchen, tritt an die Stelle der klassischen Zustandsgleichung eine modifizierte Zustands- 
gleichung. Dies kommt praktisch jedoch bei héheren Temperaturen nur in Photonen- und 
Elektronengasen vor. 

Wenn bei idealen Gasen die Partikelgeschwindigkeiten sich der Lichtgeschwindigkeit 
nahern, muf die klassische Zustandsgleichung durch eine relativistische Zustandsgleichung 
ersetzt werden. Dies tritt praktisch jedoch nur bei Photonengasen und bei Elektronengasen 
sehr hoher Temperatur auf. 

Wesentliche Eigenschaften der nichtentarteten, nichtrelativistischen, idealen Gase sind 
bekanntlich: 

1. Sie gehorchen der Zustandsgleichung pV = RT. 

2. Enthalpie und innere Energie hangen nur von der Temperatur, nicht vom Druck ab, 
d.h. auch beide spezifischen Warmen sind nur temperatur- und nicht druck-abhangig, ihre 
Differenz ist cp —¢, = R. ; 

3. Ihre Partikelgeschwindigkeiten folgen der Maxwell-Boltzmannschen Geschwindigkeits- 
verteilung. . 

4. Sie erfillen das Avogadrosche Gesetz, wonach im gleichen Volumen immer gleiche 
Partikelzahlen von jeder Gasart enthalten sind, die Gaskonstante pro Mol ist fiir alle Gase 

leich. 
3 5. Als Folge der Gaszustandsgleichung gilt die Poissonsche Adiabatenbeziehung pv* 
—konst., die allerdings auch temperaturunabhangige spezifische Warmen voraussetzt. Bei 


- 
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temperaturveranderlichen spezifischen Warmen kann man jedoch naherungsweise mit Mittel- 
werten fiir die in Frage kommenden Temperaturbereiche rechnen. 

Wieweit diese Voraussetzungen und Eigenschaften fiir technisch verwendete Gase zu- 
treffen, zeigt folgende Aufzahlung: 

1. Einfache Molekil-, Atom-, Elektronen-, Neutronen-, usw. -Gase mit Anregung der 
Translations- und gegebenenfalls Rotations-Freiheitsgrade haben konstante spezifische 
Warmen und geniigen bei geringer Dichte der klassischen Gaszustandsgleichung als Grenz- 
gesetz (praktisch ideale Gase, gasdynamische Strémung). 

2. Einfache Molekiil- oder Atom-Gase mit zusatzlicher Schwingungs- oder Elektronen- 
anregung haben temperaturabhangige spezifische Warmen und geniigen der Gasgleichung 
wie unter 1. (sogenannte halbideale Gase, thermodynamische Strémung). 

3. Gasgemische, deren Komponenten die Eigenschaften der unter 1. und 2. genannten 
einfachen Gase haben, und die nicht miteinander reagieren, sind ganz- bzw. halb-ideale Gase 
(gasdynamische bzw. thermodynamische Strémung). 

4. Gasgemische wie 3., deren Komponenten jedoch untereinander chemisch so reagieren, 
daB die Molekiilzahl und damit Molekulargewicht und Gaskonstante erhalten bleiben (z.B. 
thermisches Wassergasgleichgewicht oder Stickstoff-Verbrennung), haben bekanntlich vom 
Druck unabhingiges thermisches Dissoziationsgleichgewicht. Bei geniigend geringer Dichte 
gilt fiir Zustandsinderungen so dissoziierender Gasgemische die Gasgleichung, sie verhalten 
sich wie halbideale Gase, da die spezifischen Warmen mit der Veranderung der chemischen 
Zusammensetzung und der sie begleitenden Warmeténungen temperaturabhangig sind, selbst 
wenn die spezifischen Warmen der Komponenten es nicht sind (thermodynamische Strémung). 

5. Einfache Gase und Gasgemische bei hohen Dichten (z.B. im Van der Waalschen Be- 
reich) geniigen der Gasgleichung nicht (nichtideale Gase). Bei extrem hohen Temperaturen 
kann sich ein Gas jedoch noch bei der hundertfachen Dichte des Wassers wie ein ideales Gas 
verhalten, besonders weil die Atomabmessungen durch weitgehende Ionisierung stark ver- 
mindert sind. 

6. Gemische von im chemischen Gleichgewicht befindlichen, reagierenden Gasen haben 
im allgemeinen infolge der Verdnderlichkeit der Molekilzahl und der Warmeténungen der 
Reaktionen druckabhangige spezifische Warmen und genigen insgesamt der Gasgleichung 
nicht, wahrend die einzelnen Komponenten des Gemisches die Gasgleichung erfiillen. Die 
Gaskonstante des Gemisches ist wahrend der Zustandsanderung mit dem Molekulargewicht 
veradnderlich (praktisch nichtideale Gase, thermochemische Strémung). 

7. Thermisch ionisierende Gase sind ein Gemisch von Molekil-, Ionen-, Atom-, Elektronen- 
und meist auch Photonen-Gas (Plasma). Mit dem Molekulargewicht ist die Gaskonstante 
verdnderlich, nach der Eggert-Saha-Gleichung sind die spezifischen Warmen des Gemisches 
druckabhangig. Ebenso wie dissoziierende Gase sind also auch ionisierende Gase als praktisch 
nichtideale Gase zu behandeln (Thermoelektrische Strémung, Plasmastrémung). 

Die nun hier interessierende atmospharische Luft bei Flug-Machzahlen bis etwa 10 fallt 
nur fiir Flug-Machzahlen bis etwa 2 oder fiir auBerordentlich kurzzeitige Zustandsanderungen, 
bei denen sich die Tragheit aller inneren Freiheitsgrade gegeniiber dem Translationsfreiheits- 
grad stark geltend macht, unter den Beispielfall 1. ganz idealer Gase und folgt dann den 
Gesetzen gasdynamischer Strémung. ? 

Bei langsameren Zustandsinderungen, die die Einstellung thermischer Gleichgewichte 
erlauben, ist die atmospharische Luft bis etwa zu Flug-Machzahlen von sechs nach Beispiel- 
fall 2. als nur halbideales Gas zu betrachten, folgt also dann den Gesetzen thermodynamischer 
Strémung. 

Bei héheren Flug-Machzahlen und besonders bei geringen Luftdichten gelangt sie ohne 
weiteres in den Bereich der Falle 6. und 7. nichtidealer Gase, also thermochemischer bzw. 
thermoelektrischer Strémung. 

Da zumindest die Einfiihrung halbidealer Gase, also die Ablésung der gasdynamischen 
durch die thermodynamischen Betrachtungsweisen nicht nur fir Feuergasstrémungen, son- 
dern bei Flug- Machzahlen iiber 2 auch schon fiir die reine Luftstrémung nétig ist, zeigt Abb.2. 

Dort sind die wahren und mittleren spezifischen Warmen und die entsprechenden Isen- 
tropenexponenten der Luft mit Beriicksichtigung sowohl der nur temperaturabhangigen Frei- 
heitsgrade, wie Molekiilschwingung und Elektronentermanregung, als auch der temperatur- 
und druckabhangigen Freiheitsgrade, wie Dissoziation und Ionisation, zahlenmabig tiber der 
Temperatur aufgetragen, unter der iitblichen Annahme unendlich schneller Einstellung aller 
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inneren Gleichgewichte der Luftmolekiile, und fiir die druckabhangigen Bereiche fiir etwa 
Lat Druck und fiir den Druck unendlich. 

Gleichzeitig sind die beim Stau der Luftstrémung in der Stratosphare auftretenden Stau- 
temperaturen zu den Flug-Machzahlen eingetragen. 

Der normalerweise mit x —1,4 angenommene Adiabatenexponent der Luft kann also 
nicht nur fiir sehr rasche Zustandsainderungen, wie sie etwa in VerdichtungsstéBen auftreten, 
bis zu x =1,67 anwachsen, er kann fiir langsame Zustandsanderungen bei hohen Flug-Mach- 

_zahlen ohne weiteres auch bis gegen x =1 absinken. 

Fir den Arbeitsproze8 von sol- 
chen Uberschall-Strahlantrieben, — {%744/% 
die die Umgebungsluft mitbeniitzen, %) ‘ 
sind diese Umstinde von grund- 
legender Bedeutung. 

Die herkémmliche Behandlung 
der thermodynamischen Strémung 
in Strahltriebwerken mit Hilfe der 
Poissonschen Adiabatengleichung 

pon* = konst. (1b) 
ist also zweifellos eine sehr rauhe 
Naherung, da diese Beziehung ei- 
gentlich nur fiir ganzideale Gase, 
also fiir druck- und temperatur- 
unabhangige spezifische Warmen 
gilt und einen Grenzfall der Gas- 
gleichung fiir Zustandsinderungen 
bei konstanter Entropie, dS =0, 
darstellt. 

Demgegeniiber haben wir fest- 

-gestellt, daB sich die Feuergase 
der Strahlantriebe in der Regel wie 
nichtideale Gase, die Frischluft 
bestenfalls wie ein halbideales Gas 
verhalten. 

Beim Vorliegen entsprechender 
Zustandsdiagramme der Arbeits- 
medien, also auch der Luft, in 


7 0 
solcher Form, daB die druck- und yl 1! ale | 
a . 000 
temperatur-abhangigen x-, c)-, ¢,- 
. - z Abb. 2. Wahre und mittlere spezifische Warmen und Adiabatenexponenten 
und. Molekulargewichts Werte so der Luft und Flug-Machzahlen in der Stratosphire, tiber der Stautemperatur 
wohl die Warmeaufnahme durch (bis 2000° K druckunabhangig, dariiber fir die Driicke 1 at und unendlich), 


Dissoziation und Ionisation, als auch 
jene durch molekulare Schwingungsanregung und Elektronenanregung enthalten, kann man 
fiir das jeweilige Stromungsmedium kleine Bereiche geniigend konstanter M und geniigend 
konstanter x oder c,, c, aussondern, dort zulassige Mittelwerte innerhalb des zu erwartenden 
Druck- und Temperatur-Bereiches wahlen und dann die klassische Adiabatenbeziehung sowie 
die mit ihrer Anwendung hier zu entwickelnden Strémungsbeziehungen beniitzen, die unter 
solchen Voraussetzungen also auch fiir halb- und nichtideale Gase mit zulassiger Naherung 
giiltig sind?. 

Insbesondere lassen sich fiir den Bereich thermochemischer Strémung an Stelle der Adia- 
batenexponenten x =c,/c, Isentropenexponenten o und rt nach den Bestimmungsgleichungen 

ees konst. 


t—l1 


iy 
fiir bestimmte Strémungsmedien tabellieren, die die formale Weiterverwendung der Poisson- 
schen Gleichung gestatten. 


pV°= konst. baw. 


1 J, Sénger-Bredt, Die Rolle des Adiabatenexponenten bei chemisch verinderlichen Gasgemischen 
(im Erscheinen). Zeitschrift fiir angew. Mathematik u. Physik, Basel. 
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Die so erzielte Erweiterung des Anwendungsbereiches der klassischen Strémungsbeziehun- 
gen ist deshalb duferst erwiinscht, weil eine genauere analytische Darstellung der Zusammen- 
hange meist zu nicht expliziten und nicht integrierbaren Beziehungen fiihrt. 

Schon die Einfihrung der Temperaturabhangigkeit der spezifischen Warmen halbidealer 
Gase nur mittels der Planck-Einsteinschen Quantenbeziehung fiir die Molekilschwingung 
bringt sofort auSerordentlich unhandliche Ausdriicke* ae we Die Einfihrung der Elektronen- 
anregungs- oder gar Dissoziations- und Ionisations-Bezichungen in die Strémungsgleichungen 
erscheint noch viel umstindlicher, so daB die Naherungsrechnung hier noch wesentlichere 
Vereinfachungen verspricht, wenn man nicht zu graphischen Verfahren tibergeht. ; 

Die in geeigneten Mittelwerten beniitzten Adiabaten- oder Isentropen-Exponenten eines 
Gasgemisches, in dem alle genannten inneren Freiheitsgrade aktiviert sind, vernachlassigen 
immerhin noch die Schnelligkeit bzw. Vollstandigkeit der Gleichgewichtseinstellung zwischen 
diesen Freiheitsgraden, iiber die noch sehr wenig quantitative Kenntnisse vorliegen, und 
etwaige Kraftwirkungen der Partikel untereinander, z.B. in Plasmen. 

Trotzdem erlaubt diese naherungsweise Untersuchung der Strahlantriebe so interessante 
Einblicke in viele Einzelheiten ihres Verhaltens, und eine so groBe Ubersicht tiber ihre Eigen- 
schaften, da®B sie auch in dieser Arbeit ausschlieBlich verwendet werden soll. 

Die wahren und mittleren spezifischen Warmen und Adiabatenexponenten unter AusschluB 
der druckabhangigen Energieaufnahmen findet man fiir viele interessierende Gase und Gas- 
gemische, speziell auch fiir Luft und fiir die Feuergase der Kohlenwasserstoffverbrennung in 
Tabellen fertig vor. 

Entsprechende, wesentlich umfangreichere Tabellen mit Einschlu8 von Dissoziation und. 
[onisation diirften fiir viele Arbeitsgase und elektrische Plasmen in absehbarer Zeit dank 
der modernen. GroBrechenmaschinen verfiigbar werden. 

Fir die Untersuchung einer Zustandsanderung zwischen den Punkten 1 (Temperatur Tj, 
Druck p,) und 2 (Temperatur T,, Druck p,) des Zustandsdiagrammes Abb.3 beniitzen wir 

BS im allgemeinen dreierlei Werte der spezifischen 
hy iS & Warme c, bzw. der Adiabatenexponenten x 
\ oder Isentropenexponenten o und T. 

1. Die wahre spezifische Warme in den 
Punkten 1 und 2, namlich ¢,, =(dJ/dT),, und 
Chg =(dJ/dT),,. (z.B. zur Festlegung der je- 
weiligen Schallgeschwindigkeit a? =o, p,/0, baw. 
at —2,gR,T,. 

2. Die mittlere spezifische Warme 


dig 
J) Repel ee 
0 


zwischen dem absoluten Nullpunkt der Tempe- 
raturskala und den jeweils interessierenden 
Punkten 1 oder 2, namlich cp, und ¢Cyo (z. B. 
Abb. 3, Im Rechengang beniitzte spezifische Warmen der Luft. Aur Berechnung der Enthalpien J, — Cp OL T;, 
Jy = Cpo9T)). 


Ts 
3. Die mittlere spezifische Warme ¢p,, == [opd Ty (Ty 05) a lengeheys) 2 zwischen den 
T 


1 
Punkten | und 2, naherungsweise als arithmetisches Mittel der wahren spezifischen Warmen 


in den Endpunkten der Zustandsinderung (z. B. fiir Zustandsanderungen zwischen den Punk- 
ten 1 und 2). 


3. Der gerade Verdichtungssto8. Wenn man einem Gas erhebliche Warmebetriage zufiihrt, 
so stehen dessen Molekiilen hauptsachlich die thermischen, chemischen, optischen und elek- 
trischen Freiheitsgrade der Translation, Rotation, Oszillation, Dissoziation, Termanregung 
und Tonisation zur Aufnahme der Energien zur Verfiigung. 

Die GréBen der dem Molekiil zufiihrbaren Energiequanten wachsen im allgemeinen in der 
angegebenen Reihenfolge der Freiheitsgrade. Uber die Verteilung der gesamten zugefiihrten 


1 E, Sanger u. I. Bredt, Deutsche Luftfahrtforschung UM 3509, Ainring 1943, usw. 
_* EB, Sanger u. 1. Bredt, Deutsche Luftfahrtforschung UM 3538, Ainring 1944; Bureau of Aeronautics, 
Navy Dep. Transl. CGB,32, Washington 1946; — Ed. R. Cornogg, Santa Barbara, Calif. 1952. 
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Energie auf die einzelnen Freiheitsgrade nach dem Einspiclen. voélligen Gleichgewichtes sind 
wir fiir technische Zwecke ausreichend im Bilde, insbesondere durch den Boltzmannschen 
Gleichverteilungssatz, durch die Planck-Einstein-Funktion der Molekilschwingung, das Guld- 
bergsche Massenwirkungsgesetz fiir chemische Gleichgewichte, die Maxwell-Boltzmannsche 
Geschwindigkeitsverteilungsfunktion fiir Elektronenanregung, die Eggert-Saha-Gleichung fiir 
Tonisationsgleichgewichte usw. 

Dieser Gleichgewichtszustand, den wir unseren normalen, thermodynamischen Berech- 
nungen meist zugrundelegen, wird aber unter Umstianden erst nach vielen Millionen gegen- 
seitiger Molekilsté%e erreicht und kann von dem Gaszustand unmittelbar nach der Energie- 
zufuhr ganz erheblich abweichen. 

Die Richtung dieser anfanglichen Abweichung vom Gleichgewichtszustand wird auch von 
der Art der Energiezufuhr abhangen und oft verschieden sein, bei Aufheizung durch Kom- 
pression des Gases, durch exotherme chemische Reaktionen im Gas, durch Absorption elektro- 
magnetischer Strahlung oder durch Abbremsung schneller Partikel, wie Elektronen oder 
Nukleonen usw. 

Im Falle einer raschen adiabatischen Kompression im Gas, wie sie beispielsweise durch 
hochfrequente Schallwellen oder durch VerdichtungsstéBe erfolgen kann, wird man oft an- 
nehmen diirfen, daB die Zahl der notwendigen MolekiilstéBe zur Einstellung des thermischen, 
chemischen, optischen und elektrischen Gleichgewichtes in der Reihenfolge der obengenannten 
Freiheitsgrade stark anwachst, wobei allerdings die mit geringen Massentragheiten gekoppelten 
elektronischen Freiheitsgrade anscheinend bevorzugt werden. 

Man hat Griinde, zu glauben, daB der Vorgang des VerdichtungsstoSes selbst sich auf ganz 
wenige MolekilstéSe konzentriert, also ebenso schnell ablauft, wie die Einstellung der Maxwell- 
Boltzmann-Verteilung der Translationsgeschwindigkeiten der Partikel, was bedeuten wide, 
daB bei dem eigentlichen Sto, zumindest wenn dieser sehr intensiv ist, keine wesentliche 
Anderung des Energieinhaltes der inneren Freiheitsgrade von Rotation aufwarts eintritt, und 
deren Gleichgewichte sich erst wesentlich nach dem Stof stetig neu einstellen, wahrend 
zunachst nur der nicht gequantelte, translatorische Freiheitsgrad am Sto beteiligt ist. 

Die Tatsache, daf die Relaxation der Gleichgewichtseinstellung schon bei kurzen Schall- 
wellen gerade erstmalig beobachtet wurde, laBt solche Annahmen besonders wahrscheinlich 
erscheinen. 

In diesem Fall verhalt sich das Gas beim eigentlichen StoB wie ein einatomiges mit x = 1,67, 
beim Durchgang durch den Sto8 steigt seine Translationstemperatur gemaB diesem Adiabaten- 
exponenten, wahrend die Temperaturen aller andern Freiheitsgrade zunachst véllig unver- 
andert bleiben. 

Die Enthalpiesteigerung der Luft wahrend des Stofes kommt zunachst nur den Trans- 
lationsfreiheitsgraden zugute, die Mehrenthalpie wird sich erst nachher allmablich auf die 
ubrigen Freiheitsgrade verteilen, bei sehr hohen Temperaturen vorweg auf Elektronenanregung 
und Ionisation, weiterhin der Reihe nach auf Rotatiou — wobei der Adiabatenexponent 
der Luft von den anfanglichen x = 5/3 = 1,67 absinkt auf das tibliche x = 7/5 = 1,40 : 
weiterhin gegebenenfalls auch auf Molekiilschwingung — wobei x weiter absinken kann auf 
ye =o) (== 1,29 (siehe Abb. 2). Falls die Sto%-Machzahlen ausreichend hoch sind, kann mit 
der weiteren Einstellung von Dissoziation und Ionisation ins Gleichgewicht der Isentropen- 
exponent noch weiter bis gegen x = | sinken. 

In den bekannten gasdynamischen Beziehungen des senkrechten VerdichtungsstoBes 
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T. ae ( + 1)? (w?/a? — 1) 1+ — at/wi ae ere sth (5b) 
1 | ———__ w} /a, 
z+] 
Ts (6b) 
aa is 
a?/w? + x—1 
wily — (wx /wy)* -. hes 2 es (7b) 
w} /az Tx/ T, B w/a a a 
ace Qx wy /2 — es (8b) 
qi 0 3/2 Wy 


erscheint als Temperatur nur die Translationstemperatur, daher wird die Temperaturstei- 
gerung durch den Sto8 mit x = 1,67 wesentlich gréBer sein, als aus der iiblichen Rechnung 
mit x = 1,40 folgt, und diese ttberhéhte Temperatur wird erst nach dem Sto8 allmahlich auf 

die Gleichgewichtstemperatur 


Mp /M, ae absinken. Das nach dem Sto8 
Pfs | re ess coed jedenfalls mit Unterschall- 
CAE | = | ae | j geschwindigkeit strémende Gas 
I wird infolge der nachtraglichen 
0} x ‘Vas Ty, a Abkithlung auch noch langsamer 
| werden. Diese unmittelbar auf 
M0} 10 ieee ~~ den eigentlichen StoB folgende, 
ff verzogerte Zylinderstrémung hat 
oor 09 ane natirlich eine weitere Anderung 
5 |f MPa der Impuls- und Druckverhalt- 
Lo/Fy St 98 — a } y, nisse nach dem eigentlichen StoB 
30 i= / / zur Folge, die in Abb. 4 dar- 
0 O7 Zeitliche Feinstruktur des : 7 - gestellt ist. 
Verdichtungsstobes fir rvy/4=10 yf y Fiir die Berechnung der Gas- 
sO af zustandsgréBen unmittelbar nach 
\ / / dem eigentlichen Sto lassen sich 
or \ | / a die Gleichungen (2b) bis (8b) mit 
\ y, vi % = 1,67 verwenden. 
nes \ vy, se Die quantitativen Ergebnisse 
: 0 y a fel dieser Rechnung zeigt Abb. 4 fir 
ie Y ye den Flug- Machzahlbereich bis 10 
CA \ PSA aca Se wail EPR RS RD) in gestrichelten Linien. 
70} \ he Ne Ne Ble Man sieht, da8 z.B. bei der 
ss Noy | ie Flug-Machzahl 10 und der Stra- 
Soe tospharentemperatur 
10'- 07 yo Sane Sa) een = T,= 216,5° K 
| Ys ——|nja, die Spitzentemperatur voriber- 
ila ai Faas SaaeM | iatair Rae) MOREY EET ASO 5 A HG) gehend bis auf 
(Gesu) wad Wey nach Eunecinae arivnent pone tg nie at On elles 
ansteigt. 


Diese starke Energiekonzentration im Translationsfreiheitsgrad des Gases bei vor dem 
Sto® hohen Uberschallgeschwindigkeiten zeigt erfahrungsgemaB die Neigung, sich erst auf 
die mit nur geringen Massen verbundenen elektronischen Freiheitsgrade der Termanregung 
und Ionisation, und erst spater auf die mit den groBen Kernmassen verkniipften, trageren 
Freiheitsgrade der Rotation, Oszillation und Dissoziation zu verteilen, und zwar unter guter 
Kinhaltung der Maxwell-Boltzmannschen Verteilung der Translationsgeschwindigkeiten, und 
der zugehérigen Termanregungen und Ionisationsgrade. 

Daraus folgt eine der hohen Spitzentemperatur entsprechende, auSerordentlich starke 


fonisierung und sehr intensives Lumineszenzleuchten, beides wohlbekannte Begleiterschei- 
nungen sehr intensiver VerdichtungsstéRe in Luft. 
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Diese Spitzen der Translationstemperatur erlauben in bekannter Weise! die Intensitat des 
Lumineszenzleuchtens beim StoB und den Grad der zu erwartenden Ionisation abzuschatzen. 
Der so berechnete Ionisationsgrad bestitigt sich durch Messung des Reflexionsvermégens 
gegen kurze Funkwellen, das bei intensiveren VerdichtungsstéBen jenem massiver Metalle 
nahekommt, wodurch andererseits auch wieder die hier beniitzten Rechenvoraussetzungen 
indirekt bestatigt werden. 

. In der Gleichung (3’b) war p? der Ruhedruck nach dem Sto’. Da im Grenzfall einer 

gegen Unendlich anwachsenden Anfangs-Machzahl w,/a,—> oo die Strémungs-Machzahl wx/ax 

nach dem Stof gegen den endlichen Wert (x — 1)/2 x geht [siehe (7b)], findet nach dem Stof 
tad 


noch ein sehr merkbarer isentropischer Stau (p% [Pe)wja,> 0 = Le 1) jae _* statt. Dies 
folgt aus der Tatsache, daB bei w,/a,—> oo nach (4b) fiir x»+-1 die Dichte o, durchaus 
endlich bleibt, lediglich der Druck p, und die Temperatur T,, als auch die Geschwindigkeit w, 
werden unendlich groB, damit auch die Schaligeschwindigkeit a,, und dadurch schlieBlich 
die Machzahl w,/a, gleich Null nur fiir x = 1. 

Fir die Berechnung der GaszustandsgréBen T,, px, 0x usw. nach dem Sto8 und nach 
Einstellung aller Gleichgewichte, die wir als AusgangsgréBen fiir die weitere Strémungs- 
rechnung durch den Diffusor brauchen, mufs wegen der geschilderten Veranderlichkeit der 
x-Werte auf thermodynamische bzw. thermochemische Strémungsbeziehungen zuriick- 
gegriffen werden. 


Die bekannten Strémungsgrundgleichungen fiir 


Kontinuitat : 01 Wy = Ox Wz » 
: wi wx 
Energie: ot Cot a=" 1 pew ta 
Impuls: or ut +p = ex wh + ps, 
Zustandsgleichung : ee ee 
be eee 6 a RT, ox Rx Tx S 


ergeben die den Gleichungen (2b) bis (8b) entsprechenden Lésungen in thermochemischer 
Schreibweise fiir praktisch nichtideale Gase, wenn a? = 0, p,/0; = 6, g R, T, ist 


[ Ce ae Rz | 

= = of wt +1 (1 {- oa at/w}] (= —2) == 

x x x 
aR = R Lee i a faa 2 is ? (2) 

1 pice es) (— ai /w? +- 1) 

Cp ox O71 
2 

Px wy __ wx 3 
Pr 4 + Oy aa (. ae ( ) 
eae Ms. (4) 
0; Wx 

Tx Ry, ws we Wx 
ee Se Se eee} CO aoe Set 5 
TS Rx wy | a 1 az Wy ( ) 
sestcetiels (6) 
a? o,R,T, 

w;/ ay, AS oR, T, wr, (7) 
w/a? ox Ra Tx wi” 


1 BE, Sanger, Interavia, 3 (1948), 5.505; Brennstoffchemie, 32 (1951), 80137... 

E. cee I. Bredt, P. Goercke, C. R. Acad. Sc. 230 (1950), S. 949; Z. angew. Math. Phys. 1 (1950), 
S. 209; NACA, T. M. 1805, 1951; Z. physik. Chem., 199 (1952), S. 46; Ann. Télécomm. 6 (1951), S. 250. 

E. Sanger, Z. Naturforschg., 5a (1950); 5. 467; Z. phys. Chem. 197 (1951), S. 265; Z, Sprengtechnik 1, 
1952), 5S. 47. 
- ease Z. Naturforschg. 6a (1951), S. 302; C. R. Acad. Se. 235 (1952), 5. 460; Z. Sprengtechnik 1, 
(1952), S. 151; Z. Naturforschg. 8a (1953), S. 204; Z. Weltraumfahrt 5 (1954). 


386 Singer: Zur Thermodynamik des Uberschall-GeradstoBdiffusors. Ingenieur-Archiy 


bzw. fiir halbideale Gase, wenn also R,= Ry= ¢px— Cyx = Cpy— v1 = Cpox — vox = “poi vor 


rie, WHOL (oy Ege 
Pore, ; 2 .— | 
14 we at/w? ( s 


Hy X— 1 : <0) x ; (2a) 


2 


Wx be 1 £ \ 

— =— 4 — a?/w? + i 1+ 1 ; 
i a? /w? + 1) 

a 


Ww, iy 
usw. 

Verwendet man nun an Stelle des fiir die Spitzenwerte gebrauchten Adiabatenexponenten 
der Luft x = 1,67 die Werte ~, = %,; = 1,40 und x, baw. x), nach Abb. 2, so erhalt man den 
in Abb. 4 voll eingezeichneten Verlauf der GaszustandsgréBen nach dem Sto und nach 
Gleichgewichtseinstellung. 

Der zum Teil erhebliche Unterschied der vollen und gestrichelten Kurven zeigt zugleich, 
da®B die Vernachlassigung des wirklichen Luftverhaltens mit veranderlichen spezifischen 
Warmen bei héheren Flug-Machzahlen auch keine ganz grobe Naherung mehr darstellen wiirde, 
sondern zu vollig falschen Ergebnissen fiihren kann. 

Die kleine Nebenfigur zeigt die zeitliche Feinstruktur des VerdichtungsstoBes hinsichtlich 
der Temperatur bei der Flug-Machzahl 10 in der Stratosphare. Wahrend der Stof selbst eine 
Temperatiurspitze von T= 6950 °K aufweist, sinkt diese Temperatur nach Einspielen aller 
Gleichgewichte im Gas auf weniger als die Halfte, namlich T, = 2890° K, also etwa auf den 
Schmelzpunkt von Molybdan. Der Vorgang ist ganz analog den bekannten Verhaltnissen 
im Detonationsstof. 

Durch die Abkithlung sinkt die Strémungsgeschwindigkeit auch im ahnlichen Verhaltnis 
ab, wahrend Druck und Strémungs-Machzahl durch die Gleichgewichtseinstellung weniger 
berithrt werden. Die Energieverluste durch abdiffundierende Elektronen und Photonen aus 
dem Zustand T,, sind dabei vernachlassigt, die mit diesem Zustand verbundene hohe Ioni- 
sation wird als rechtzeitig riickbildbar angenommen. 


4, Isentropische Unterschall-Anstrémung. Bei der isentropischen Unterschall-Anstrémung 
nach dem VerdichtungsstoB zwischen den Querschnitten F, und dem Einlaufquerschnitt F,, 
gelten folgende thermochemische Strémungsbeziehungen fiir 


Kontinuitat: Fy 02 We = Fy Oje Wips 
E Oe we, ale T __ wi, 
nergie: a, 1 “ponte te tr poly 
ie} ARC) 
] c\ Cu 1? 
Isentrope: a eg 
Py Oy 
. ; Px is Py 
Zustandsgleichung: = ap 


ox Rx Tee oy Ry Ty, 


Daraus folgt die leider nicht explizite Bestimmungsgleichung fiir p,,/p, in thermochemischer 
Schreibweise zu 


he (pulps) — Ios (9) 
F, / ? e asi) ae 
Cc 5 [ery , my : 
ye ae POY (a8 ft) 1 ee ened o x1’ 
) a0 Gas (45/ Wx) Reece (Py|Px) 


Mit Hilfe ausreichender Zustandsdiagramme der Luft bis in den thermochemischen Strémungs- 
bereich und der Gleichungen (7), (5) und (2) kann man die Funktion (9) in der Form 


Pir iy, Ww) 

——— -——_ Jf. —, ae! / 

Px ie i “Gh (9) 
tabellieren oder in Diagrammform bringen, wie es in Abb. 6 in rudimentarer Form versucht 


wurde, so da die Zahlenwerte ohne umstindliche Rechnung jederzeit verfiigbar sind. 
Damit ergeben sich dann die weiteren ZustandsgriéBen im Querschnitt F,, 


OATS Es B, Falke 
"as ’ 


Wy. Fy, (10) 


ale 


Gir (Pe \ Neat? 
Ox Px ; (11) 
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Ox y/—1 
g — Rx Pi) cian (12) 
Te Whe Ws Z 
Cg peel 
wy / ay, = ry Py\ a Owl? 
ee ae et px Ge) 
‘in thermodynamischer Schreibweise fiir halbideale Gase 
_—1fx , 
F, ae ay 7 5 (PrrlPx) = up a (9a) 
PF, / “ a Hey) rT 
/ 2 x Wan eel ; ; 
Ete ek ea ss WON, ||] ees Se ihioy ee 
| eee A ge Sa ne ee } 
in gasdynamischer Schreibweise fiir ganzideale Gase 
: — 1/% 
F,, = 7 (P1//Px) fe % : ¢ (9b) 
F, a mere Lo a Poke 
| 1+ ees (a3/w%) cE —(Ppy/px) ™ 
Nach (13) wird in F,, die Schallgeschwindigkeit erreicht, wenn 
Pi = Wx Ox FP? | on 1/41 (14) 
Px ee Op Be, 


ist, womit wieder alle kritischen Zustandsgréf$en in F,, mittels (10) bis (13) folgen. 


Insbesondere ergibt sich aus den Gleichungen (14) und (9) eine Beziehung fiir das bei 
stationdrer Strémung iiberhaupt gréS’tmigliche F,/F,, in Abhangigkeit von der Anfangs- 


Machzahl w,/ax zu 


Fy, s Ve Ws 
F, so Oy ax 


ae 


is 06 67 a7 40 
Unterschall —Machzah! Mae 7074/7 
L00 O15 165 405 Shh 40 


Yberschall-Flugmachzah! 74/a7 


Abb. 5. GréBtmégliche Verhiltnisse von F,/Fy 
iiber der Unterschall- -Machzahl w,/a, oder 


w,/a, und tiber der Uberschall- Flug-Machzahl 

w, a,, mit Riicksicht auf das Erreichen der 

Schallgeschwindigkeit im Einlaufquerschnitt 
Fy 


Ox I/F] 
01,/6. = 2 Cpov/ Ov Ry, 2( 6, 7/—1) (15) 
w? a2 + 2 cpox/Ox Rx 
Pifflas — \fuaifftep 
Mines | Fh=0, Ye Hepy'=hb; Gps =Ho7'=1046 
[f hi|hi-O, Hay-Hepy! 13, Hoge ~Hop'-4337 ee) 
fe us Eh hi=0, 3 ta “Uant'Hox Hora 
‘ HOI, | »— 
hifi '=030, ” De) 
18 1 Le Ff, y rn 
a | | fff'=050; Fi Be 
7 = 
fi/hf'=0,60; 
16 
110 45)— 
id 
= 
105 
1 —— 
ih 
400- 0 
Abb. 6. Druckinderung Tre = Py! , [px isentropischer Strémung in Ab- 


hangigkeit vom Flachenverhaltnis F/F, , und der Ausgangs-Machzahl w, ./4 Ss 
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In gasdynamischer Schreibweise vereinfacht sich die letzte Beziehung zu 


x+1 
Up ail & x—1%\2 1) 
Fy, oS Wx (= 1 | Bop eae : (lez) 


Bei der Anfangs-Machzahl w/a, = 1 wird F,/F, also héchstens Eins. Dieser Fall tritt offenbar 
nur bei der Flug-Machzahl 1 ein. Aber auch schon bei der praktisch sehr wichtigen Anfangs- 
Machzahl w,/a, = 0,8; entsprechend etwa einer Unterschall-Flug- Machzahl w,/ ay = 0,8 baw. 
einer Uberschall-Flug-Machzahl 1,25 und dem gesamten dazwischenliegenden Schalldurch- 
gangsgebiet, kann die Eintrittsflache, wie Abb. 5 zeigt, des Luftzylinders F, praktisch kaum 
mehr gré$er werden, als die gebaute Eintrittsflache F',, des Diffusors. 

In Abb. 6 und 7 sind die Glei- 
jee chungen (9) und (10) zahlenmaSig 
Lie co aia mal 1 etic wal A dargestellt fiir den beim Schalldurch- 

| | | | gang wichtigen Bereich F,/Fy,< 1. 
ba | 7 7 iar L Saati Diese Diagramme gelten daher auch 
Oe H; Mzp/Qg_g=0 \ fir w,/a,< 1. 


08|\ 1. —— #ar=Hany! Hepp =tg= U4; Mep/[ Ue = 95 4 
» » » =95.. \ 
» » ” =Q7. 
y » ” =Q8. 

in ae | ae » ” y =99.N 


» ” ” =40 \ : 
Hp =He1'= 15, Hog =Hoy'=133, Up/ep= 40. OK K\ 

95 rey dery'= 14; Haar=Hor'~%246; M/~40 KX 
SS 


05 


O4}— 


Q3}- 


02) 
Of Me 
Alp 
0 Of 02 G3 04 OS G6 O7 08 G9 40 
Abb. 7. Geschwindigkeitsinderung w,,//w, isentroper Strémung in Ab- Abb. 8. Isentropische Unterschall- 
hangigkeit vom Flichenverhiltnis F,/F,- und der Ausgangs-Machzahl Wee / Oxe Anstrémung. 


5. Schub der isentropischen Unterschall-Anstrémung. Bei der Herleitung der Strémungs- 
beziehungen isentropischer Unterschall-Anstrémung wurde der Impulssatz nicht beniitzt, so 
da® dieser zur Berechnung des Schubes dieser Strémung noch zur Verfiigung steht. 

In Abb. 8a und b ist die hier betrachtete Strémung durch strichliert angedeutete, rotations- 
symmetrische Kafige abgegrenzt, fiir die beiden Falle F,/F,, <1 und F,/F,,>1. Nach dem 
Impulssatz muff der aus dem Kafig austretende Impuls, vermindert um den eintretenden 
Impuls, gleich sein der geometrischen Summe aller auf die umgrenzte Gasmasse wirkenden 
Driicke. 

Die unbekannte Drucksumme auf die gekriimmten Seitenwande des Kafigs, die aus 


Symmetriegriinden mit der Diffusorachse zusammenfallen mu8, wollen wir mit P,,, bezeichnen. 
Wir kénnen uns diese Seitenwande, die von keiner Stromlinie durchsetzt werden, zunachst 
materialisiert denken, d. h. der Diffusor erscheint in der gezeichneten Form nach vorne durch 
eine diinne Blechschale so verlingert, daB sein Einlaufquerschnitt genau gleich F, wird. 
Es gilt dann der Impulssatz in der Form 


Oy Fy, w; — Ox F, w= px Fy+ Pie Dg gas 


Da der tatsachliche Schub P,,, auf den Kafigmantel jedoch nur die geometrische Summe 
der Uberdriicke oder Unterdriicke gegeniiber dem Druck p, der ungestérten Umgebungs- 
atmosphare ist, wird weiterhin 


Pry = Pry — pi (Fy — Fi) = Fy (py — py + Oy wt) — F, (pe— py + Oe wz). (16) 
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Mit Hilfe der Beziehungen (3), (4), (8), (10), (11) und P;/01 w? = (1/0) (a2/w?) 1aBt sich diese 
Gleichung noch in die Form bringen 

Jey 1 six , D6 Diy 
as Frat ~ oy (att) [pe Crue sy 1) — rey — )] 4+ 2 (Se — 1) 


| w, \ wx 
=e FRG) i 1) (1 a ae ales Maes =| 
1% Px W, | Px FY} \ px Wy; es 
die den Schub der isentropischen Anstrémung als Funktion nur der vorgegebenen GréBen 
w,/a, und F,/F,, angibt, wenn man fiir die Werte wx/w, baw. py/px noch die Beziehungen (2) 
baw. (9) benitzt. 
Fir F,/F, > 1 wird dieser Schubanteil positiv, im andern Fall negativ. Dieser andere 
Fall, F,/F\, > 1 ist also nicht nur nach Abb.5 im Schalldurchgangsbereich nur in sehr he- 
schranktem Umfang realisierbar, er ist auch fiir die Schuberzeugung ungiinstig insofern, als 
der Hauptschub des Diffusors durch jenen der Anstrémung vermindert wird. Wir werden noch 
ae daB dieser Fall allerdings im héheren Uberschallbereich die einzige Regelméglichkeit 
offen laBt. 


In thermodynamischer Schreibweise lautet (17) 
Py. Hy, 2 Py Wx Fy \?/py\— Ve wx 
ss = = fatint\le 1 Yee tse Ie a Rede Gt SB pct 
Q Fy wi FB, E ( i/ y be )+ wy) px * \Fy} \px Wy Me Cee 


wobei fiir w,/w, baw. p,,/p, nun (2a) baw. (9a) zu beniitzen sind. SchlieBlich in gasdynamischer 
Schreibweise etwas ausfihrlicher nach Einfithrung von (2b) 


(17) 


9 


[Pre [Peal De 2 Py \ 
ed ere aaa 


2 5 
0, F, wi F; 25 


1 eee ee CN 20 2) eg 
be TS Be 108} VP) \ ps ; J 


Fir F,,/F, = 1 wird der Schub P,,, naturgema8 gleich Null. 


(17b) 


6. Regelung des Diffusors. Bei der Ableitung des Schubes der isentropischen Unterschall- 
Anstrémung war der tatsachliche Diffusor nach vorne durch diinne Blechwainde von der Form 
der Seitenwande des Impulskafigs verlangert gedacht, auf die sich die Strémung abstiitzen 
und den errechneten Schub iibertragen kann. 

Tatsachlich bestehen diese gedachten, materiellen Zusatzwande aber nicht, und der Schub 
mu in Wahrheit an den AuSenwanden des wirklichen Diffusors iibertragen werden, worauf 
die Méglichkeit beruht, den Uberschall-GeradstoSdiffusor bei fester geometrischer Gestalt in 
gewissem Umfang regeln zu kénnen. Daher muf die auBen um den Kafigmantel strémende 
Luft jene Driicke auf die innen langs des Kafigmantels strémende Luft ausiiben, die wir bisher 
dem materiellen Kafigmantel zugeschrieben hatten. Die Masse des auBeren Luftringes um den 
Mantel erfahrt dadurch eine auBere Kraft, die im Fall a von Abb. 8 nach vorne, im Fall 6 nach 
rickwarts in die allgemeine Strémungsrichtung weist, und deren Lage mit der Diffusorachse 
zusammenfallt. 

Da diese AuSenumstrémung, abgesehen von Strémungsverlusten, jedoch insgesamt keiner- 
lei axiale iuBere Krafte erfahren kann, mu der ringférmige Strémungskérper im weiteren 
Strémungsverlauf eine genau so grobe, entgegengesetzte axiale Kraft durch Driicke der 
wirklichen DiffusorauBenwande auf ihn erfahren. Im Fall a miissen diese offenbar Unter- 
driicke, im Fall 6 Uberdriicke auf die Diffusornase scin, die wirkliche DiffusorauBenwand 
erfahrt im erstern Fall einen Sog nach vorne, im zweiten Fall einen Formwiderstand nach 
riickwarts, abgesehen von den eigentlichen Luftwiderstandsvorgangen, die uns hier nicht 
beschaftigen. 

Dieser Ubertragungsmechanismus des Schubes der isentropischen Unterschall-Anstrémung 
auf die AuBenwande des Diffusors setzt natirlich voraus, daB geeignete, derartige Wande 
iiberhaupt vorhanden sind. 

Wenn der eigentliche Diffusor beispielsweise lediglich aus einem unendlich diinnwandigen, 
zylindrischen Rohr vom Querschnitt F',, bestiinde, ware dies offenbar nicht der Fall. 

Da dann der Schub nur durch Uber- oder Unter-Driicke auf die scharfe Kante tibertragen 
werden kann, die physikalisch itberhaupt méglichen Uber- oder Unterdriicke aber durchaus 
begrenzt sind, kann man bald eine kleinste Grenzdicke der scharfen Einlaufkante angeben, 
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unterhalb der der Schub jedenfalls kleiner bleibt, als er zu den vorgegebenen Strémungsver- 
haltnissen infolge der stationiren Strémungsmechanik unabanderlich gehort. 

Die notwendigen Strémungskrifte kénnen also infolge fehlender, achsparalleler Wand- 
druckkomponenten nicht hergestellt werden. Da anderseits durch die Vorkehrungen am riick- 
wartigen Ende des Diffusors jedes F,/F',, insbesondere jede beliebige Verengung F/ Fy aa 
des Querschnittes F', erzwungen werden kann, bleibt in diesem Fall als Ausweg nur eine grund- 
sitzlich andere Strémungsform iibrig, z. B. Riickstrémung eines Teiles der Luftmassen in die 
Bewegungsrichtung des Diffusors, oder der Ubergang zu pulsierender Strémung. 

Dabei beginnt die Gassiule mit ihrer Eigenfrequenz zu schwingen und stellt das Impuls- 
gleichgewicht durch nach vorne oder nach riickwarts weglaufende Impulswellen her!, ohne dafB 
es zu Schub auf den Diffusor kommt, dieser wirkt mehr als Orgelpfeife denn als Triebwerks- 
organ. 

nn physikalischen Grundlagen des Vorganges sind sehr ahnlich denen der bekannten 
thermischen Verstopfung beheizter Strémung, doch stellt diese Art der Strémungsverstopfung 
einen viel allgemeineren Fall dar. 

Praktisch werden nun die Verhiltnisse selten ganz so ungiinstig liegen, der Eintrittswulst 
des Diffusors wird immer eine endliche Dicke d haben, so daB er eine Stirnflache von der Gréfe 
(F,, — F,,) den Unter- oder -Uber-Driicken darbietet, welch letztere im Grenzfall absoluten 
Vakuums den héchstméglichen Sog 


Slee Pen Dee ian 
PTTL = (ate) : (tee) 


beziehungsweise im Grenzfall vollen Aufstaues nach dem VerdichtungsstoB den héchst- 
mioglichen Druckwiderstand 


D if Felis Py 
pee eS ee Pe che TiS hs of x 
aah woh oe ry Pe 1) (* 1 (19b) 
iibertragen kénnen, wenn man sich hier — angesichts des approximativen Charakters der 


ganzen Betrachtung — mit gasdynamischer Schreibweise begnigt. 

Der nach (17b) tatsachlich zu tibertragende Schub der isentropischen Unterschall-An- 
strémung lautet, bezogen auf die Bedingungen vor dem StoB und bezogen auf den im Grenz- 
fall auf die Diffusorwande itibertragharen Schub, 

Py 


S 


2% Hepes 2 a ee 
(Pr [p21 + | eta = 1) pyle + (F/R) (eis) ( 


<1l= 


x— | 


w/a? + i — # wi/a? - B/E, 


(Fy-[Fy —1) ; 
(20) 
fir F,/F,, < 1, baw. 
IBS = 
WY s1l= 
pe Saatin ee — fx 47 — 1 
(Py/px — 1) += Ty |‘ wy /at—- 1) pylpe+ (Fy/Fy)* (py-/px) ( wt /at + i —x w/a? - F,/F, 


(F,,/Fy — 1) ( P2/Pi— 1) 
(21b) 
fir F\/F;,>1. Diese beiden Gleichungen geben die Mindestwerte von (Fy,/Fy-)min in Ab- 
hangigkeit von F,/F,, und w,/a, an. Abb. 9 zeigt eine numerische Auswertung fiir x = 1,4. 
Wir verfolgen die Regelbarkeit des Luftdurchsatzes an einem Diffusor mit der allerdings 
geringen Nasendicke d von 4,25% des Durchmessers des gebauten Einlaufquerschnittes F',,, 

also mit F’,,/F',= 1,2, der in den Nebenfiguren der Abb. 9 skizziert ist. 
Im inkompressiblen Unterschall-Flug-Machzahlbereich bestehen offenbar praktisch keiner- 

lei Einschrankungen der Regelbarkeit. 

Aber schon nach Uberschreiten der Flug-Machzahl 0,5 beginnen diese merkbar zu werden, 
insofern, als bei zu groBen F, der Durchflu8 beschrankt wird durch Erreichen der Schall- 
geschwindigkeit in F\,, dagegen bei zu kleinen F, die notwendige Sogkraft an der Diffusornase 


1 J. Ackeret u. N. Rott, Schweizer Bauzeitung, 67, (1949), S. 40, 58. 
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wegen der auf absolutes Vakuum begrenzten Unterdriicke an derselben nicht mehr ubertragen 
werden kann. 

Bei Anniherung der Fluggeschwindigkeit an die Schallgeschwindigkeit der ungestérten 
Atmosphire (w,/a;—> 1) kann schlieBlich F, nicht mehr gréBer als F,, werden, wogegen fiir 
kleinere F’, noch ein gewisser, beschrankter Spielraum bleibt. 

: Nach Uberschreiten der Flug-Machzahl 1 6ffnet sich der Spielraum fir groBe F, nach oben 

_ wieder rasch, weil der VerdichtungsstoB fiir geniigend tiefe Unterschall-Strémungs- Machzahlen 

vor dem Diffusoreintritt sorgt, doch sté®t man mit zunehmender Flug-Machzahl bald an die 

_ Grenze, wo der Staudruck.auf die Nasenflache nicht mehr zur Ubertragung des Schubes aus- 

reicht. Dadurch wird fiir die betrachtete Nasendicke eine Grenze fir groBe F, schon zu 

1,28 F, gezogen, die weiterhin fast unabhangig von der Flug-Machzahl bleibt, da der Schub 
der Anstrémung ahnlich wie der Staudruck wiichst. 
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Abb. 9. Maximale Regelbereiche des Uberschall-GeradstofSdiffusors, abhangig von der Nasendicke d/VFy = (VR, [Py — Tiler 


Dagegen schrumpft nach Uberschreiten der Flug-Machzahl 1 der Bereich kleiner F, sehr 
rasch vollig gegen Null, da hier dem anwachsenden Anstrémschub der konstante Vakuumsog 
gegeniibersteht. Bei héheren Flug-Machzablen kann daher F, praktisch nie kleiner als F,, 
werden. 

Aus Abb. 9 1aBt sich ablesen, in welchem Umfang sich diese engen Regelbereiche durch 
dickere Nasen erweitern lassen, insbesondere sieht man, da8 Fy/F,, >1,5, d. h. Nasendicken 
tiber 12% des gebauten Einlaufdurchmessers, sinnlos werden, da sie im Bereich groBer F,/F,, 
wegen Erreichens der Schallgeschwindigkeit in F, unwirksam sind, im Bereich kleiner F,/F,, 
praktisch kaum eine Erweiterung des Regelbereiches bewirken, und auBerdem sehr ungiinstig 
auf den hier nicht betrachteten reinen Luftwiderstand des Diffusors wirken, besonders wenn 
F,/F,, = 1 ist. Die tatsachlich ungleichmafige Druckverteilung auf die Diffusornase wird die 
zu 12% gefundene Grenze der Nasendicke vielleicht noch etwas weiter hinausschieben. 

Die tatsdchlichen Regelgrenzen des Luftdurchsatzes werden aus diesem Grunde aber 
sogar noch etwas enger sein, als Abb. 9 angibt, da die wirkliche Strémung nicht die voraus- 
gesetzten extremen Driicke auf der ganzen Stirnflache der Diffusornase herstellen wird. 


Pa 
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Eine Luftdurchsatzregelung des Diffusors fester geometrischer Gestalt wird daher in den 
angegebenen Grenzen im Fluggeschwindigkeitsbereich des Schalldurchganges zwischen 
F,/Fy <1 und F,/F,, = 1 méglich sein, im Bereich hoherer Uberschall-Machzahlen dagegen 
zwischen F,/F,, = 1 und F,/F,,> 1. 

Beim Staustrahltriebwerk mit fester Triebwerksgeometrie kann diese Variation des Luft- 
durchsatzes bekanntlich! durch Anderung der Aufheizungstemperatur der Luft erzielt werden, 
indem bei sinkender Aufheizung der Luftdurchsatz steigt und der Schub sinkt und umgekehrt. 

Man wird beim Schalldurchgang daher mit hoher Aufheizung zu fahren haben, bei héherer 
Fluggeschwindigkeit dagegen mit niedrigerer Aufheizung, was auch den flugmechanischen und 
thermischen Anforderungen entspricht. 

Die in Abb. 9 dargestellten Begrenzungen engen also den Bereich der automatischen Quer- 
schnittsanpassung an den Durchflu® beim Uberschall-Fangdiffusor sehr merkbar ein, dieser 
ist weit weniger elastisch, als der entsprechende Unterschall-Fangdiffusor. 


7. Adiabatische Einstrémung. Wahrend die Unterschall-Anstrémung vor dem Diffusor, 
zwischen F, und F,, ohne weiteres als isentropisch angenommen werden durfte, ist dies fiir 
die eigentliche Diffusordurchstrémung zwischen F', und F, wegen der dort merkbaren Grenz- 
schichtverluste an den festen Wanden erfahrungsgema8 nicht mehr der Fall, sie bleibt indes 
naherungsweise adiabatisch. 

Der Druckdiffusorteil hat einen von Eins verschiedenen ,,Diffusorwirkungsgrad*: yp, 
worunter man traditionsgem48 das Verhaltnis der im Druckdiffusorteil wirklich erzeugten 
Druckerhéhung zu der bei isentropischer Strémung theoretisch erzeugbaren Druckerhéhung 
versteht, und der empirisch ermittelt wird. Dieser traditionelle Wirkungsgrad stellt bei kom- 
pressibler Strémung allerdings weder ein Energieverhaltnis, noch im allgemeinen ein Schub- 
verhaltnis dar. 

Die tatsachlich erreichten Diffusorwirkungsgrade np bei den iiblichen Durchstrémge- 
schwindigkeiten hangen hauptsachlich von der Diffusorbauart ab und gehen praktisch von 
Werten iiber 7p = 0,95 bei groBen, sorgfaltig geformten Diffusoren, bis herunter zu den 
Werten der Borda-Carnotschen Strahlausbreitung bei vollstandiger Strémungsablésung 
Hoe= 2/(1 + F,/F,), vorausgesetzt, daB der Diffusor lang genug ist, um véllige Strahldurch- 
mischung zu erlauben, andernfalls kann der Diffusorwirkungsgrad auch noch kleiner werden. 

Er hangt aber auBerdem von der Durchstrém- 
geschwindigkeit selbst ab, wie man schon daran er- 
kennt, daB bei w, = 0 jedenfalls der volle Staudruck 
erreicht werden muf, d.h. bei F,/F, =0 muB yp 
gegen Eins gehen, und ebenso wachst er mit an- 
wachsendem Diffusorverhaltnis F,,/F,, indem er dem 
Wirkungsgrad der schublosen, Borda-Carnotschen 
Strémung ausweicht, die ja bei Diffusorverhaltnissen 
uber F,/F, = 0,75 selbst schon die Wirkungsgrade 
guter Diffusoren iibersteigt. 

Abb. 10 zeigt dieses Verhalten eines groBen, gut- 
geformten Diffusors von etwa konstantem Offnungs- 
winkel im inkompressiblen Strémungshereich — 
der mit anwachsenden Uberschallgeschwindigkeiten 

“jg wieder gut zutrifft — itber dem gebauten Diffusor- 

0 Os 40 verhaltnis F',,/F,. Die Angabe eines 90 bis 95 %igen 
ane Bs otlent Aes Web aaetersg ial D eee a Diffusorwirkungsgrades bezieht sich naturgemaB auf 
verhaltnis F.)/F,, im Vergleich zurschublosen Borda, wittlere F',/F,, in deren Grenzbereichen (Null bzw. 

Carnotschen Verdichtungsstrémung 7pp. Eins) steigt der Wirkungsgrad notwendig gegen Eins, 

eimerseits weil der Staudruck fir F,,/F,—> 0 voll 
erreicht wird, anderseits weil, wie gesagt, der Wirkungsgrad der schubfreien Verdichtungs- 
strémung mit voller Strémungsablisung hei F,,/F,—> 1 nicht unterschritten werden kann. 

Diese Komplikation riihrt yon der unzweckmafBigen traditionellen Definition des Diffusor- 
Nees peace ae eo ene pees Sees ee des Diffusors nicht unter 
Laken hevat Wisk oie TSM oe aoe or a-Carnotschen Strémung sinken 
e iebselement wird natiirlich in Wahrheit 
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‘ E, Sanger, Schw. Arch. f. angew. Wiss. u. Technik, 16 (1950), S. 341. 
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nicht erst bei der Drucksteigerung Null gleich Null, sondern — auf die Verhaltnisse in Fi, 
bezogen — schon dann, wenn der Schub auf die Diffusorwinde infolge der Strémungsver- 
zogerung zwischen F, und F,, also die Diffusorwirkung verschwindet. Das ist aber bereits 
beim Erreichen der schublosen Borda-Carnot-Strémung der Fall, die selbst noch erhebliche 
Drucksteigerungen aufweist. Man wird daher eher damit rechnen diirfen, daB ein Ausdruck 
von der Art (yp — ypp)/(1 — ps) aber weitere F,,/F,-Bereiche konstant ist. 


Dain Abb. 10 die Ordinate 7p = 1 dem vollen, verlustlosen Diffusorschubgewinn der Innen- 
durchstrémung entspricht, erkennt man, da8 eine Wirkungsgradverminderung um einen 
bestimmten Bruchteil nur im Falle F,,/F,— 0 dem gleichen Bruchteil an Schubgewinn-Ver- 
minderung entspricht, mit anwachsendem F’,,/F, jedoch immer griéferen Schubverlusten. Dies 
gilt allerdings nur fiir jenen Schubanteil, der durch die Druckerhéhung infolge der Innen- 
durchstrémung allein bewirkt wird, z. B. nicht aber fiir den Anteil an Druckerhéhung infolge 
des VerdichtungsstoBes oder infolge der isentropischen Unterschallanstrémung. 


Dieses Wirkungsgradproblem verliert infolgedessen bei dem hier zu betrachtenden Uber- 
schalldiffusor, besonders bei héheren Machzahlen, sehr an praktischer Bedeutung durch den 
Umstand, dai der gréBte Teil der niitzlichen Drucksteigerung und Schuberzeugung schon 
durch den geraden Verdichtungssto8 und zum Teil auch durch die Unterschall-Anstrémung 


erzeugt wird, denen gegeniiber die Drucksteigerung der Einlaufstrémung selbst oft eine 
geringere Rolle spielt. 


Zur Berechnung der Einlaufstrémung mit Verlusten kann man folgende Beziehungen 
verwenden, wenn die Dissipation der entstehenden Turbulenz in Enthalpie unendlich schnell 


angenommen wird. An Stelle des isentropischen Enddruckes p; erhalt man nur den wirklichen 
Enddruck 


Pe = Py + HD (P; = Piz) 


Bei der nichtisentropischen jedoch adiabatischen Strémung zwischen F',, und F, stehen daher 
die finf Beziehungen zur Verfiigung, fiir 


Kontinuitat: Oy Fy, Wye S Os Fy Wa. 
; we, we 
Energie: 5 2 + Cpoy ly = re = Cneg lo > 
: = Rite eae P2 
Zustandsgleichung: RT, ents” 
. : [Boe ONG 
Diffusorwirkungsgrad: YpD=—> ; 
(DRE TY 
Ik —lfoy+6 

/ ie (peeps) 12 

Mie RROD ESC Hee) Cite QUE i AU Ge ne ; 
a / Ve en tee nl Poe 
Pe ee aes a pp ae 
Ce Ry 1// Vy R, panne! 


Aus den beiden letzten Beziehungen folgt ohne weiteres der Zusammenhang des p,/py mit 
F,/F,, und. w,,/a,, zu 


( Po/ Py —— a 


By = ie Gib 
; py» —1\ oy79 

|/ IPs; 2 pose (at/wy) }1— ee (14 ae aa 
oO, Ry 1a Cp yy \ "ip | 


Mit Hilfe ausreichender Zustandsdiagramme der Luft bis in den thermochemischen Strémungs- 
bereich und mit Hilfe der Gleichungen (13), (9), (7), (5) und (2) kann man die Funktion (22) 
in der Form 


D,, F, Rae uw \ 99/ 
Dae (Fe Bee io} (22’) 
af 


tabellieren oder in Diagrammform bringen, so da die Zahlenwerte cin fitr alle Mal verfiighar 
sind. Gleichung (22’) ist offenbar lediglich cine Erweiterung von (9'). 
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Wenn damit p,/p,, bekannt ist, folgen aus den itbrigen der fiinf Beziehungen die restlichen 
ZustandsgréBen des Gases in Fy zu 


i 2 \2 ) 0 j 2 
Les Cpo2 Fy ry ai a 2 Cpoir Air Cpor Fy Pz Air (23) 
Py. = ) 

ee oy R, Fy py w}, Oy, Ry, wi, Oy Ry Eye Pye wi, 

ARs Ry, Fy py wy (24) 
ee a P) 

Pe R, Fy, py Wye 

0 Fw, (25) 
aS 2 ee ean 

Ow F, w, 

W,/d, | / ye By, Py Wy (26) 
Wyr/ Ay, 0, F, p, Wy 


Eine Komplikation entsteht erst wieder, wenn man in F’, so hohe Massenfliisse fordert, 
daB in F,, die Schallgeschwindigkeit erreicht wird. Damit ist die bekannte Grenze fiir den 
Flu8 0, w, erreicht und der Druckdiffusor mite weiterhin als divergente Uberschalldiise 
wirken, also mit Uberschallgeschwindigkeit bei fallendem Druck durchstrémt werden. Da 


bei den praktischen Anwendungen jedoch meist p, > p, ist, wird die Strémung im Diffusor | 


mittels eines oder mehrerer VerdichtungsstéBe wieder zu Unterschallgeschwindigkeit zuriick- 
kehren. 

Die Lage eines einzelnen VerdichtungsstofSes ist dadurch bestimmt, da8 die Drucksteigerung 
durch ihn und durch die anschlieBende Verdichtungsstrémung in dem dann wieder als Diffusor 
wirkenden, divergenten Strémungskanal auf den héchstméglichen Druck p, zuriickfithren. 

Eine einfache Nachrechnung zeigt, da p, um so héher wird, je weiter vorne der Verdich- 
tungssto8 im Diffusor sitzt, und am héchsten, wenn der dann verschwindende Verdichtungs- 
stoB in F,, saBe. Je weiter der StoB in den Diffusor hineinwandert, desto geringer wird po, 
der Diffusorwirkungsgrad sinkt in diesem Fall auBerordentlich ab. 

Die adiabatische Einlaufstrémung mit VerdichtungsstéBen im Innern des divergierenden 
Diffusors bietet daher kein praktisches Interesse und wird tunlichst zu vermeiden sein. 


8. Schub des Diffusors. Der Gesamtschub P,,. der Einlaufstrémung, d. h. die geometrische 
Summe aller Uber- oder Unter-Driicke gegeniiber dem ungestiérten, atmosphiarischen Um- 
gebungsdruck p, auf die festen Innenwinde des Diffusors, folgt ohne weiteres aus dem Impuls- 
satz 

Qo Fy wy — Oy Fy wh = py Fy, + Pry — Pp Fy 
Zu 
Pug = Pry — Pi (Fe — Fy) = F, (Po — Py + Oo w3) — Fy (Pv — Pa + Ow?) 
oder, in Analogie zu Abschnitt 5 


Py ae 1 I F, L F, v3 OSI Rit i, C1 Py, wy 
0, FF, wt guy FF, 0, F, w} Qo, wy #F, a1 Fy wt 
] (a2/w2) Fs ( 1 1 Wy, Ux [| Wy (27) 
= — (A4/ we) |= (Mo My + Hy — — (21,7, 1) - a il 
6, Aiea F,* 2 Vx? Uy ) F, (Fr. Mer ) ae We Ww, \ wy ? 


worin noch fiir w,,/wy, wx/w,, w,/w,, beziehungsweise die Ausdriicke (10), (2) und (23°) zu 
setzen sind. 

In dem auch praktisch wichtigen Sonderfall, daB der Wirkungsgrad "p der adiabatischen 
Einlaufstrémung sich eins nahert, also auch dieser Strémungsteil isentropisch wird, braucht 
bei der Schubberechnung kein Unterschied zwischen den beiden Strémungsteilen zwischen F, 
und Ff’, einerseits und F,, und F, anderseits gemacht zu werden. Wir setzen in diesem Fall 
naherungsweise F,,/F', = w/w, =. =1 und erhalten eine vereinfachte Form der Glei- 
chung (27), die der Gleichung (17) entspricht: 


Py, = 1 2 Jays F, Wx | Ws 
0, F, w? _ O1 (at/w?) F, (7p x? Te 1) (7.4 1) w, ee 1 P) (28) 


die in diesem Fall den gesamten Schub, einschlieBlich jenem der Einlaufstrémung umfaBt. 
Nebenbei bemerkt, erhalt man durch Nullsetzen dieser Gleichung eine Beziehung fiir die 
Borda-Carnotsche Verdichtungsstré mung in nichtidealen Gasen. 
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Der tblicherweise auf die Hauptspantfliche F, des Diffusors bezogene Beiwert des ge- 
samten Diffusorschubes wird mit (28) 


2 P eee F, ee Peery nape 
ard = a MD (ta —Y) —pea— 0) + 2p SS —1). es) 
Bei F,/F, > 0 wird der Durchflu8 null und 
2 
Cs = — (at /wi) (7%, %,1 — 1). (30) 


Oo; 


Dies ist also der einzige Fall, wo Schub und Drucksteigerung (p, — p,) im Diffusor einander 
streng proportional sind, P,. = F, (pz — p,), also wenn (p, — p,) proportional yp ist, wie 
im Unterschallbereich, auch P und 7p proportional werden. 

Dieser Schub entspricht genau dem Schub auf den Brennkammerboden einer Rakete. 
Sein Aufbau mu8 natirlich dadurch gesichert sein, daB die vorausgesetzte, geordnete Diffusor- 
strémung vor dem Zusammenbruch bewahrt wird, z.B. durch eine den Diffusor in F, abschlie- 
Bende Membran, auf die dann genau derselbe Gegenschub wirkt, oder durch die tblichen 
Triebwerkseinrichtungen, wie Aufheizung der Luft und Absenkung des Druckes bis zum 
AuSendruck oder wenigstens bis zum kritischen Druck. 

Der Fall des Aufstaues auf die Geschwindigkeit w, 0 ist zugleich der einzige, wo die 
Drucksteigerung p,/p, =7)x.'7x1 in geschlossener Form anschreibbar ist. Mit Hilfe von (3) 
und (22) fir F,,/F,—0 wird 


Ry cpox Ox Rx sty en = | 4 apes we 
Opa ( 95 Gee wz ax 2 il | 1+ O74 (w?/a®) 1 i 

(31) 
in thermochemischer Schreibweise, wobei fiir w,/w, Gleichung (2) zu beniitzen ist, fiir w;/a; 
die Gleichung (7). 


Im thermodynamischen Bereich, also fir Flug-Machzahlen w/a, <6, vereinfacht sich die 
Beziehung zu 


Be. Pa Pe 
Pi Px P 


See | no | 


Ming 2 
Po __ [[¥oo —1 ox Leper |e Ny | eee Beal |e 
eS [' : tte Hon  %gx—l (1 ; 2 Hox Bele %e3 | +s (uried) W; 
(31a) 


wobei fiir w,/w, Gleichung (2a) zu beniitzen ist, fiir w;/a; Gleichung (7a). 
Endlich vereinfacht sich (31) im Seed aren enon Bereich, also fiir Flug-Machzahlen 
w,/a, <2, zu 


= Man Me . ae + 4D 


Ww 
T ] 


acaiess: ie | nacace te) 


oder wenn man hier (2b) und (7b) gleich einsetzt, 


Po == My, Toys he -+- Wp 


D, Dee 21a O17 
a see =) 31) 
P, 2 x 1 hE tae (nt ) ity) 


Setzt man schlieBlich auch noch yp = 1, so folgt die bekannte Lord Rayleighsche Beziehung 
(3’b) 


[ Beate Ma as ees 
- w/a —1 7 « i 
Om ee is 2 Mwai aot 31k 
r= Man Me = \ae eae | (Atle ne ( Dp) 
1 ON Tie eis § 
lei | 


mit ihrem Analogon im Unterschallbereich: 


fre x—1 wt 
Sele ( += 2 ett) : Cri 
1 
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Fiir w,/a, = 1 geben die beiden letzten Gleichungen natiirlich denselben Wert. Bildet 
man fir den nach (30) mit (31’’’b) berechneten Schubbeiwert die Ableitung nach dex Flug- 
Machzahl, dc,/0(w,/a,) und setzt diese gleich null, so folgt an der Stelle w,/a, =0 ein Minimum 

der Héhe c,=1. * : 
Sa sate Vorgang mit den Uberschallbeiwerten aus (30) und (31 b) fihrt zu einem Maxi- 
mum von c, =1,85 bei w,/a4, = 00, d.h. auBerhalb des Giltigkeitsbereiches der gasdynamischen 
Beziehung. ee. 

Da mit unbegrenzt wachsender Machzahl die Isentropenexponenten aller wirklichen 
Stoffe physikalisch gegen x1 konvergieren, geht der tatsaichliche Grenzwert der Schub- 
beiwerte gegen den alten Newtonschen Wert c,=2, der aus der elementaren Mechanik ebenso 
folgt, wie aus (30) und (31’’b). 
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Abb. 11, Schubbeiwerte c, = P/q F, des Uberschall-GeradstoBdiffusors, ttber der Flug-Machzahl w,/a,. 


Den Gesamtverlauf der c, tber alle drei Strémungsbereiche fiir F,/F,—0 zeigt Abb.11. 
Im gasdynamischen Bereich stimmt er itberein mit der bekannten Auftragung der Beiwerte 
des Druckanstieges vor einem Staupunkt im Uberschallbereich. 

Insgesamt beginnen die Beiwerte mit der Flug-Machzahl Null bei c, 1 und streben mit 
der Machzah] Unendlich gegen c,=2. Bei der Flug-Machzahl 10 in Bodennihe bedeutet 
der dargestellte Beiwert einen Schub von iiber 1300 Tonnen je Quadratmeter Hauptspant- 
flache. 

In Abb. 11 sind auch die c,-Verlaufe fiir einige endliche F,/F, + 0 angegeben, deren Zahlen- 
werte sich leicht aus (29) ergeben, wenn man die 7, fiir die jeweils vorliegenden Verhaltnisse 
von w,/a, und F/F, aus Abb.6 abgreift, die a,, = px/p, und w,/w, aus Abb.4 und die w,/w, 
aus Abb.7. 

Der EinfluB des Diffusorwirkungsgrades 7p auf den Schubbeiwert ist am gréBten im Unter- 
schallbereich, um mit zunehmender Uberschall-Flug- Machzahl in dem Ma8’ abzusinken, wie 
die Druckerhéhung zunehmend durch den senkrechten VerdichtungsstoB und nur mehr zu 
immer geringeren Anteilen durch die verzégerte Strémung im konischen Diffusor stattfindet, 
auf welch letztere allein der Diffusorwirkungsgrad bezogen ist. 
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Die Gegeniiberstellung dieser durchwegs auBerordentlich hohen Schubbeiwerte mit den 
bekannten Widerstandsheiwerten von Uberschall-Flugkérpern. zeigt die hohe Wirksamkeit 
des Uberschall-GeradstoBdiffusors als Antriebsorgan. 

Die dargestellten hohen Diffusorschiithe werden immer einen wesentlichen Anteil des 
Schubes eines vollstandigen Luftstrahltriebwerkes bilden. 


Folgen hinter dem Diffusor noch Kompressoren — eine Anordnung, die in Form von 
_ Turbinenstrahltriebwerken und Turbinenstaustrahltriebwerken bis etwa zur Flug-Machzahl 2 
sinnvoll sein kann — so mag der Triebwerksschub noch wesentlich ither den Diffusorschub 
allein hinaussteigen. 
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Abb. 12. Drucksteigerung p,/p,; am Ende eines GeradstoB-Diffusors in Abhingigkeit von der Flug-Machzahl w,/a,, vom 
Diffusorverhaltnis F,/F, und vom Diffusorwirkungsgrad 7p. 


Wird der Luftdruck jedoch nach dem Diffusor nicht mehr weiter erhéht, sondern die Luft 
nurmehr erwarmt — eine Anordnung, die in Form der reinen Staustrahltriebwerke je nach 
der Aufgabenstellung ab Flug-Machzahlen von 1 bis 2 vorzuziehen ist —, dann folgt auf die 
meist zylindrische Heizkammer oft wieder cine Verengung des Strémungskanals zur Aus- 
strémdiise. Diese Verengung gibt naturgema8 einen Gegenschub, so daB der Gesamtschub 
des Staustrabltriebwerkes dann kleiner als der dargestellte Diffusorschub wird. 

Wenn man jedoch an den Diffusor nurmehr ein zylindrisches Rohr als Heizkammer und 
Ausstrémdiise anschlieBt — wie die Nebenfigur in Abb.11 dies andeutet — dann bleibt der 
Gesamtschub des Staustrabltriebwerkes gleich dem dargestellten Diffusorschub. 

Diese Anordnung gibt also die héchsten Schiibe, jedoch héheren Warmeaufwand je An- 
triebseinheit, als eine eingezogene Diise. Sie ist also besonders fiir Gerate geeignet, bei denen 
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der spezifische Treibstoffverbrauch nur eine sekundare Rolle spielt, z.B. fir Kurzzeitgerate 
oder bei Anwendung von Atomheizung. ‘ 

Diese Geriate erlauben auch wegen ihrer Sonderverhiltnisse am chesten die Beniitzung 
des gesamten behandelten Flug-Machzahlbereiches bis v/a = 10. 

Im Falle der rein zylindrischen Abstrémung mu8 offenbar die durchstrémende Luft sowelt 
aufgeheizt werden, da®B sie die Verstopfungsgeschwindigkeit (~ Schallgeschwindigkeit) er- 
reicht, das Feuergas wird dann mit Uberdruck gegenithber dem ruhenden Umgebungsdruck 
entlassen. Die Hohe der notwendigen Aufheizung wachst mit fallenden P/F, iitberschreitet 
etwa bei F,/F,~0,25 die mit chemischen Reaktionen erreichbaren Temperatursteigerungen. 
und geht mit F,/F,>0 gegen unendlich. Heizt man jedoch weniger, als bis zur thermischen 
Rohrverstopfung auf, bleibt also 
mit der Heizkammerstrémung im 
Unterschallbereich, so gleicht sich 
der Druck nach riickwarts mit dem 
AuBendruck aus, die Diffusorstré- 
mung bricht zusammen und die 
berechneten Diffusorschiibe kénnen 
nicht entstehen. 

Anderseits fiihrt die Aufheizung 
bis an die Verstopfungsgrenze auch 
in die Nahe instabiler Strémungs- 
verhaltnisse, die sehr sorgfaltige 
und gleichmaBige Heizungsvorkeh- 
rungen erfordern, um zu verhindern, 
daB das gesamte Strémungssystem 
ins Schwingen gerat. 

Diese Probleme der Ofenstré- 
mung fallen indes in den Bereich 
thermochemischer Strémung und 
werden daher in dieser thermo- 
dynamischen Studie nicht weiter 
behandelt!?. 


9. Zustandsgrofen am Diffusor- 
ende. Die wichtigste quantitative 
Eigenschaft des Uberschall-Gerad- 
stoBdiffusors ist zweifellos sein in 
malt nja; Abb.11 dargestellter Schub. 

0 7 2 ae, 5 6 7 Cee ae 10 Daneben kann es jedoch fir 
Abb. 13. Absolute Strémungsgeschwindigkeiten w, am Diffusorende (relativ manche Zwecke wiinschenswert sein, 
zur Schallgeschwindigkeit a, der Umgebungsluft) in Abhangigkeit yon der Flug- : ci ° * : 
Machzahl w,/a,, vom Diffusorverhiltnis F,/F, und yom Diffusorwirkungsgrad Ap: einen Uberblick iiber die Luftzu- 
standsgréBen am Diffusorende, also 
im Querschnitt 2, zu haben, etwa zur Berechnung der nach dem Diffusor folgenden Stré- 
mung, oder zur Beurteilung von mechanischen Festigkeitsverhaltnissen, Anblaseverhaltnissen 
der Feuerung, kritischen Strémungs-Machzahlen, Warmeibergang an die Diffusorwande usw. 

Abb.12 zeigt die Drucksteigerungen p,/p, am Ende des GeradstoB-Druckdiffusors in 
Abhangigkeit von Eintrittsverhaltnis F,/F,, von Flug-Machzahl w,/a, und vom Diffusor- 
wirkungsgrad 7p. Bei héheren Flug-Machzahlen erfolgt der Druckaufbau iiberwiegend durch 
den VerdichtungsstoB selbst und nur zum geringen Teil durch die Verdichtungsstrémung im 
Diffusor, so da dort der Einflu8 des Diffusor-Eintrittsverhaltnisses und des Diffusorwirkungs- 
grades sehr zuricktreten. 

Die zum Vergleich noch eingezeichnete Kurve des isentropischen Druckanstieges laBt das 
unmaBige Anwachsen der isentropischen Staudriicke bei héheren Machzahlen erkennen, dic 
keineswegs von einem entsprechenden Anwachsen des Diffusorschubes selbst begleitet sind 
und zu baulichen Schwierigkeiten fiihren miiBten. 

Abb.13 zeigt die absoluten Strémungsgeschwindigkeiten w, am Ende des Geradsto8- 
Druckdiffusors (relativ zur Schallgeschwindigkeit a, der Umgebungsluft) in Abhingigkeit 


10;,/2 


1 EB. Sénger, Schweizer Arch, f. angew. Wiss. u. Technik, 16 (1950), S. 341, 
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ae der Flug-Machzahl w,/a,, vom Eintrittsverhiltnis F,/F, und vom Diffusorwirkungs- 
rad 7p. 
4 Bemerkenswert erscheint, daf die absoluten Strémungsgeschwindigkeiten nach dem geraden 
an zunichst merkbar abnehmen, besonders bei groBen Eintrittsverhiltnissen, 
2 daB sie auch bei den praktisch wichtigen kleinen F,/F, nur sehr schwach zunehmen. 
ieser Umstand hat fiir den Bau der Heizungen der Uberschall-Staustrahlantriebe grund- 
legende Bedeutung. 
Interessant ist ferner, daS der Einflu8 der Kompressibilitit der Luft bei der Verdichtungs- 
strémung im Diffusorinnern, die das Abweichen der w,/w, baw. w/w, von F/F, bewirkt, in 


70;,/; 


Abb. 14. Strémungs-Machzahlen w,/a, am Diffusorende in Abhangigkeit von der Flug-Machzahl w,/a, und yom 
Diffusorverhaltnis F,/F,. 


der Nahe der Flug-Machzahl Eins am gréBten ist und sowohl nach kleineren, als auch nach 
gréBeren Flug-Machzahlen abklingt, und schlieBlich, da® sich auch ein von eins verschiedener 
Diffusorwirkungsgrad besonders in der Nahe der Flug-Machzahl Eins am stiarksten auf die 
Strémungsgeschwindigkeit am Diffusorende auswirkt. 

Abb.14 zeigt schlieBlich die Strémungs-Machzahlen w,/a, am Ende des GeradstoBdiffusors 
in Abhangigkeit von der Flug-Machzahl w,/a, und vom Eintrittsverhiltnis F,/F’,, die im Uber- 
schallbereich mit wachsender Flug-Machzahl monoton absinken. 


(Hingegangen am 19. Dezember 1953.) 


Anschrift des Verfassers: Dr.-Ing. Eugen Saénger, Forschungsinstitut fiir Physik der Strahlantriebe, E.V. 
Stuttgart-Flughafen. 
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Ein Beitrag zur Stabilitatstheorie der Wirbelstrafen unter 
Beriicksichtigung endlicher und zeitlich wachsender 
Wirbelkerndurchmesser.* 


Von U. Domm. 


1. Einleitung. Von besonderer theoretischer und experimenteller Wichtigkeit sind ein- 
reihige und doppelreihige WirbelstraBen. Bei der einreihigen WirbelstraBe sind die Wirbel 
auf einer Geraden in konstanten Abstaénden voneinander angeordnet. Die experimentelle 
Erfahru zeigt, da® einreihige WirbelstraBen instabil sind. Die hydrodynamische Theorie 
hat die Instabilitat der einreihigen WirbelstraBe fiir die ideale Flissigkeit durch Annahme 
kleiner Stérungen, die an den Wirbeln angebracht werden, nachgewiesen. 

Von denkbaren zweireihigen WirbelstraBen sind im besonderen die Anordnungen studiert 
worden, bei denen die Wirbel auf zwei parallelen Geraden vom Abstande h so angesetzt sind, 
daB ein Wirbel dem vorhergehenden im Abstande /, der Teilung, folgt. Die Zirkulation der 
Wirbel verschiedener StraBenseiten ist dem Betrage nach gleich, dem Vorzeichen nach ent- 
gegengesetzt. Fordert man, daB die WirbelstraBen sich als Ganzes parallel zu den Geraden 
translatorisch bewegen, dann verbleiben nur zwei Konfigurationstypen der Wirbel: Entweder 
stehen sich die Wirbel der verschiedenen StraBenseiten gegeniiber, oder die Wirbel stehen auf 
Liicke, d.h. die eine StraBenseite ist um die halbe Teilung gegen die andere Strafenseite 
verschoben. 

Im Experiment beobachtet man nur die auf Liicke stehende StraBe. Den ersten Hinweis 
findet man bei A. Mallock+, die ersten eingehenden Untersuchungen bei H. Bénard*. Die 
theoretische Deutung der WirbelstraBe mit auf Liicke stehenden Wirbeln verdankt man 
Th. v. Karman’. 

Nach v. Kdrmdan ist in idealer Flissigkeit die WirbelstraBe mit gegeniitberliegenden Wirbeln 
instabil, wahrend die WirbelstraBe mit auf Liicke stehenden Wirbeln unter der Bedingung 
h/l =0,281 stabil ist. Spater zeigte Bl. Dolaptschiew*, daB diese Stabilitatsbedingung bei 
Bericksichtigung von Stérungen zweiter Ordnung erhalten bleibt. 

Von Kdrmdn und Rubach® geben Mefergebnisse an, die die Theorie von v. Kdrmdn gut 
bestatigen. H. Bénard konnte die Stabilitatsbedingung von v. Kdrmdn experimentell nicht 
verifizieren. Er findet Werte von h/l in einem Bereich von 0,24---0,56. S.C. Hooker® zeigte, 
da8 man bei experimentellen Untersuchungen zur Bestimmung von A/l eine Korrektur 
anbringen mu, die aus dem realen Kernprofil der Wirbel resultiert. Die Verbreiterung 
von WirbelstraBen in einigem Abstand hinter dem Widerstandskérper ist nach Hooker 
nur ein scheinbarer Effekt. R. Wille’ findet, daB h/l eine Funktion der Zeit ist. Fir 
kleine Zeiten geht h/l gegen 0,28. Fir groBe Zeiten kann h/l erheblich von 0,28 abweichen. 
Diese Abweichungen sind sicher zum Teil auf die scheinbare Verbreiterung nach Hooker 
zuriickzufihren. Fage und Johansen’ geben h/l-Werte von 0,4 bis 0,5 an und stellen gleich- 
zeitig fest, daB die Wirbelkerndurchmesser in einem Abstand von 8 Kérperdurchmessern 
hinter dem Kérper von der Gré®enordnung der StraBenbreite waren. 

Die Arbeiten von Hooker, Wille und Fage und Johansen legen es nahe, in der Stabilitits- 
theorie der WirbelstraBen auf das reale Geschwindigkeitsprofil der Wirbel, das sich unter der 
Wirkung der Reibung ausbildet, Riicksicht zu nehmen. Wenn man im ibrigen das v. Kérmén- 


* Gekiirzte Fassung einer von der Fakuitat fiir Maschinenwesen an der Technischen Universitat 
en ne genehmigten Dissertation (Berichter: Prof. Dr.-Ing. R. Wille, Prof. Dr. phil. 
‘, Mohr). 
i a pene Proc. Roy. Soc. A. 79 (1907), S. 262. 
* H. Bénard, C. R. 147 (1908), S. 839 u. 970; 156 (1913), S. 5) 0 pias, isp y 
1523; 183 (1926), $. 20 u. 184 ee a eee ea 
* Th. v. Kérmén, Géttinger Nachrichten 1911 u. 1912. 
Bl. Dolaptschiew, Z. angew. Math. Mech. 17 (1937), S. 313; 18 (1938) an 200s 
Th. v. Karman und H. Rubach, Physikal. Z. 13 (1912), $. 49. 
S. G. Hooker, Proc. Roy. Soc. A. 154 (1936), S. 67. 
R, Wille, Jahrbuch S. T. G. 46 (1952), S. 176. 
A. Fage und F. C, Johansen, Proc. Roy. Soc. A. 116 (1927), S. 170. 
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sche Modell der WirbelstraBe beibehalt, dann ist, falls unter den abgewandelten Bedingungen 
noch Stabilitat miglich ist, der Quotient h/l als eine Funktion der Zeit zu erwarten, Diese 
Funktion mu8 die Bedingung erfiillen, dafB eingeleitete kleine Stérungen der Wirbel nicht 
bestandig anwachsen. 


2. Das Geschwindigkeitsprofil des realen Wirbels. Nach Lamb! gilt fiir das Abklingen eines 
Wirbels in der ebenen realen Strémung mit der kinematischen Zahigkeit », der Zeit t, der 
Zirkulation zur Zeit Null I, dem Radiusvektor r und der azimutalen Komponente der Ge- 
schwindigkeit c, die Beziehung 
a Lig a: 
= fal, =e : (1) 


Fir geniigend grofe r weicht cy bei konstanten t beliebig wenig vom idealen Profil c, = 1/2ar 
ab. Fiir Kleine r gilt die Entwicklung 
Lo 
Were ) 


Das bedeutet eine lineare Zunahme von c, mit dem Radius. Der Kern des Wirbels rotiert 
also annahernd wie ein starrer Kérper. 

Fir die folgenden Stabilitatsuntersuchungen soll das Geschwindigkeitsprofil (1) verwendet 
werden. Dabei wird es erforderlich, die Strémungsgeschwindigkeit in den WirbelstraBen durch 
Summierung der induzierten Geschwindigkeiten der Einzelwirbel zu berechnen. Diese Uber- 
lagerung ist streng wegen der quadratischen Glieder der Navier-Stokes-Gleichungen nicht 
erlaubt. Diese verschwinden fiir den isolierten Wirbel in der unendlich ausgedehnten Flissig- 
keit, sie verschwinden nicht mehr, wenn die Kreissymmetrie gestért ist. Fiir kleine Ver- 
haltnisse von Kernradius zu StraSenteilung wird man das Wirbelkerngebiet niherungsweise 
als isoliert in der unendlich ausgedehnten Flissigkeit betrachten dirfen. In der Umgebung 
der Wirbelkerne mu man sich auf die mit 1/r abnehmende Geschwindigkeit und die damit 
verbundene Vernachlassigbarkeit der quadratischen Glieder berufen. Man darf deshalb 
Naherungsaussagen fiir kleine Zeiten erwarten. Fiir grofBe Zeiten ist die Superposition nur 
dann erlaubt, wenn die quadratischen Glieder in den Navier-Stokes-Gleichungen klein gegen 
die Reibungsglieder werden, d.h. bei kleinen Strémungsgeschwindigkeiten. 


3. Die Stabilitatsuntersuchung der einreihigen Wirbelstrafe. a) Die Aufstellung der 
Bewegungsgleichungen. Man denke sich die WirbelstraBe zur Zeit t = 0 in der realen 
Flissigkeit entstanden. Die Wirbel seien nach Abb.1 auf der «-Achse eines kartesischen 
Koordinatensystems angeordnet. Fir die kartesischen Geschwindigkeitskomponenten wu, v, 
die ein Wirbel, der im Ursprung des Systems angeordnet ist, induziert, gilt im Aufpunkt (x, y) 


Tie Y 
ele ie Aries) ay 
u (x, y) i ) (; e js ? ve 
: aa) -0}- 5 Fa >i = 
tone Bho 8): | 
: 2a7 } Abb. 1. 


Die Einzelwirbel seien um kleine Strecken aus ihrer Normallage verriickt. Wenn die Stérungs- 
komponenten mit x, und y, bezeichnet werden, wobei fir n = — o~,..., — 1,0, + 1,..., + 00 
ist, so liegen die n Wirbel in den Punkten (nl + x,, y,). Legt man den Aufpunkt in einen 
der betrachteten Wirbel, etwa denjenigen mit dem Index x, so erhalt man als den Geschwin- 
digkeitsbeitrag aller Wirbel auf den Aufwirbel 


- ; Tun 
d Xx 14 > Yu — Yn 1 e 498 
= 5 ; , : 
cai nt ran 


dt 

d Toe = (n oi = (4) 
Si ny Ce a aed as UU \, 4vt 

die eee es re : : 
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1 H, Lamb, Lehrbuch der Hydrodynamik, 2. Auflage, Leipzig 1931, S. 669. 
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Dabei ist die Summierung iiber alle n ausgenommen n = x durchzufiihren, was durch den 
Strich am Summenzeichen angedeutet wird. Da x alle ganzen Zahlen durchlauft, liegt ein 
System von unendlich vielen Differentialgleichungen vor. 

b) Die Linearisierung des Differentialgleichungssystems. Die Stérungen 
(X%n»¥n) mobgen so klein gegen die Teilung / sein, daB eine Linearisierung des vorliegenden 
Systems (4) erlaubt ist. Die linearisierten Bewegungsgleichungen lauten 


(n—xz)P 
dxy pe al Yu —= Yn € < Ayvt 
== a es a - ee)? BE \ —e > 


dt 2G (= 23) 


= (n=l aes _(@—xP 
dyx ry Xn — Hn ( -* aot np Wiest alge Baa | 


dt Dh Gs SS r? 


(5) 


Bei der Linearisierung tritt die Summe 


; elie LE: 
> 7 ii al mete 4vt (6) 


n—x)1 


auf. Sie ist im Cauchyschen Sinne konvergent und hat den Wert Null. Physikalisch bedeutet 
dies, daB man sich bei der Betrachtung auf die mittleren Teile einer beiderseits unendlich 
langen StraBe beschrankt. 

Das System (5) von unendlich vielen Differentialgleichungen wird nach dem Vorgange 
von Th.v. Kérmdén' durch Annahme von Phasenbeziehungen zwischen den Stérungen der 
einzelnen Wirbel auf ein System von zwei simultanen Differentialgleichungen zurickgefiihrt. 
Man setat 

Re SHO OTRO) wy, ge vere? (CSO a2 az)e (7) 
Dann erhalt man, wenn die dimensionslose Zeit t =(4y/I*)t eingefibrt wird — sie bedeutet 
physikalisch bis auf einem konstanten Faktor das Quadrat des Quotienten aus Kernradius 
und StraBenteilung — und f(t), g(t) nach (8) definiert werden, 


I) ya2 SD aeeeels <4), 


ate) = 2 eel a gm aa | 
= 


fi) = Flat S)G) 0) 


G8)-G) mera 


c) Die Integration des simultanen Systems. Durch Integration von (9) unter 
Beriicksichtigung der Anfangsbedingung gewinnt man 


) (i) | atey 0") (ie) 2 % 


und indem man fiir Gey) (10) einsetzt, gelangt man zu 


B(t) 


i) = Be | (een 0") Ui) * la) i] ee ica ree) eats mm 


00 


das zu lésende System 


mit der Anfangshedingung 


1 Vel. FuBnote 4, von S. 400. 
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und schlieBlich durch wiederholte Fortsetzung des Verfahrens zu 


pen) is | (eyo) (Rt) + (ey | ae, (p) drt 


14% é ae 
; [om | | | eve 0 a acne t= 


0 


Wenn die Reihe (12) absolut konvergiert und der Differentialgleichung (9) geniigt, so ist sie 
die verlangte Lisung des Problems, da die Anfangsbedingung sicher erfillt ist. 
Die Reihe (12) konvergiert so gut wie die Exponentialfunktion. Angenommen, die Reihen 
fiir f(r), g(t) konvergieren, dann geniigen f(r), g(t) in einem endlichen Intervall 0 <1 <1, 
_ der Bedingung 
fOI<R, la)|<R, RDO. 
Das n-te Glied der Reihe (12) 1a8t sich nunmehr abschatzen: 


| CP te ot eee 


0, 


1 | X [enor 
| ( Po) een ee) 
Damit ist die Konvergenz bereits bewiesen. 
Es bleibt zu zeigen, daB® (12) der Differentialgleichung (9) geniigt. Durch Differentiieren 
von (12) erhalt man 


BLOG) +N] OH aaa) 


= Flay HR) aay) Get 


Die letzte Zeile zeigt, daB die Lésung (12) der Gleichung (9) geniigt. Da die Anfangsbedingung 
ebenfalls erfillt ist, ist die Reihe (12) die vollstandige Lésung des Systems (9)!. 

d) Stabilitatsdiskussion der einreihigen WirbelstraBe. Die & 
Lésung (12) wird je nach Wahl der Integrationskonstanten von wesent- 
lich verschiedenem Charakter. Tragt man die Anfangswerte der Sté- 
rungen nach Abb.2 in einem rechtwinkligen Achsensystem auf und oe 
wahlt man solche Anfangswerte, die durch Punkte des ersten und dritten 
Quadranten dargestellt werden, dann sind die Reihenglieder von (12) an W 
von alternierendem Vorzeichen. Im ibrigen unterscheiden sich Liésungen, 
die aus am Ursprung des 0, fy-Systems gespiegelten Anfangswerten her- 
vorgehen, nur durch das Vorzeichen. 

Wahlt man Anfangswerte der Stérungen, die durch Punkte des zweiten baw. vierten 
Quadranten dargestellt sind, so gelangt man zu Reihen (12), deren Glieder von einerlei Vor- 
zeichen sind. 


Abb, 2. 


Es ergibt sich zum Beispiel fiir die Anfangswerte «(0) = —|ao|, B(0) == |(B5)| 
eee po T ae } 
ae = | F(t) |Bol ars + (5) | | f (cs) g(t) [00 ; dt, + 
: les) || [4 8(t2) f(T) |Bo| dy dt, dt + °° | ? | aa) 


1 Eine anschauliche Deutung dieses Verfahrens findet man bei FE. Mohr im Zusammenhang mit 
der Kreiselbewegung in Math. Nachr. 6 (1950/51), S.1. 
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Lt Rt 


B =4|Bo| ame = i (v1) || de, + (eal | [ (ea) £0) [Bo | dy dt + 


00 
Ty 


: eal i Jao (t1) f (te) @(tg) ao | dr, dt. dt, +--+: ; 


Wiederum gilt, daB die ee (14) bei Spiegelung der Anfangswerte am Ursprung des a), [i,- 
Systems nur das Vorzeichen vertauscht. 
In (14) liegt eine Lisung vor, die, falls f(r), g(r) gréBer als Null sind, zu bestandig wachsen- 
den Stérungen fihrt. Es geniigt zum Nachweis der Instabilitat der einreihigen WirbelstraBe, 
, da8 mindestens eine Lisung zu 
pg bestandig wachsenden Stérungen 
a ae ee fiihrt. Man mu8 deshalb die Funk- 
tionen f(r), g(v) (8) naher betrach- 
ten. Die Konvergenz der Reihen, 
durch die f(r), g(z) dargestellt wer- 
den, ist sofort ersichtlich, da man 
die vorhandenen Reihen in solche 
aufspalten kann, die entweder 


co 


; : 1 ‘ 
mindestens wie > —, konvergie- 
mW 
T 


| 
=| Sie z ie eae ee 
0 Gi We Os 10 12 14 16 18 20 ren oder die wie > e 7 kon- 
Abb. 3. Zur Stabilitétsdiskussion der einreihigen WirbelstraBe. Verlauf ‘ ; 1 , 
der Funktionen f(z), g(z) tiber 7. vergieren. Die Funktionen f(t), 


g(r) hangen vom Parameter , der 

Phasenbeziehung zwischen den Stérungen, ab. Es geniigt, zu zeigen, daf fiir einen er- 
laubten Phasenwinkel g mit 0<qm<2a die Funktionen f(t), g(z) gréBer als Null sind. 
Dabei sind fiir m die Grenzen 0 und 27 ausgeschlossen. Die Werte yp =0; 22 liefern im Sinne 
y der Theorie keine echten Stérungen, 


= S eS . . 5 da sie einer Verschiebung der Wirbel- 


straBe als Ganzes in der (x, y)-Ebene 


- entsprechen. 
Zs ce 5 . 5 | , Fir die zahlenmaBige Auswertung 
i Hees, von f(t), g(t) wurde p=a gewahlt. 
In Abb.3 ist der Verlauf von f(r; 2), 
g(t;a) tber 7 dargestellt. 
Aus den dargelegten Ausfiihrungen und dem Verlauf von f(t; 2), g(r;2) ist bewiesen, daB 
die einreihige WirbelstraBe auch bei Bericksichtigung endlicher und zeitlich wachsender 


Wirbelkerndurchmesser instabil ist, da sie jedenfalls fiir den Phasenwinkel gy =z instabil 
ist. Damit ist der Anschlu8 der Theorie an die experimentelle Erfahrung weiterhin gesichert. 


Abb. 4. 


4. Die Stabilitatsuntersuchungen der zweireihigen WirbelstraBen. a) Die Aufstellung der 
Bewegungsgleichungen. 

Wiederum denke man sich die WirbelstraBe zur Zeit t =0 in der realen Flissigkeit ent- 
standen. Die Anordnung der Wirbel in der WirbelstraBe ergibt sich aus der Abb.4. 

Wenn man die Komponenten der Stérungen der Einzelwirbel in der oberen StraBenseite 
Mit Xm, Ym, diejenigen der unteren Strafenseite mit x), y, bezeichnet, erhalt man die Koor- 
dinaten der Einzelwirbel 


fiir die obere Reihe zu (mI => "Gm ts =a Yn) ° 


fiir die untere Reihe zu (nt ao = so Yn] 


und fiir den Betrag des Radiusvektors zum Aufwirbel (*. us _ ¥0) ergibt sich 


in der oberen Reihe in V(x ml Xm)? + (Yo — ¥m)* » 
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in der unteren Reihe Ton = | (onl Pe Eh ey, — 7). 


it I | 


Mit u,(t) als gemeinsamer Translationsgeschwindigkeit des Wirbelsystems in Richtung der 
Parallelen ergeben sich die Bewegungsgleichungen 


2% r2 22 2 
Om " (15) 
1 Le ah ee 4 n Sarnast 1 Th 
Stee ene ast A ee ee 
dt 2% m hes 2% n Ton 


Der erste Term der beiden Gleichungen stellt den Einflu8 der Wirbel der gleichen StrafBen- 
seite dar. Die Summation lauft iiber alle m, m—0 ausgeschlossen. Der zweite Term beider 
Gleichungen liefert den Beitrag der Wirbel der unteren StraBenseite. Die Summation lauft 
uber alle n. Die entgegengesetzten Vorzeichen der zweiten Terme gegeniiber den ersten 
Termen riihren von der dem Betrage nach gleichen, aber dem Vorzeichen nach entgegen- 
gesetzten Zirkulation her. 

b) Die Linearisierung des Differentialgleichungssystems. Wiederum sei an- 
genommen, daB die Stérungen (%m, Vm), (%,, ¥n) klein gegen die Teilung | seien. Die Lineari- 
sierung des Systems (15) fiithrt zu den Bewegungsgleichungen 


mo 
dx» Ve 1 ag 4yvt 
di 7 De am m2 12 (Yo Vn) 
272 at h? 
n? 12 + h? path i 
; lf n2 [2 — h2 ( ae > he é 4vt 1 
Og a Ge PL ep iw Re Saye Jo Jn 
Shik She mens P+ h? 
ibe > 2nlh ( ae | Bolle e om | , 
Sas (n? 1? + h?)? a n2 12 1 Ay t ame 16 
m? |? x m? |? ; ( ) 
er a meer, 
dy I to 4Ayt ‘ 
dt aa 25 Pei m2 [2 (X5 Rea) 
22 4 h2) n®]? + h? 
ee dh? ee ree 
r, n2 [2 — h?2 ( ear 2 he e | , 
nD et + h?)? oe M\2 P+ he 2 Tey i (%o Xn) 
212 4+ he n? [3 + h® 
f = 5 ange: 
Le > 2nlh ( Ave PAT IN | ; 
as 1 272 2 - n) e 
20 fe Pompey ° rape at (en 
Fiir die gemeinsame Translationsgeschwindigkeit der Wirbel u,(t) erhalt man 
n? 1? + he 
IB h ( “7 4 | 
: BS 1 17 
Ne gy an a PB? ; (17) 


Die Summe 


( n? [2 + i 
nl Te eye 
ty ay 1 18 
2 ape as 
die bei der Linearisierung auftritt, konvergiert im Sinne von Cauchy und hat den Wert Null. 
Wenn man in dem Gleichungssystem (16) den Summationsindex n durch n +- 1/2 ersetat, 
besitzt man die Bewegungsgleichungen fiir den Fall der auf Liicke stehenden WizbelstraBe. 
c) Vereinfachung der Differentialgleichungen fiir den Fall, daB nur der 
Aufwirbel gestoért wird. Es soll zunachst ein einfacheres Problem gelést werden. Die 
Storungen aller Wirbel bis auf den Aufwirbel mégen verschwinden. Dann gilt mit 


(Ese ia) = (0, 0) ? fi) SW ¢ 
(Xn ¥n) = (0, 0), 
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und unter Beriicksichtigung, daf diejenigen Summen, die in nm ungerade sind, verschwinden, 


das vereinfachte System 
Xp ae, ed Ly 0 F(t) Xo\ (19) 
yo) BNC) 0 } Aye) 


mit den Anfangsbedingungen 
(210) (8) tare =o. 
¥o(0) No 


F(z), G(r) sind folgendermaBen erklart: 


2.1 2 
m?* n + 


wae, He n? + k? a le | 
LOTS ees @ We = 12 = . 2 k2 T 
cea = es ae + elh : a a meer (20) 
mak a | 
ety) rt) +212 =e pea J 


d) Integration der Bewegungsgleichungen und Stabilitatsdiskussion der 
beiden Konfigurationsméglichkeiten der doppelreihigen WirbelstraBe. Die 
Lésung der Gleichung (19) ist nach Abschnitt 3, c) (13) bekannt. Man erhalt 


alien | (ote) 0°) Ga) at Gane. ota) Feed) (a) ae 


ioe L ey 
(EY LL (ogy ten cteg 95? PO) ante test | 


00 0 


Die Diskussion der Lésung ist fiir die beiden Konfigurationstypen der zweireihigen Wirbel- 
straBen gesondert durchzufiihren. Zunachst soll die auf Liicke stehende WirbelstraBe betrach- 
tet werden. Indem man in F(z), G(r) n durch n + 1/2 ersetzt, gelangt man zu F*(r), G*(r). 
Diese Funktionen lassen sich folgendermafSen darstellen: 


i? ee ee bial ee 
Pe =a *( 
(n+) +H 


Pea 
2S 2k? e G bing (22) 
"(m4 =p ke i 
penis P 
6*(2) = F(t) +2 [2 —_ ee 


Die Konvergenz der Reihen, durch die F*(r), G*(r) dargestellt werden, ist gesichert. Die 
Reihenglieder der Reihen in F'*(r) verschwinden wie 1/n?, die zusatzlichen Reihen in G*(z) 


ni 


aoe cakes da die Reihenglieder wie e * gegen Null gehen. Die Reihen besitzen reelle 
ummen, 


Die Funktionen F*(r), G*(t) hangen von dem Parameter k = h/l ab. Um eine Ubersicht 
tiber den Charakter der Lisungen (21) zu gewinnen, ist es zweckmaBig, nach den durch F*(r,k’) 
= 0 und G*(r, k’’) = 0 definierten Funktionen k’(r) und k’'(t) zu fragen. DaB®B solche Funk- 
lionen existieren, ist zu erwarten, da fiir 7 > 0 G*(r) > F*(r) geht und F*(0, k,) = 0 eine 
Lésung besitzt, namlich die alte Stabilitatsbedingung nach Th. v. Karman. Der 


; Verlauf von 
ki = k(t) und k’’ = k''(t) ist in Abb. 5 dargestellt. 
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Die beiden Funktionen unterteilen die (k, t)-Ebene in drei Teile. 
Oberhalb der Kurve k’ = k’(r) gilt F*(r, k) > 0 und Cree kes 0, 
unterhaib der Kurve k” = k'’(r) gilt F*(z, k) < 0 und G*(r, k) < 0, 
zwischen beiden Kurven ist F*(r,k) <0 und G*(r, k) > 0. 


Bei der Diskussion der Lésung (21) sind die drei geschilderten Méglichkeiten zu beriick- 
sichtigen. Angenommen F'*(r,k) > 0, G*(r,k) > 0; fiir diesen Fall geniigt es, cine Stérung, 
die durch einen Punkt des zweiten bzw. vierten Quadranten eines &,, NHo- Systems dargestellt 
ist, zu wahlen, um zu zeigen, 
daf} die WirbelstraBe fiir zu- 
gehérige (k,t)-Werte instabil 
ist. Der Nachweis kann hier fort- 
gelassen werden, da alles wie im O5 
Abschnitt 3, d) verlauft. 

Nimmt man F*(r,k) < 0, ue 
G*(t,k) <0 an, befindet man 
sich also unterhalb der Kurve ae 
k'’ =k’'(r), dann fiihren Anfangs- 
werte der Stérungen, die durch 02 
Punkte des ersten bzw. dritten 
Quadranten dargestellt werden, ae 


a 
os |" ———, ——-———- 


zu instabilen Lésungen. Es sei 


T 
x(0) =&>0, 0 
: Abb. 5. Stabilitats- und Instabilititsbereich der WirbelstraBe mit ,,auf Liicke“ 
¥(0) SSA e: 0 2 stehenden Wirbeln. 
dann gilt 
Tt DW 
Ae (Ty\? ff 
% = 8 +5 rai |F*(t)| 0 dt, + eee | |F*(z,)| |G*(v_)| & dry dr. +°-- 
0 0 


8a 


ae 
T oe 2 | (23) 
Yo = No + oa | |G*(t4)| 0 dey + (ex) | | |G*(t1)| |F*(t2)| mo dey dt, + °°: | 


Da die Summen nur positive Glieder besitzen, ist die Behauptung bewiesen. Interessant sind 
die’ Grenzfalle F*(7,k’) =0, G*(t,k5)> 0 und F*(ryk’) <0, C*(c,k’) =0. Auch sie 
stellen keine stabilen Anordnungen dar. Betrachtet man die allgemeine Lisung (21), so er- 
kennt man, daf im ersten Fall die Stérungskomponente x, konstant bleibt, jedoch die Sté- 
rungskomponente y, dem Betrage nach immer groBer wird. lm zweiten Fall ist es genau 
umgekehrt: Die Stérungskomponente y, bleibt konstant, wahrend die Komponente x) an- 
wachst. 

Die Kurven k’= k’(r) und k’’= k''(r) miinden fir kleine Zeiten mit gemeinsamer Tan- 
gente in die gestrichelt gezeichnete Stabilitatsbedingung der WirbelstraBe in der idealen 
Flissigkeit ein. (Die zahlenmaBige Differenz von 0,363 gegen 0,281 ergibt sich aus der Ver- 
einfachung, da® nur ein Wirbel gestért wurde.) Das setzt in Evidenz, dafi fir kleine Zeiten 
die Kdérmdnsche Stabilitatsbedingung in der Tat die einzig mégliche Stabilitatsbedingung ist. 
An dieser Stelle wachst weder die x)-Komponente noch die y)-Komponente an. 


Falls es unter den vorausgesetzten Bedingungen dieser Arbeit stabile Lésungen gibt, so 
sind sie nur noch im Gebiet F*(t, k) < 0, G*(t, k) > 0 méglich. Es sei 7 ein konstanter Fak- 
tor, dann fiihrt die Bedingung 


22 F*(c, ky!()) = — G*(x, ky'"(0)) (24) 


zu stabilen Lésungen. Man setze 


F*(r) =— A(t), G*(t) = P An); (25) 
28* 
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dann schreiben sich die Komponenten der Lésung (21) 


& Date 
pate 
Ne == Co Ho 3 =e | A(t;) dt, — & ieee i? al A(t) A(t,) dt, dt, 
) 0 0 
ety ty Ted a Ts 
7; bea | { [4e A(t) A(t) dt dt_ dt + & ie A | és | A(x,).. A(vq) dry. . dtq 
000 0 0 
r ig. I 4 t Py 
Yo =m — So? | Al) dey — np ter ze | A(t) A(ty) dry dry 
0 w 2 
ji rh eds ae t Ta 
Aten is | | A(t,) A(t) dt, dt3 +70 (a) is | - | A(t)... A(t) dt. . dt, 
0 0 0 0 


als doppelt alternierende Reihen. Man zerlegt sie in zwei je einfach alternierende Reihen und 


kann alsdann summieren. Die Zerlegung ergibt die Komponenten | 


hoa Egle 
Rene, p — 72 & [ A(t,) A(ts) dry dry + 28 Gea [. , | Ale): Alea) dey dry ie: a 
0 0 0 0 | 


Tons 
0 


1 pe ie r 
eae nf 8a | AG — 25 =a | : fA  A(ts) dty- dt, + — +-: | > 
0 


0 0 


ees ; = [4@) dr, — 23 eer | Aces) MAG yer ie ee | 
0 Oe 0 


Pie eres awe ye 
a h ae ie | | SaaS R a a iP * | A(t) ++ A(t4) dty--dty— 4+ —.- 
‘6 9 0 0 
Die Reihen (27) lassen sich summieren, wenn man beriicksichtigt, das die Beziehung 
% ™m—1 Tt 
. ra 1 mm 
| aaa | Aryans Ate dn Tipps dr, =n({ A(t,) | (28) 
3 6 0 


gilt. Die Reihen erkennt man als die Reihen des Sinus und Kosinus vom Argument 


t 
1 


a A f 

Age | A(t) dt, ; 
0 

und man schreibt (27) geschlossen in Matrizenform 


«| “| ft i A(t,) | = sin me i A(t) eal ih e 


(29) 


| —A sin ) et | A(t) dt, | cos I, Hy [ 4 dt, 
0 


ae 


oe Gesamtheit der durch (29) gegebenen stabilen Lésungen iiberdeckt den durch die Kurven 

= k'(z) und k’! = k (tr) abgetrennten Raum zwischen den Instabilitatsgebieten, wenn man 

Ain 0 <A < o variiert. Die Grenzfalle 27 > 0 und A—>oo sind bereits als F*(c k') =0 
, = 


und G*(z, k’’) = 0 diskutiert worden. Studiert man die sich 
so findet man fiir alle erlaubten 1 Ellipsen. emma crocs 
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Ks gilt 
2 2 
Xo Yo 
—+agti=1. 30 
a t ne Eo + 8 ee? 
50 T Fe lo ; 


Insbesondere ergibt sich fiir 7 =1 ein Kreis vom Radius r — (€ + 72)1/2. Die Stabilitat 
der Wirbelstrafe mit auf Liicke stehenden Wirbeln wird durch die geschilderte Gestalt der 
Bahnkurven anschaulich klar. Auch der Ubergang zur Instabilitat wird anschaulich, wenn 
man bedenkt, daB mit 1 > 0 und) — oo, also bei Annaherung an die Kurven k’’(r) und k’(z), 
jeweils die eine Achse der Ellipse immer gréSer und schlieBlich unendlich groBh wird, wahrend 
die andere Achse sich beliebig der gegebenen Anfangsstérung &) bzw. 7) nahert. 

Die Auswertung der Bedingung “ 
(24) ist fiir spezielle Werte von 72 aol? Re Sai, 7 
in Abb. 6 dargestellt. : | | 

Die Stabilitatsdiskussion der | 
WirbelstraBe mit gegeniiberliegen- 
den Wirbeln ist nach den Aus- 
fiihrungen iiber die WirbelstraBe 
mit auf Liicke stehenden Wirbeln 
schnell zu erledigen. 2A 

Die Konvergenz der Funk- 
tionen F(z), G(r) (20) ist ebenso 97 
gesichert wie diejenige von F*(r) 
und G*(r). F(z), G(r) hangeneben- 47 
falls vom Parameter k = hil ab. 
Es ist daher wiederum zweck- 
maBig, nach Lisungen von 
ick j= 0-sand G(r, hk’) = 0 
zu suchen. Die numerische Untersuchung zeigt, daB G(r, k) im untersuchten Bereich 0 <7 <1 
und 0,1< k <0,6 stets gréfer als Null bleibt. F(z, k’) = 0 definiert eine Lésungskurve, dic 
in Abb. 7 dargestellt ist. h 

Sie ist fir k— 0, nicht erklart. Ap u =. eee = 
k’ = k'(r) teilt die (k,r)-Ebene in 0<Fell 
zwei Teile auf. Im oberen Teil ist 
ea.ky> 0) “und —F(r, k) >-0; 
wahrend unterhalb der Kurve 


Ca.k)> 0 und Fi(r,k) <0 


ist. Nach den Ausfihrungen der 
vorausgegangenen Kapitel kénnen 
im Raum F(r,k) >0, G(r, k) > 0 
keine WirbelstraBen mit gegen- 
uberstehenden Wirbeln existieren, 
weil beliebige Anfangsstérungen 
schlieBlich bestandig anwachsen 
wiirden. Insbesondere sind solche Abb. 7. Stabilitats- und Instabilitatsbereich der WirbelstraBe mit gegenitiber- 
WirbelstraBen fiir kleine Zeiten, stehenden Wirbeln. ‘ 

also beiihrer Entstehung, instabil. 

Im Raume G(r,k) >0, F(t,k)>0 sind nach der hier dargestellten Theorie stabile Wirbel- 
straBen mit gegeniberstehenden Wirbeln méglich, allerdings werden solche WirbelstraBen 
diesen Raum durch die vorgelegte Instabilitatsbarriere nicht ohne weiteres erreichen kénnen. 


05; — 


OY 


| | cr 
0 G7 G2 G3 OY 45 06 O7 48 G9 a 


Abb. 6. Spezielle Lésungen der Wirbelstrafe mit ,,auf Liicke‘‘ stehenden Wirbelp 


| 


te 
0 O17 42 03 OF OS 06 G7 48 49 10 


5. Bemerkungen zu der vorgelegten Theorie. Eine quantitative Nachpriifung der vor- 
gelegten Theorie stoBt im Augenblick noch auf zwei Schwierigkeiten. Es liegen noch keine 
expliziten Ergebnisse fiir die Stabilitatsbedingung in dem Fall vor, daB die Stérung aller 
Wirbel beriicksichtigt ist. Diese Schwierigkeit diirfte jedoch bald behoben sein. Die ent- 
sprechenden Rechnungen sind in Arbeit und werden hoffentlich bald vorgelegt werden kénnen. 

Die zweite Schwierigkeit liegt darin, daB bei den bisherigen Versuchen wohl die Bedeutung 
des Kernradius angedeutet worden ist, da aber noch keine ausfiithrlichen MeBergebnisse iiber 


410 


Domm: Ein Beitrag zur Stabilitatstheorie der WirbelstraBen usw. Ingenieur-Archiv 


dieses Problem vorliegen. Im engen Zusammenhang damit steht die tatsachliche Zirkulation 
des Einzelwirbels in der WirbelstraBe. Diesbeziigliche Untersuchungen werden dariiber Aus- 
kunft geben miissen, ob zur Berechnung des Kernradius die kinematische Zahigkeit der 
Fliissigkeit ausreicht, oder, wie wir es vermuten, turbulente AustauschgréBen Beriicksichtigung 
werden finden miissen. Dennoch diirfte die Theorie auch in diesem Falle in der oben an- 
gegebenen Deutung der dimensionslosen Zeit Stabilitatsauskiinfte geben kénnen. 

Die Theorie liefert nach unserer Ansicht aber bereits jetzt neue Einblicke in das Stabili- 
tatsproblem der WirbelstraBen. Zur Erklarung der bei den bisherigen Versuchen gefundenen 
Nichtreproduzierbarkeit der MeBergebnisse bei der Untersuchung der WirbelstraBen mit auf 
Liicke stehenden Wirbeln wird durch die hier dargestellte Theorie die Méglichkeit aufgedeckt, 
das WirbelstraBen nicht fiir nur einen Wert von h/l stabil sein kénnen. 

Fir kleine Werte von rt baw. fiir verschwindende Zahigkeit geht die dargestellte Theorie 
in die Kdrmdnsche Theorie iiber, wie es auch verlangt werden muf. 


6. Zusammenfassung. Nach dem Bericht tiber die bekannten Arbeiten iiber die Stabilitats- 
theorie der WirbelstraBen wird die Vermutung ausgesprochen und begriindet, da® der in den 
bisherigen Theorien nicht beriicksichtigte endliche Wirbelkerndurchmesser und seine zeitliche 
Veradnderung einen neuen Einblick in das Stabilitatsproblem der WirbelstraBen verspricht. 

Es werden fiir die einreihige WirbelstraBe Bewegungsgleichungen aufgestellt, linearisiert 
und durch eine unendliche Reihe gelést. Die Stabilitatsdiskussion zeigt die Instabilitat der 
einreihigen WirbelstraBe in Ubereinstimmung mit der experimentellen Erfahrung. 

Im Anschlu®B hieran werden die Bewegungsgleichungen fiir die beiden Typen der zwei- 
reihigen WirbelstraBe angegeben, linearisiert und fiir den Fall, da®B nur ein Wirbel gestért 
wird, gelést. Die Diskussion der Lésungen ergibt, daB fiir die WirbelstraBe mit auf Liicke 
stehenden Wirbeln fiir endliche Zeiten ein Stabilitatsgebiet existiert. Es werden spezielle 
stabile Lisungen angegeben. Fiir verschwindende Zahigkeit schrumpft das Stabilitatsgebiet 
auf eine Gerade zusammen, die die alte Stabilitatsbedingung nach Th. v. Kérmdn darstellt. 

Die Diskussion der WirbelstraBe mit gegeniiberstehenden Wirbeln zeigt das iiberraschende 
Ergebnis, da neben der Instabilitat dieser Anordnung fiir kleine Zeiten bzw. fiir grobe Werte 
des Parameters k = h/l auch ein Stabilitatsbereich existiert, der jedoch nur durch Instabilitats- 
gebiete hindurch erreicht werden kann. 

Zum Schluf werden einige Bemerkungen iiber die Giltigkeit der Theorie und iiber die 
Méglichkeiten der experimentellen Verifizierung gemacht. 


(Eingegangen am 25, Januar 1954.) 


Anschrift des Verfassers: Dr.-Ing. U. Domm, Berlin-Charlottenburg, Technische Universitit, 
Hermann-F6éttinger-Institut. 
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Die Hohlkugel bei einem nichtlinearen Elastizitatsgesetz. 
Von F, Jindra. 


1. Einleitung. Im folgenden soll die Spannungsverteilung in einer homogenen isotropen 
dickwandigen Hohlkugel untersucht werden, die durch gleichmaBigen AuBen- und famen drivel 
beansprucht wird. Zur Beschreibung des elastischen Verhaltens des Materials beniitzen wir 
vein allgemeines nichtlineares Elastizitatsgesetz fiir kleine Verzerrungen.! In rechtwinkligen 
Koordinaten «,¥,# lauten die Spannungs-Dehnungsbeziehungen dieses nichtlinearen Elasti- 
zitatsgesetzes 


1 ] 
Puy = G8) tey> Pos = GB) tys, Yas = 78 (1B) Ten 


Hierin bedeuten ¢€,, &y, Es. Yay. Pyzs Pex die Verzerrungskomponenten, Gi Oyn Cal enn Vy eos 
die Spannungskomponenten, K den Kompressionsmodul und G den Schubmodul; s, und t, 
die dimensionslosen Spannungsgréfen 


ee ogc.) ee | 
(2) 


=F [ple — a +; — oo + (0. 09) ety te + eh). | 


Die Kompressionsfunktion k(s)) und die Schubfunktion g(i?) werden als Potenzreihen der 
SpannungsgréBen (2) angesetzt: 

k(s)) = 1 ky So t ky s?-+ hk, 58 + --- , | 

(i) =1ltegtestset---. f 
Der Verlauf der Funktionen (3) ist experimentell durch einfache Versuche zu ermitteln. 

Wir leiten hier zuerst allgemein die Differential- 

gleichung der Radialspannung her. Unter Beschrankung 
auf bestimmte einfache Formen der Funktionen (3) 
geben wir Lésungen in geschlossener Form und durch 
sukzessive Naherung an. Aus den Zahlenergebuissen 
kann man schlieBen, daB der in der Praxis auftretende 
Abbau der Spannungsspitzen durch kleine Abweichungen 
des Elastizitatsgesetzes von der Linearitat erklart werden 
kann. 


2. Die Differentialgleichung fiir die radiale Spannung. 
Wenn eine Hohlkugel nur durch gleichmafigen Auben- 
und Innendruck belastet wird, so mu der in ihr auf- 
tretende Spannungs- und Formanderungszustand sym- 
metrisch um den Kugelmittelpunkt sein. Das Krafte- 
gleichgewicht an irgendeinem Element der Kugel (Abb. 1) Abb. 1. Kugel unter gleichmaSigem Aufen- 
im Abstand r vom Mittelpunkt liefert die Beziehung EEN SU LS 
zwischen der Radialspannung o, und der Tangentialspannung o, 

r do, (4) 


ay alee Ta 


(3) 


Op 


1 H. Kauderer, Ing.-Arch. 17 (1949), 8.450; F. Jindra, Ing.-Arch, 22 CODA, Tio A 
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Mittels der radialen Verschiebung u lassen sich die Dehnungen in radialer und tangentialer 
Richtung darstellen durch 


ae se bo = (5) 
Der Zusammenhang der Dehnungen mit den Spannungen folgt aus dem nichtlinearen Elasti- 


zitatsgesetz (1). Mit den Hauptspannungen o,, 04, 0» bildet man zunachst gemaB (2) den 
Mittelwert 


1 
on = (6, + 2 4) : 


Dann wird aus (1) 


du Ae — 

& 7 = sR k(Sp) (co, + ey Oy) + 3 G b(t) (c, Oy) ° (6) 
u 1 1 9 

Ey = 7 = 7K Eso) (Or + 2 6) + GZ GSlh) (% — 9) 


mit 
9 


1 a 2 
Ty == OnKA (o, Se 2 Oy) ? is a 9 G3 (Or -) Oy)" ; 


Indem man mit der zweiten dieser Gleichungen du/dr bildet, diesen Ausdruck mit dem 
Wert der ersten Gleichung vergleicht und (4) einfiihrt, findet man fiir die radiale Spannung 
die Differentialgleichung 


seen + eta] (E+ 68g) + he a aa) Eye 
mit 
Sega oR" aa sas v1). = re ie (Z): (7a) 
Hierzu treten die Randbedingungen: 
(CG) 28 Pe PLU hare | 
(Go) Se pres fare er reel (8) 


wenn r, den inneren und r, den auferen Halbmesser bezeichnet. 

Die Lésung der Aufgabe hangt also von der Integration der Gleichung (7) bei vorgegebenen 
Randbedingungen ab. Wegen der auf erordentlich komplizierten Form dieser nichtlinearen 
Differentialgleichung sollen hier nur zwei Sonderfalle beachtet werden. 


3. Geschlossene Losung. Da sich die Kompressionsfunktion bei vielen Stoffen! als kon- 
stant erweist, so wollen wir zuerst den Sonderfall behandeln, daB die Funktionen (3) in der 
Form 


K(s)=1, g(t) =1lt+aey (9) 
darstellbar sind. Dann vereinfacht sich die DESOTEDEN Cine (7) der Radialspannung zu 
(a? r2 of? + 1) roy +2 (a? of? 4+ 2)o; =0 (10) 
mit der Stoffzahl 
oe Suh eee ee Gee fan IS 
YFG SKE EG Gi Ke ae) oe) 


die hier allein die Abweichung vom Hookeschen Gesetz kennzeichnet. Striche bedeuten jeweils 
Ableitungen nach r. 

Da o, in (10) nicht explizite vorkommt, geht (10) fiir p = o, in eine Differentialgleichung 
erster Ordnung iiber, die nach der Multiplikatorenmethode? lésbar ist. Der integrierende 
Faktor lautet 

1 


TP (RTS ps) 


M = 
ie Vel. H. Kauderer, a. a. O. S. 468. 


: 2 Vee Kamke, Die eaeni sec eee Lésungsmethoden und Lésungen, 3. Aufl., Bd. LS 28e 
eIpZig 
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und somit 1aBt sich (10) in der Gestalt 
a2 r? p2 + 1 x E: 2 (or? p? + 2) 
p(a®r® p? + 3) P " r(atr pt + 3) 


als eine exakte Differentialgleichung schreiben. Die Gleichung (11) besitzt eine Stamm- 
funktion 


a) (11) 


Ep) = = Ings pin 5? p> 42 3)". (11a) 


Die Lésung von (11) ist daher F(r, p) = konst. Fihren wir wieder o/ ein, so entsteht 
Oro. 42 37 o. = 2B =.0, (12) 


wobei 2 B die Integrationskonstante bezeichnet. Mit Hilfe der Cardanischen Formel! liésen 
wir (12) nach o; auf. Die einzige reelle Wurzel ist 


Ih Wey ———__—__- —__- 3 ———— 
=, Alla B+ Vo? B+ ro + Va B— Vo? BP 18). 


Q 
~~ 


Mit einer weiteren Integrationskonstanten A ergibt sich dann 


7s 2 3 ee dr 
=A+ | (ig B+Ve B+ r+ Va Blo Bp) yn. (13) 
Zwecks Auswertung von (13) fiihren wir die Substitution ein: 
Las — 2 art ——s 
0=—(in B+let BP tla B_le BL) (14) 
oder 
Tw -- 3 rw —2B=0. (14a) 
Damit wird (13) in das leicht lésbare Integral 
6, =A — lee a 
a2 w + 3 
verwandelt. 
So erhalt man endlich die eae adi 
Or, = A—w+ are feo (15) 


a6 ia 


Durch Grenziibergang «— 0 geht (15) in die bekannte Lisung o; fiir die lineare Theorie 
itber? 2 
BS Te (16) 
Zum Beweis muB man das letzte Glied in (15) in eine Potenzreihe entwickeln und den Aus- 
druck (14a) beachten. , 
Zur Bestimmung der Konstanten A und B haben wir die vorgeschriebenen Randbedin- 
gungen. Mit Hilfe von (15) erhalten wir aus (8) die Beziehungen 


ope wy + pee arc tg - otk, = Jie | 
x 3 3 (17) 
A — ww +— : are tg ot = — pas | 
a A /: 3 / 3 
wenn wir aa Pe 
= 3 . oa Fi ae Ee 
oe Be a B x? BE 
am =P +9 Pais |e Vat} 
Vo ieee LG I+ C@ 
und 


uw —/CR +1 LOR +VCR 1+ CR 


!Vel.z. B. Hiitte, 27. Aufl., Bd. I, S. 77, Berlin 1942. 
2 Vel. ZS: Picence Theory oe Elasticity, 5. 325, New York 1934. 
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mit den Abkirzungen 


Cpe ‘Rees (18) 


setzen. Aus (17) gewinnen wir dann durch Substraktion die Bestimmungsgleichung fiir die 
Integrationskonstante C bzw. B 


2 & (wy — wa) V3 


& (Po — Pa) = % (W) — Wa) — Fs are tg (19) 


3 + a? wo Wy 
Sobald die Lésung dieser in C transzendenten Gleichung bekannt ist, kann die Integrations- 
konstante A aus (17) leicht ermittelt werden. 

Die Tangentialspannung o, berechnet sich gemaf (4) mit (15) zu 


it 
Oo, =A — WwW 


are tg oe (20) 
und hat insbesondere am Innen- und Auenrand die Werte 
if i 
(Garr, SS + 2 Wo > (Ge \rars aa | i iy Wa: (21) 


Die tangentiale Spannung besitzt einen Extremwert 


1 iA il 
= (kg oe 
(Cg)r=r,, as A I Xx ( 5 /3 >] 
an der Stelle 
Le ee a B 
m 2 9 


welche natirlich nicht notwendig in den Bereich r, 
bis r, fallen muf. 
Der Vollstandigkeit halber sei hier noch die Tan- 


gentialspannung oj, fiir die lineare Theorie (x = 0) 
hinzugefigt 
0, =A = sd (22) 


Wir behandeln jetzt zwei Anwendungsfalle: 

a) Es sei py = py. In diesem Fall mu8 von vorn- 
herein die Konstante B gleich Null gesetzt werden, 
wie man sich aus (19) iiberzeugen kann. Man erhalt 
dann an jeder Stelle auch fiir das nichtlineare Elasti- 
zitatsgesetz 


Abb, 2. Verlauf der Radial- und Tangentialspannung. 


OK 1 Lint 
b) Es sei p, = 0, aber p) 4 0. Dann folgt die Integrationskonstante C gemaf (19) aus 
(Ly ee me X (Wy — Wa) V3 
a (w, — w,) — 3 are Daren rae = 


Fur die Spannungen ergibt sich 


2 a > ——— 9g) “3 
Orne (w Wa) ——— arc ie a (10 ea) 3) 6 
al 3 3 + 0? wy w 
CC. = i (w hag ) se 2 aioe a (w —- wy) 3 (23) 
— ole is 4, oem 


und somit fiir die Randwerte 


1 I 
(os) z= To = Po 9 Wo ? Cop ee aad 5 Wa S (23a) 
In Abb. 2 ist der Verlauf der Spannungskomponenten co, und o, langs eines Halbmessers 
fiir py = 500 kg/cm?, x = 10-3 cm?/kg und R =2 dargestellt. Man erkennt den starken Abbau 
der Spannungspitze (o,),_,, infolge der Nichtlinearitat. 
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Um einen Uberblick iiber die Abweichungen der Lisungen (23) gegenitber den Werten 
der linearen Theorie (22) zu erhalten, wurden die Randwerte (23a) der Tangentialspannung 


fiir verschiedene Werte der Stoffzahl « und R — 1,25 berechnet. In Abb. 3 und 4 sind die 
Anderungen 


Adm  %—%, 
ea) ek 
a) ae) 
Ap 1Op 
Op ox 
~ 04 ; 5 o4}-© —);——— - 
—O3 03 - =I a 
0=8-10 
G2 OC = 
: o 
| a\ 
neg 
tt 
“G7 07 yee 
=i 
[0 a= [to 
0 700 £00 JOO 400: kg/em* 500 0 100 £00 300 400 kg/em? 500 
Abb. 3. Anderung der Tangentialspannung am Innenrand Abb, 4. Anderung der Tangentialspannung am Aufenrand 
(fir R = 1,25). (fiir R = 1,25). 


am Innen- und AuBenrand iiber dem Innendruck p,faufgetragen. Wie zu erwarten, sind die 
Anderungen der Tangentialspannung am Innenrand negativ und am Aufenrand positiv. Die 
bezogenen Anderungen werden um so gréfer sein, je gréBer die Stoffzahl x ist, die die Ab- 
weichung des Elastizitatsgesetzes von der Linearitat angibt. AuBerdem wurden die Ande- 
rungen Ao,/o; am Innen- und AuBenrand fir 


46; : 
me x = 10-* cm?/kg und verschiedene Annahmen 
06|* von R ermittelt. Die Ergebnisse sind in Abb. 5 
| und 6 dargestellt. Wie man aus Abb. 5 ersieht, 
| wird die Abnahme der Tangentialspannung an 
~05 : | der Stelle r = r, fir groBe Werte von p, um so 
‘ gréBer sein, je gréBer das Verhaltnis der Halb- 
SY messer ist, d. h. je gréBer die auftretende Span- 
eee » nungsspitze am Innenrand ist. 
A0p 
| OF 
-03 1 a3” 
Ra105 
-02 a8 
R=11 
=U i 
| 
| 
| [to [tg 
— Z + 4 7 = 
0 700 200 300 400 xglon’ 500 0 700 200 300 400 Xgjone 500 
Abb. 5. Anderung der Tangentialspannung am Innenrand Abb, 6. Anderung derjTangentialspannung am AufSenrand 
(fiir « = 10-8 cm?/kg). (fir « = 10—% cm?/kg). 


4. Lésung durch sukzessive Niherung. Jetzt wollen wir uns dem praktisch wichtigen 
Sonderfall zuwenden, da® die Reihenentwicklungen (3) wenigstens nach den quadratischen 
Glicdern abgebrochen werden kénnen, also die Kompressions- und Schubfunktion in der 
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Form 


k(sy) = 1 + hy 59 + he sf , (4) = 14 gt (24) 
darstellbar sind. Aus (7) wird dann 


[a2 r® of? + By (rot + 30,)? + By (ro, + 30,) +1] roy + 2 [o? ro? 


+ 2 By (ro, + 30,)? + 2B, (ro; + 34,) + 2]o,=0, (25) 
wenn folgende Stoffzahlen eingefiithrt werden: 
s fe), BREE AG Bo K 
Se Ga K © 2OGR. 46)’ 
2h 3K 446 _ 8k, G 
box a6. ~TROK S4CC (25a) 
ye ee 
Po =o Kp? Wee CON.G ew ar 


Striche bedeuten wiederum Ableitungen nach r. 

Die nichtlineare Gleichung (25) lésen wir durch sukzessive Approximation. Durch Division 

mit dem ersten Klammerausdruck bringen wir (25) auf die Form 
PAR a elke Psd 


a® r’o!? + 8, (ro! + 3,)? + B, (ro! + 36,) +1 ~ 


ro, +40, = (26) 


Wir integrieren zuerst die homogene Gleichung, indem wir die rechte Seite von (26) gleich 
Null setzen. So erhalten wir die nullte Naherung, die bekannte Lésung der linearen Theorie 


B B 
Die Randbedingungen (8) liefern dann bekanntlich 
Pa v2 — Po Ts ii ike 
A= 3 ae: ef B= 3 se (Pa — Po) + (27a) 


Den Lésungswert (27) fihren wir jetzt in die rechte Seite von (26) ein und gewinnen somit 
fiir die erste Naherung die Gleichung 


(1)r7 Gby? 54 a? B 
NOs Oe (A275 + 9 a2 B2) r4 (28) 
mit der Abkirzung 
A®@ = I 3 py A498, AAs (28a) 
Nun hat aber die nichthomogene Differentialgleichung 
ay Phy f(x) 
die partikulare Lésung 
il 1 1 i 
yay [Se) de — 575 | ASle) de =a | | fle) de dia 
Somit erhalt man die erste Naherung der radialen Spannung o™ in der Form 
(ay D B ; 9 a? B? DAG Aor, 
Os OSes pam (4" r r& | on Sonn (29) 


und gemaf (4) die zugehérige tangentiale Spannung 


(2) eee | \ 1 ( 7) 9 o? B? 2A Aor 
Op EKG 
; ne In{ A’? 4 + Oa atc tg az: (30) 


2r3 ' 78 


r6 


D 
(1 


x 
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C und D sind Integrationskonstanten, deren Bestimmung aus den Randbedingungen (8) 
zu erfolgen hat. Unter Beachtung von (27a) findet man 


joe IEE 
Pes hice 2A’ 1’ 73 4/73 
C= A(rz — r8 Wie enti aly gus 3 se 
eat) pee) Oa aa BOOT ge) |: 
: 5 
i 


G 
Ta 


B 3 19 9 x2 B ies ( 2 2 
ee eg ent eo) ana 2s 


gatas Pacts AGE Meet) il 
TG gy 1810 BNO NB A ere + Oo BAYS 


Die Anderungen der Spannungskomponenten infolge der nichtlinearen Zusatzglieder in (25) 
berechnen sich dann zu 


Ao, =o — 6, Ao, = oy = o'®) ; (31) 
Als Anwendungsbeispiel soll wieder der Fall behandelt werden, da8 die Hohlkugel nur 


durch den Innendruck p, belastet wird. Die Naherungslésungen nullter Ordnung nehmen die 
Werte an 


mit 


wobei wieder zur Abkirzung R = r,/r, gesetzt ist. Die Anderungen der Spannungskompo- 
nenten werden dann nach (29) bis (31) 


I DYE Lae igs ro\ 2VAG Ao? 
/ == Liga? OT 0 ‘a 12 2 Dus 2 eee ee Ae Rea 
Ao, = Ch + 24 = sin(4 ee 
5 Ds Aare weal is Siete 2A“ A’ 7 
Ag, = Co ae (1 a In (4 + 9a An eg ue rey ee. 
mit 
Ay 51, 424902 ABR? | 2A’ fo, = A’ 
po = 3 (RP) IR Pires soda eee, ee 
und 
3 
yh 3G i In (A'2 + 9 92 A? RS) — In (A’2 + 9 a? 42) 
2 A! peer ge 
ey 8g oigt eRe |S 


woraus die Randwerte fiir die Anderung der Tangentialspannung folgen: 


A, R? A? +. 902 AB RE | A’ R? 3.0 4’ Ay (R? —1) 


= ~— A, R?° 
(A¢,),-r, = o7Rs— 1) @ a2 4908 BT Ba(Ri— 1) B42 + 908 BR! Ureree 
Ay R3 , A? + 902 A? RE oa 34.4’ A (R'—1) 
(Aoe)-=r, = 3 (RI—-1) "A? + 90? B 3 (Re) oe eo AR ee 
Setzt man in diesen Audsriicken A’ = 1 (d. h. B, = B. = 0), so erhalt man eine Naherungs- 


lésung fiir den Sonderfall (9). Wie die Durchrechnung einiger Zahlenbeispiele zeigte, ergeben 
die hiermit erhaltenen Werte eine gute Ubereinstimmung mit den exakten Liésungen von 


Ziff. 3. 
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5. Zusammenfassung. Mit einem allgemeinen nichtlinearen Elastizitatsgesetz fir kleine 
Verzerrungen werden die Spannungen in einer dickwandigen Hohlkugel berechnet. In allen 
Zahlenbeispielen zeigt es sich, daB die infolge der Nichtlinearitaét des Elastizitatsgesetzes 
errechnete Spannungsverteilung erhebliche Abweichungen gegeniiber den Werten der linearen 
Theorie ergibt, und zwar insbesondere fiir die Randwerte. Die Ergebnisse legen somit die 
Richtigkeit der Vermutung nahe, wonach bereits kleine, bei den meisten Werkstoffen vor- 
kommende Abweichungen vom Hookeschen Gesetz einen starken Abbau der Spannungs- 
spitzen zur Folge haben. 

Zum SchluB méchte ich Herrn Prof. Dr. R. Grammel sowie der Deutschen Forschungs- 
gemeinschaft meinen aufrichtigsten Dank fiir die Unterstiitzung bei der Durchfiithrung dieser 
Arbeit aussprechen. 


(Eingegangen am 5. Februar 1954.) 


Anschrift des Verfassers: Dr.-{ng. F. Jindra, Stuttgart-O, Gerolstr. 70. 
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‘Einftthrung in die theoretische Gasdynamik. yon pz, rover 


Sauer, o. Professor fiir Mathematik und analytische Mechanik an der Technischen Hoch- 
sechule Miinchen, Zweite Auflage. Mit 107 Abbildungen. VIII, 174 Seiten Gr.-8° 
1951. . : : DM 16,50 


»:++ Das Buch gibt eine ausgezeichnete und gut verstandliche Einfdhron in die Theorie der kompressi j 
3 Bez ; essiblen Stré- 

ine Bh einer kurzen Erliuterung der Grundbegriffe der Gasd Seestie werden zunichst die stationiiven Unteee 
ae ees caeih eke 9 in linearisierter Form behandelt. Dann folgt ein Abschnitt tiber Uberschallstrémungen 
eee ee gel und die Theorie des VerdichtungsstoBes. Ein weiterer Abschnitt iiber nichtlineare ebene und dreh- 
oh se rische Strémungen befaBt sich zunichst mit den verschiedenen Iterationsverfahren fiir kompressible Unter- 
at allstromungen. AnschlieBend werden die Transformationen von Molenbroek und Legendre zur Linearisierung der 
al ee Differentialgleichungen der Gasdynamik und das Prandtl-Busemannsche Charakteristikenverfahren fiir 
erschallstrémungen besprochen. Zum SchluB dieses Abschnittes ist in der neuen Auflage noch ein Paragraph ein- 


gefiigt, der die Erweiterung des Charakteristikenyerfahrens auf nichtwirbelfreie Strémungen bringt, wie sie insbesondere 


hinter VerdichtungsstéSen auftreten. Der neueren Entwicklung der G i i 
u BSSt6 treten. ‘ g der Gasdynamik folgend ist am Ende des Buches noch 
a ibe eek peasy itber_dreidimensionale Uberschallstrémungen eingefiigt.. Hier findet man die Theorie des Trag- 
ie ee oe wate shea das Singularititenverfahren wie die Theorie der kegelsymmetrischen 
2 I usemann besprochen. Auferdem werden in diesem Kapitel h di é 6 
Severin seam Cate em Kapitel noch die Strémungen an Drehkérpern 


Allen ‘denjenigen, die sich in die Theorie der kom i * ; ; 
. ‘ : pressiblen Gasstrémungen einarbeiten wollen, kann das Buch nur 
bestens empfohlen werden,**. "(Die " Naturwissenschaften) 
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Gasdynamik. Von Dr. Klaus Oswatitseh, Dozent an der Kénigl. Technischen Hoch- 
schule in Stockholm, friiherer wissenschaftlicher Mitarbeiter am Kaiser-Wilhelm-(Max- 
Planck-)Institut fiir Stromungsforschung in Géttingen. Mit 300 Textabbildungen und * 

3 Tafeln. XIII, 456 Seiten Gr.-8°. 1952, - . 


Inhaltstiibersicht: Thermodynamik. — Stationare und instationire Fadenstrémung. —. Allgemeine Glei- 


 chungen und Satze. — Spezielle Anwendungen der Integralsitze. — Allgemeine Gleichungen und spezielle, exakte’ 


Lésungen fir stationare reibungslose Strémung. — Stationire, reibungsfreie, ebene und achsensymmetrische Unter- 
schallstromung. — Stationiire, reibungsfreie, ebene und achsensymmetrische Uberschallstrémung. — Stationire, 
reibungsfreie schallnahe Str6mung. — Spezielle stationare und instationiire riumliche Strémuagen. — Strémungen 
mit Reibung. — Uherblick iiber die Versuchstechnik, Analogien. — Namen- und Sachverzeichnis. 


Die Tafeln sind auch einzeln zu beziehen 


: Bei Einzelbezug Bei Abnahme von 20 und mehr - 


Tafel I: DM 2,80 DM 2,40 
TafelII: - DM 2,80 ~ . DM 2,40 
Tafel III: DM 1,50 DM 1,30 


+>. Der Verfasser, z. Z. Dozent an der Kgl. Techn. Hochschule in Stockholm und bis zum Ende des zweiten Welt- 
krieges langjahriger Mitarbeiter von Prof. Prandtl am Kaiser-Wilhelm-Institut fiir Strémungsforschung in Gottingen, 
hat sich der dankenswerten Aufgabe unterzogen, das sehr umfangreiche, abex in vielen, zum Teil nur schwer zu- 
ginglichen, in- und auslandischen Zeitschriften und Abhandlungen niedergelegte Yorschungsmaterial der ilteren und 
besonders auch der neueren Zeit zu sichten und in einer zusammenfassenden Darstellung niederzulegen. Dabei wird 
das Hauptgewicht nicht zu sehr auf die mathematischen Fragen gelegt, die dabei eine Rolle spiclen, sondern vielmehr 
auf die klare Herausarbeitung der physikalischen Gesichtspunkte, die ihren Niederschlag in vielen technischen An- 
wéndungen finden .. ,“* , (Zeitschrift fiir angewandte Physik.) 


Grundlagen der Aeromechanik und Flugmechanik. __ 
"Von Dr.-Ing. Arthur Proll, Professor i. R. an der Technischen Hochschule Hannover. Mit 
278 Textabbildungen. XVI, 612-Seiten Gr.-8°. 1951. Gayzleinen DM 48,— 


Inhaltsibersicht: Grundlagen der Aerodynamik. Physikalische Grundlagen. — Mathe- 
matische Grundbezviehungen der Aerodynamik. Kinematik der Strémungen. Dynamik der 
Stromungen. Abbildungsaufgaben und Aerodynamik der Profile. Strémungen mit Widerstinden, Tragfliigel von 
endlicher Spannweite. Aus der Gasdynamik. — Flugmechanik. Theoretische und praktische Aufbauelemente. Vlug- 
leistungen. Die Geschwindigkeitsleistungen am Flugzeug. Steigleistungen eines Flugzeugs. Sonderleistungen von 
Flugzeugen. Nichtstationaére Flugzustinde. Flugeigenschaften. Stabilitit der Flugzeuge. Dynamische Stabilitit 
des Flugzeugs- —Namen- und Sachverzeichnis. “; 


,,..- Das Buch enthalt auBer den theoretischen Betrachtungen auch reichhaltiges Exfahrungsmaterial sowie eine 
ganze Reihe durchgerechneter praktischer Anwendungsbeispiele. Die Ausstattung ist bervorragend. Alle Studieren- 
den, Ingenieure und Physiker, die sich mit dem Problem des ,,Fliegens‘‘ beschiftigen wollen, werden aus ihm reichen 
Nutzen ziehen.‘* 3 (Zeitschrift fiir angewandte Physik.) 
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~ Technische Strémungslehre. ; 
_ besserte Auflage. Mit 407 Abb. X, 422 Seit 


Inhaltsibersicht: 1. Hydrostatik: Der statische D 
“ spire. — HM. Bewegungsiehre: Vergleich mit der Bewegung fester KO 
- senkrecht ‘zur Strémungsrichtung. Drehung von Fliissigkeitsteilchen. — 
der Flissigkeit als Kontinuum. Drehungsfreie Bewegung. — 
strémung (der Wirbel).} Zirkulation. Das “Potential. Energi betrachtun 
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- reibungsfreien Pliissigkeiten. Quellen und Senken. Gra 
‘Verfahren. -Besonderheiten der rotationssymmetrisehen Str 
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mung. Impulssatz. Reibungsfreie Strémung in offen 
- ‘$trémungen: Innere Reibung infolge Zahigkeit, Ahnli 
-_- Schmiermittelreibung. Strémungsformen. Strémungswide 
- bulenz. Druckverlust in ruhenden Schiittungen. Ro 
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Flissigkeit. 

_ (Aerostatik). —Dritter Abschni 
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